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1 Introduction 


1.1 Units 

We will mostly work with natural units where Planck’s constant divided by 2i r, 
the speed of light and Boltzmann’s constant are set to one, i.e. 

h = c = ks = 1 • (1.1) 

As a consequence, for instance, the mass of a particle (m) is equal to its rest energy 
(me 2 ) and also to its inverse Compton wavelength ( mc/h ). More generally, we 
have that the units of 

[energy] = [mass] = [temperature] = [length] -1 = [time] -1 

may all be expressed in terms of the basic unit of energy where we take the 
electronvoIt. (eV), i.e. the amount of energy gained by the charge of a single 
electron moved across an electric potential difference of one volt. Recall that 

GeV = 10 3 MeV = 10 6 KeV = 10 9 eV. 


The conversion to other units, such as length, can be obtained from h.c/MeV ~ 
197.33T0 -15 m = 197.33 fm with c = 2.9979-10 8 m/s and h = 6.5822-KT 22 MeVs. 
For later use we have for the conversion to length, time, temperature and mass: 


MeV -1 = 197.33 fm, 
MeV -1 = 6.5822-10 -22 s, 
MeV = 1.1605 ■ 10 10 K , 
MeV = 1.7827- 10 -2 V 


1.2 Isolated quantum systems in extreme conditions 

In recent years we have witnessed a dramatic convergence of research on 
nonequilibrium quantum systems in extreme conditions across traditional lines 
of specialization. Prominent examples include the evolution of the early universe 
shortly after a period of strongly accelerated expansion called inflation, the initial 
stages in collisions of ultrarelativistic nuclei at giant laboratory facilities, as well 
as table-top experiments with degenerate quantum gases far from equilibrium. 
Even though the typical energy scales of these systems vastly differ, they can 
show very similar dynamical properties. Certain characteristic numbers can even 
be quantitatively the same and one may use this universality to learn from table- 
top experiments with ultracold atoms something about the dynamics during the 
very early stages of our universe. 

Here the focus lies on isolated quantum systems, whose dynamics is governed 
by unitary time evolution. This sets them apart from other condensed matter 
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Figure 1: Schematic evolution towards thermal equilibrium as exemplified in 
sections u and 11.41 Starting far from equilibrium, the system is attracted 
towards a nonthermal fixed point, thus acquiring its universal properties before 
relaxing towards thermal equilibrium. (K. Boguslavski) 


systems where the influence of the environment is largely unavoidable. Isolated 
systems offer the possibility to study fundamental aspects of quantum statistical 
mechanics, such as nonequilibrium instabilities at early times and late-time 
thermalization from first principles. A quantum many-body system in thermal 
equilibrium is independent of its history in time and characterized by a few 
conserved quantities only. Therefore, any thermalization process starting from 
a nonequilibrium initial state requires an effective loss of details of the initial 
conditions at sufficiently long times. 

An effective partial memory loss of the initial state can already be observed 
at earlier characteristic stages of the nonequilibrium unitary time evolution. 
The corresponding prethermalization is characterized in terms of approximately 
conserved quantities, which evolve quasi-stationary even though the system is 
still far from equilibrium. An extreme case occurs if the nonequilibrium dynamics 
becomes self-similar. This amounts to an enormous reduction of the sensitivity 
to details of the underlying theory and initial conditions. The time evolution 
in this self-similar regime is described in terms of universal scaling exponents 
and scaling functions associated with nonthermal fixed, points , which is similar 
in spirit to the description of critical phenomena in thermal equilibrium. While 
the notion of thermal fixed points describing different universality classes is well 
established for systems in thermal equilibrium, the classification of far-from- 
equilibrium universality classes is a rapidly progressing research topic. 

In the following, we will outline these concepts for the examples of relativistic 


4 





inflaton dynamics in the early universe in section 11.31 and for the nonrelativistic 
dynamics of an ultracold Bose gas in section 11.41 These sections are meant to 
provide an overview of the different stages, starting far from equilibrium such that 
the system approaches a nonthermal fixed point before eventual thermalization. 
A schematic view is given in figure CD 

Introducing the basic ingredients for nonequilibrium quantum held theory 
in more detail in section [2j we analyze the notion of thermalization in isolated 
quantum systems using functional integral techniques in section [3] Section [4] 
explains to what extent nonequilibrium quantum held theory can be mapped onto 
a classical-statistical held theory problem, which provides powerful simulation 
techniques. Section [5] is devoted to nonequilibrium instabilities and the 
subsequent evolution towards nonthermal hxed points is explained in section |6j 

1.3 Heating the universe after inflation 

Inflation: It is one of the most striking developments in our understanding of 
the universe, that its properties at smallest length scales described by quantum 
held theory are tightly connected to its properties at largest scales up to the size 
of the observable universe today. In order to discuss some aspects of this, which 
will be relevant for these lectures, we recall that the universe is expanding with 
cosmological time t. The change is described by the Hubble parameter 

tt “(t) 
a(t ) ’ 

which is defined in terms of the scale factor a(t ) for a homogeneous and isotropic 
universe and the dot denotes the time derivative. The scale factor relates the 
proper distance between a pair of objects such that space itself is expanding 
with time. At the present moment to we take a(to) = 1, where roughly 
t 0 — 1/H 0 ~ 14 Gyr ~ 4.4 • 10 17 s. 

The expansion has not been the same at all times. Inflation denotes a 
period at very early times, where the expansion is strongly accelerated with 
a > 0. During this period a tiny region of space expands to a huge size. As a 
consequence, spatial curvature is decreased and the universe acquires a practically 
flat geometry. Furthermore, the accelerated expansion becomes so fast that 
points on the sky which are not in causal contact today were indeed correlated 
before inflation. This explains why we can see cosmic microwave radiation with 
practically the same temperature from all directions of the sky, thus providing a 
solution of this so-called horizon problem. Most strikingly, quantum fluctuations 
stretched to macroscopic sizes by the expansion can explain the origin of initial 
perturbations in the matter density. The gravitational instability of these initial 
perturbations then leads to the observed large-scale structure formation into 
clusters and superclusters of galaxies along with the fluctuations of the cosmic 
microwave background radiation, all in remarkable agreement with observations. 


( 1 . 2 ) 
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Knowledge about the amplitude of the perturbations can give important 
information about the energy scale of inflation. A typical characteristic energy 
density may be around e lat ion (10 15 — 10 16 GeV) 4 , which is reminiscent of 
the grand unification scale above which the electroweak and strong forces may 
become similar in strength. However, many different scenarios still exist and in 
the following we will concentrate on some generic aspects. The energy density at 
a given time, e, determines the Hubble parameter by the Friedmann equation 


H 2 = 


8nG 

-e 


(1.3) 


Here G denotes Newton’s gravitational constant and we define the reduced Planck 
mass M P = (8v tG )~ 1 ' 2 ~ 2.435 • 10 18 GeV. 

In a simplest model of inflation, the dominant contribution to the energy 
density at the time is described in terms of a spatially homogeneous scalar field 
0(f), which can be decomposed as 


e <t> ~ 2 ^ + (l- 4 ) 

into a kinetic energy part 0 2 /2 and a potential energy part K(0). Likewise, the 
pressure is given by the difference 

= (1.8) 

As a consequence, a time-independent held has negative pressure: P^ = —e^ with 
constant energy density. According to (1 1.3 f) . in this case H = const such that 

- = const , i.e. a(t) ~ e Ht . (1.6) 

a 

To be consistent with the observed properties of the cosmic microwave radiation, 
one may need about 60 e-folds of growth before the end of inflation. This requires 
the specification of a suitable potential K(0) along with initial conditions that 
allow for such dynamics. It should be stressed, however, that there is a large 
amount of models of inflation with different degrees of freedom and levels of so¬ 
phistication. Since this will not be relevant for our purposes, we continue for the 
moment with our simple single-held model. 

Coherent field evolution: Generically, the previous matter (with pressure 
P m = 0) and radiation (P r = e r /3 describing relativistic particles) content of the 
universe is diluted away during inflation according to 

Gn ~ -o , e r ~ -7 • (1.7) 
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The suppression factor for the radiation can be understood from the fact that 
its energy density e r ~ n r u) r is given by the number density scaling with inverse 
volume as n r ~ 1/a 3 times the mean energy per particle being red-shifted as 
ui r ~ 1/a, since the wavelength is stretched by the expansion. As a consequence of 
the dramatic dilution, after inflation all the energy is stored in the large coherent 
field amplitude of the inflaton. Inflation ends once the time-dependence of the 
initially slowly moving field 0(f) can no longer be neglected. In this case, there is 
no constant energy density dominating and the solution (II.6ft is no longer valid. 

To discuss the subsequent time evolution, we first note that the classical 
equation of motion for the scalar field is given by 

. dV 

<f> + 3 H(f> -|——— = 0 . ( 1 - 8 ) 

dcp 

This may be seen as the equation of motion for a relativistic scalar field with a 
“friction” term 3 Hcf> due to the expansion of the universe. For the following, we 
also specify to a quartic potential of the form 

V{<t>) = \ ™ 2 <f? + ^ 0 4 • (1-9) 

Furthermore, we note that for relativistic particles dominating the energy density 
the scale factor evolves as 

apt) ~ t 1//2 , H ~ —, (1-10) 

t 

in agreement with H 2 ~ e r ~ 1/a 4 according to (11.31) . In this case the equation 
of motion for the field can be rewritten in a form that is very useful for later 
analysis: We introduce conformal time t c and a field rescaling as 

/ d£ 

— , 0c = a 0 > (1-11) 

a 

where the integration boundaries run from some initial time to the final time 
(today). Then the field equation (11.81) for the potential (11.91) becomes 

0c + (mV - ^ 0c + ^0 3 = 0 , (1.12) 


where primes denote f c -derivatives. For radiation domination we have apt ) rsj t l/2 
and with (II ,11ft that t c ~ 21 1 / 2 . Consequently, a ~ f c leads to a" = 0. From (II. 12ft 
we observe that for massless fields, where m 2 = 0, the dynamics in expanding 
space is conformally equivalent to dynamics in Minkowski space-time with no 
expansion in this case, i.e. 



(1.13) 
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For a given non-zero initial field amplitude, this equation describes an oscillating 
field in a quartic potential. For notational simplicity, in the following we will 
drop index labels and employ 0(t) to denote the rescaled field in conformal time. 

Nonequilibrium instabilities and preheating: As a classical homogeneous 
field, the solution of (11.13j) would continue oscillating. However, in a quantum 
field theory the macroscopic inflaton field amplitude 0(i) corresponds to the 
expectation value of a (real) scalar Heisenberg field operator <£>(£, x): 

0(f) = ($(f, x)) = Tr {ft, $(f, x)} , (1.14) 

which involves the density operator qq as will be explained in more detail in 
section [2j After inflation ends at some time, which we may set to f = 0, the 
initial state for the subsequent evolution is approximately described by a vacuum- 
like, pure-state density operator £> 0 = Q(t = 0) with trace Tr £> 0 = Tr gl = 1. 
Typical physical normalizations of the inflaton model employ an initially large 
field amplitude 0(£ = 0) ~ 0.3 Mp with a very weak coupling of about A ~ 10~ 13 . 
We note that the fluctuating quantum field $(f, x) explicitly depends on the 
spatial variable x, while its expectation value 0(f) is homogeneous in our case. 

ft is a central topic of these lectures to show how the transition from a 
pure state described by a coherent macroscopic field amplitude proceeds towards 
radiation, thus connecting to the thermal history of the universe at later times. 
This can involve dramatic far-from-equilibrium phenomena at intermediate times. 
Indeed, the situation of the classically oscillating macroscopic field 0(f) turns out 
to be unstable in a quantum world, which will be shown in detail in section 0 
Here we will outline main aspects of the dynamics in order to introduce some 
important notions of nonequilibrium quantum field theory. Seeded by quantum 
fluctuations, a parametric resonance instability occurs. This instability is signaled 
by an exponential growth of correlation functions, which corresponds to very fast 
particle production called preheating. 

While the macroscopic field is described by a one-point function ($(t,x)), 
particle production can be extracted from two-point correlation functions. 
More precisely, one considers the symmetrized correlator ($(f, x) < f>(t', x') + 
$(f', x , ) < F(f, x)), which is the anti-commutator expectation value of two fields 
evaluated at equal times f = t!. The spatially homogeneous correlation function 
may be used to define a time-dependent occupation number distribution /(f, p) 
of particle modes with momentum p after Fourier transformation: 

/(f A? + , 1/2 + (27T)V(p)^(«) = \ fd 3 x <T"“ (t(i, x)4 (t, 0) + #(«, 0)4(1, x)> 

(1.15) 

for a given (in general time-dependent) dispersion relation ca(t, p). Because of 
spatial isotropy, the functions f(t, p) and cu(f,p) depend on the modulus of 
momentum, and we will frequently write /(t, |p|) etc. Here the term ~ 0 2 (t) 





coming together with the Dirac ^-f un ction denotes the coherent held part of 
the correlator at zero momentum. We note that in a finite volume V, we have 
(27 t) 3 <5(0) —y V and the scaling of the held term with volume just rehects the 
presence of a coherent zero mode spreading over the entire volume. For the 
following, all integrals will be considered to be suitably regularized, say, by some 
high-momentum cutoff if necessary. 

After inflation, all previous particle content is diluted away and the absence 
of particles corresponds to an initial state with f(t — 0 , p) = 0 for the subsequent 
evolution. Apart from the coherent held zero-mode, in this case only the 
“quantum-half” coming from the vacuum in (11.151) contributes. These quantum 
huctuations seed the parametric resonance instability, which leads to a decay of 
the macroscopic held 0 together with an exponentially growing mode occupancy 
after inflation: 

(1.16) 

Here, 7 q denotes a real and positive growth rate for some range of momenta 
around a fast growing characteristic mode Q. The far-from-equilibrium dynamics 
of an instability with exponentially growing modes may be seen as the opposite 
of what happens during the relaxation of some close-to-equilibrium modes, which 
decay exponentially. It is important to note that for transient phenomena the 
presence of instabilities is as common as relaxation is for the late-time approach 
to thermal equilibrium. 

The exponential growth (11.161) leads to an unusually large occupancy of modes 
at the characteristic momentum scale Q. Most importantly, this over-occupied 
system becomes strongly correlated despite the presence of a very weak coupling 
A C 1. Qualitatively, this may be understood from a “mean-held” approximation 
for quantum corrections. I 11 this approximation, the quartic interaction term in 
the potential (11.91) is replaced in the quantum theory by 


f(t,Q) ~ exp ( 7 Q t) 


-> j 6 (<* 2 (i,x))<* 2 (i.x). (1.17) 

Here, ($ 2 (f, x)) means the anti-commutator expectation value appearing also 
in (j 1.15 f) . now evaluated at equal space-time points. Applying proper 
renormalization conditions as in standard vacuum quantum held theory, we 
concentrate on the finite part of the term A(<F 2 (f, x)) given by 


A 


d 3 P f(t, p) 
(27r) 3 U)(t, p) 



(1.18) 


according to (11.151) . 

Strong correlations may be expected when interaction effects become 
comparable in size to the contributions from the typical kinetic energies of 
particles. For our case of particles with characteristic momentum Q, this happens 
at the time f* when their occupation number /(f*,Q) grows so high that their 
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contribution to the mean-field shift (II . 18jl is of the same order than their kinetic 
term ~ Q 2 . The latter is characteristic for the relativistic theory, which is second- 
order in space-time derivatives. Altogether we thus have 



d 3 P f(U, p) ^ 2 

(2vr) 3 o;(A,p) 


(1.19) 


Parametrically, we can estimate the integral for a relativistic dispersion u ~ |p 
as 



( 1 . 20 ) 


This is of order ~ Q 2 if the occupancy is as large as 


/(V Q) ~ t 


( 1 . 21 ) 


Taking into account (1 1.16 f) . such an extreme nonequilibrium condition is expected 
to occur at the time 



t, ~ — In - (1.22) 

IQ W 


after inflation. 

The maximally amplified mode turns out to be given by Q r-u V\<Kt = o) 


and approximately Q/^/q ~ 0(10). It is striking that the quantum correction 
in (11.1 7 jl cannot be neglected no matter how small the coupling A is. Though 
this is a robust feature of the preheating dynamics, we will see in section [5] that 
for a quantitative analysis one has to go beyond a mean-field analysis. In par¬ 
ticular, strongly nonlinear quantum corrections have to be taken into account 
even before t*, which are not included in the mean-field approach. Apart from 
the primary growth (II. 16jl starting early, nonlinear effects lead to an important 
secondary amplification period of fluctuations with enhanced growth rates for 
higher momentum modes. As a consequence, a wide range of growing modes lead 
to a prethermalization of the equation of state, given by an almost constant ratio 
of pressure over energy density, at the end of this early stage while the distribu¬ 
tion function itself is still far from equilibrium. 

Nonthermal fixed points and turbulence: The exponentially fast transfer 
of the zero-mode energy into fluctuations during preheating stops around £*, 
when the energy densities stored in the coherent field and in fluctuations become 
of the same order. Quantum-statistical corrections far beyond the mean-field 
approximation govern the subsequent dynamics of the strongly correlated system. 
If everything is strongly correlated, one might naively expect a fast evolution. 
However, the opposite turns out to be true and the evolution slows down 
considerably after t*. In contrast to the earlier instability regime, which exhibits 
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a characteristic scale for the growth of fluctuations, the subsequent evolution is 
governed by a nonthermal fixed point. This is an attractor solution with self¬ 
similar scaling behavior of correlation functions. 

A somewhat similar phenomenon is well-known close to thermal equilibrium, 
where strong correlations near continuous phase transitions lead to scaling 
behavior and a corresponding critical slowing down of the dynamics. However, 
the far-from-equilibrium situation after preheating is rather different. The 
typical momentum of a relativistic system in thermal equilibrium is given by the 
temperature T, and the Bose-Einstein distribution is of order one for momenta 
| p | ~ T. Moreover, thermal equilibrium respects detailed balance , where each 
process is equilibrated by its reverse process as will be discussed in section [3j In 
contrast, the strong overoccupation (11.211) of typical momenta ~ Q represents an 
extreme nonequilibrium distribution of modes, which leads to a net flux of energy 
and particles across momentum scales, thus violating detailed balance. 

To get a first understanding of this point, we may compare the characteristic 
energy density e /«dV(2ir) 3 |p|/(p) Q A /A of the far-from-equilibrium 

state to a Stefan-Boltzmann law ~ T 1 in thermal equilibrium. One observes 
that the nonequilibrium distribution exhibits a smaller characteristic momentum 
Q ~ A 1/,4 T than the corresponding thermal system for the same energy density. 
Therefore, we expect a transport of energy from lower to higher momentum scales 
on the way towards thermalization. Indeed, the underlying mechanism for the 
slow dynamics is related to the transport of conserved quantities, which leads to 
power-law behavior as will be discussed in detail in section [6j 

The situation is further complicated in the presence of additional conserved 
quantities, apart from energy conservation. Though total particle number is a 
priori not conserved for the relativistic system, an effectively conserved particle 
number emerges during the transient nonequilibrium evolution. To get a simple 
illustration of the consequences of both energy and number conservation, one 
may think of an idealized three-scales problem: We consider particles that are 
“injected” at the characteristic momentum scale Q as a consequence of the decay 
of the coherent held as described above. One then asks for the particles and the 
energy transported to a characteristic lower momentum scale K and to a higher 
scale A with K <C Q <C A. We let fiQ denote the injection rate for the number 
density at the scale Q , while hx and h\ characterize the corresponding “ejection” 
rates at K and A. Then number and energy conservation in this simple setup 
correspond to 

n Q = h K + n x , Quq = Kn K + A h A . (1.23) 

These can be solved for hx and h\ as 

A -Q Q-K Q 

nK = A _ K n Q- n Q , ri A = A _ K n Q ~-n Q) (1.24) 

where the latter approximate equalities exploit the separation of scales. 
According to the first equation the particles are transported to lower scales, while 
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Figure 2: Schematic behavior of the occupation number distribution near the 
nonthermal fixed point as a function of momentum |p| for two times t = t\ with 
f 2 > t\. The scaling exponents a and f3 characterize the rate and direction of 
the inverse particle cascade, while a' and /?' describe the direct energy cascade 
during the self-similar evolution. (A. Pineiro Orioli) 

the second says Aua ~ Quq such that the energy transport occurs towards higher 
scales. 

To illustrate this further, figure [2] shows the schematic behavior of the 
distribution function f(t , p) in the vicinity of the nonthermal fixed point. The 
distribution is given as a function of momentum for two subsequent times t — t\ 
and f 2 > ti. One observes that the particle transport towards low momenta is 
part of a dual cascade , in which energy is also transfered towards higher momenta. 
It turns out that there is no single scaling solution conserving both energy and 
particle number in this case. Instead, a dual cascade emerges such that in a given 
inertial range of momenta only one conservation law governs the scaling behavior. 

The inverse particle cascade occurs in a highly occupied infrared range of 
momenta with |p| <C Q. In this nonperturbative regime with /(f, p) > 1/A, 
the effectively conserved particle number receives a dominant contribution at the 
time-dependent scale K(t), where the integrand for particle number |p| 2 /(*,p) 
has its maximum. The cascade exhibits an approximate power-law behavior 
in an inertial range of momenta. This power-law is similar to the one known 
from stationary wave turbulence in the presence of external sources or sinks. 
However, we are dealing with an isolated system. Instead of simple power laws, 
the transport in isolated systems is described in terms of the more general notion 
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of a self-similar evolution , where the distribution function acquires a scaling form: 


(1.25) 

Here, all quantities are considered to be dimensionless by use of some suitable 
momentum scale. The values of the scaling exponents a and P, as well as 
the form of the time-independent fixed point distribution fsi^p) are universal. 
More precisely, all models in the same universality class can be related by a 
multiplicative rescaling of t and |p|. Quantities which are invariant under this 
rescaling are universal. 

The exponents for the infrared particle cascade are well described by 

a = fid , ^, (1.26) 

where d denotes the spatial dimension and here d = 3. The positive values for a 
and fi determine the rate and direction of the particle number transport towards 
low momenta, since a given characteristic momentum scale K{t\) = K\ evolves as 
K(t) = Ki(t/tf)-P with amplitude f(t, K(t )) ~ t a according to (11.251) . The fixed 
relation between a and P reflects the conservation of particle number density 
n = Jd d pf(t,p) ~ t a d/3 in this inertial range by using the self-similarity (ll.25j) . 

The emergence of an effectively conserved particle number for infrared modes 
can be explained by the dynamical generation of a mass gap for the (massless) 
relativistic theory, which arises from an intriguing interplay of condensation 
and medium effects. Indeed, the inverse particle cascade towards the infrared 
continuously (re-)populates the zero-mode, which leads to a power-law behavior 
also for the condensate held fi(t). In particular, if we would start the evolution 
from overoccupation (II.21)1 without an initial coherent held part (0(t = 0) = 0), 
the inverse cascade would lead to the formation of a Bose condensate in this 
far-from-equilibrium state, which is shown in section O 

The particle transport towards low momenta is accompanied by energy 
transport towards higher momenta. In figure [2] we indicate that in the direct 
energy cascade regime other scaling exponents a' and P' with a different scaling 
function f' s are found than for the inverse particle cascade. In the perturbative 
higher momentum range we have 1/A /(p) 1 and the energy is dominated 

by the scale A(f). For the hard modes with |p| ~ A the negative exponent 
P = —1/5 determines the evolution of characteristic momenta and a = —4/5 of 
their occupancy. The latter exponent determines the parametric time t ~ f*A~ 5 / 4 
at which the occupancies of the hard modes /(f, A(f)), which are order 1/A 
at t*, dropped to become comparable to the “quantum-half”. At this stage 
a subsequent approach to a thermal equilibrium distribution with the help of 
elastic and inelastic scattering processes sets in. 


f(t, P) = t a f s {t p p). 
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1.4 Ultracold quantum gases far from equilibrium 

Experimentally, the investigation of quantum systems in extreme conditions is 
boosted by the physics of ultracold quantum gases. Because these gases are very 
cold and dilute, they can be used to obtain clean realizations of microscopic 
Hamiltonians built from low-energy elementary scattering and interactions with 
applied fields. By using optical or atom chip traps, they provide a flexible testbed 
with tuneable interactions, symmetries and dimensionality, with connections to 
a wide variety of systems. Setups employing ultracold quantum gases can be 
largely isolated in contrast to many other condensed matter systems, where 
significant couplings to the environment are mostly inevitable. This offers the 
possibility to study fundamental aspects of far-from-equilibrium dynamics and 
the thermalization process in a highly isolated environment. 

We have seen in section 11.31 that the dynamics of the early universe after 
inflation can have different characteristic stages. After a very rapid evolution 
caused by nonequilibrium instabilities, the dynamics becomes self-similar. The 
time evolution in this self-similar regime is described in terms of universal scaling 
exponents and scaling functions. The universality observed opens the exciting 
possibility to find other systems in the same universality class, which may be 
easier accessible experimentally. Ultracold quantum gases at nanokelvins can 
differ in their characteristic energy scale by more than 38 orders of magnitude 
from inflationary physics. Nevertheless, they can show very similar dynamical 
behavior and universal properties can even be quantitatively the same. Therefore, 
one may learn from experiments with cold atoms aspects about the dynamics 
during the early stages of our universe. 

We will consider here interacting bosons with s-wave scattering length a, and 
leave the inclusion of fermions to later sections. For a gas of density n the average 
interatomic distance is n -1 /' 3 . Together with the scattering length a, this can be 
used to define a dimensionless “diluteness parameter” 

( = vW. (1.27) 

For a typical scattering length of, e.g., a ~ 5 nm and bulk density n ~ 10 14 cm -3 
the diluteness parameter ( ~ 1CU 3 is very small, and in the following we will 
always assume ( <C 1. The density and scattering length can also be used to define 
a characteristic “coherence length”, whose inverse is described by the momentum 
scale 

Q = y/l6nan. (1-28) 


To observe the dynamics near nonthermal fixed points for the interacting Bose 
gas, an unusually large occupancy of modes at the inverse coherence length scale 
Q has to be prepared. More precisely, for a weakly coupled gas of average density 
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n = f d 3 p/(27r) 3 / nr (p) this requires a characteristic mode occupancy as large as 



(1.29) 


This represents a far-from-equilibrium distribution of modes, and we note that 
the diluteness parameter ( plays the corresponding role of the coupling in (1 1.21 f) . 
Such an extreme condition may be obtained, for instance, from a quench or 
nonequilibrium instabilities similar to the previously discussed relativistic case. 

Most importantly, the system in this overoccupied regime is strongly 
correlated. These properties may be understood from a Gross-Pitaevskii equation 
for a nonrelativistic complex Bose held y: 

i<9 t y(t,x) = g\x(t,x.)\ 2 ^j y(t,x). (1.30) 

Here the coupling g is not dimensionless and determined from the mass m and 
scattering length as g = Ana/m. The total number of particles is given by 
f d 3 x\x{t, x)| 2 and is conserved. 

In the mean-held approximation the effect of the interaction term in the Gross- 
Pitaevskii equation is a constant energy shift for each particle, 

= 2g(|y| 2 ) = 2 gn = 2g J ^^/ n r(p) , (1-31) 

which can be absorbed in a redefinition of the chemical potential. However, we 
note that for the very high occupancy (11.2911 of the typical momentum Q the 
shift in energy is not small compared to the relevant kinetic energy Q 2 /2m, i.e. 
2 gn rs-/ Q 2 /2m. Parametrically, this can be directly verihed using HI,29jl : 

9 [ d 3 P /nr (p) ~ 9 Q 3 fnr(Q) ~ 9^T ~ 9 ~VT ~ ~ ■ (1-32) 

J Q mgQ m 

Here we have used that with a = mg/(An) equation (11,28ft implies Q = 2^/mgn 
and (II,27ft gives £ = mgQ/(16 tt 3 / 2 ). Most importantly, the energy shift 2 gn is 
of the order of the kinetic energy Q 2 /2m irrespective of the coupling strength g. 
This already hints at a strongly correlated system, where the dependence on the 
details of the underlying model parameters is lost. 

To study the nonequilibrium evolution of such a system, one may define a 
time-dependent occupation number distribution / nr (t, p) from the equal-time 
two-point correlation function: 

/nr(tp) + ^ + (2 tt) 3 5(p)|xo| 2 (*) = \ J d3 ^e _ipx (y(f,x)y*(f,0) + y(f, 0)y*(f, x)). 

(1.33) 
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Rescaled momentum: x = Pp 

Figure 3: Shown is the normalized fixed point distribution fs(t^ p) = t a f(t, p) in 
the infrared regime of the inverse particle cascade. The simulation results for the 
relativistic (circles) inflaton model with quartic interaction and the nonrelativistic 
(squares) Gross-Pitaevskii field theory agree to very good accuracy. 


The term | Xo | 2 (^) coming together with the Dirac (5-function determines the 
condensate fraction at zero momentum. These definitions are in complete 
analogy to equation (1 1.15 j) for the relativistic theory. The only major difference 
is the appearance of the dispersion in the definition for the relativistic case, 
which is a consequence of the second-order differential equation (j 1.13 jl while the 
nonrelativistic equation (11 .30j) is first-order in time. 

Since the far-from-equilibrium dynamics of the strongly correlated system is 
expected to become insensitive to the details of the initial conditions, we may 
choose the initial condensate fraction to be zero. Moreover, it turns out that 
in the overoccupied regime, where / nr (t, Q) is much larger than the “quantum- 
half”, the quantum-statistical dynamics is essentially classical-statistical such 
that the quantum and the corresponding classical field theories belong to the 
same universality class. The classical aspects of nonequilibrium quantum theories 
will be addressed in section [4j 

Figure [3] shows a striking example of universality far from equilibrium: 
Displayed is the normalized fixed point distribution /s,nr(^p) = £ _ “/ n r(^p) 
in the infrared, after evolving from overoccupied initial conditions for the 
nonrelativistic theory (circles) using the simulation techniques of section SI The 
corresponding fixed point distribution can be extracted for relativistic theories, 
such as discussed in section [I73l which we can generalize to include ^-component 
inflaton models as done in section [5j If the relativistic and the nonrelativistic 
theories belong to the same universality class, then the normalized fixed point 
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distribution has to agree. That this is indeed the case in the infrared scaling 
regime is exemplified in figure |3] for the example of the relativistic N — 2 
component field theory (squares). 

A crucial ingredient for the universality observed is the presence of strong 
correlations, which can occur even for weakly coupled systems if they are 
overoccupied. These extreme occupancies of modes can be found in a variety of 
isolated systems in extreme conditions and the identification of nonequilibrium 
universality classes represents a crucial step for their understanding. This is of 
particular relevance also for nonequilibrium gauge theories, such as employed 
in the description of the initial stages of ultrarelativistic heavy-ion collision 
experiments. Though we cannot cover this in depth in these lectures, we suggest 
some further reading at the end of each chapter. 
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2 Nonequilibrium quantum field theory 

2.1 How to describe nonequilibrium quantum fields 

There are very few ingredients for the description of nonequilibrium quantum 
fields. Nonequilibrium dynamics typically requires the specification of an initial 
state at some given time to- This may include a density operator qq = g(to) in a 
mixed (Tr of < 1) or pure state (Tr of = 1). Nonequilibrium means that Qo does 
not correspond to a thermal equilibrium density operator: go f with, for 
instance, g ^ ~ e~ l3H for the case of a canonical thermal ensemble with inverse 
temperature (3 for given Hamilton operator H. 

Completely equivalent to the specification of the initial density operator p 0 is 
the knowledge of all initial correlation functions: the initial one-point function 
Tr {^o$(to, x)}, two-point function Tr {p 0 $(to 7 x)<f>(to, y)}, three-point function 
etc. To be specific, here <f)(a:) may denote a scalar Heisenberg field operator 
depending on time and space, x = (x°,x). Typically, the “experimental setup” 
requires only knowledge about a few lowest correlation functions at the time to, 
whereas complicated higher correlation functions may build up at later times. 
From the behavior of the correlation functions for times x° > to one can extract 
the time evolution of all other quantities. 

Once the nonequilibrium initial state is specified, for closed systems the time- 
evolution is completely determined by the Hamiltonian. Equivalently, for a given 
classical action S the dynamics can be described in terms of a functional integral. 
From the latter one obtains the effective action T, which is the generating 
functional for all correlation functions of the quantum theory. There are no 
further ingredients involved concerning the dynamics than what is known from 
standard vacuum quantum field theory. It should be stressed that during the 
nonequilibrium time evolution there is no loss of information in any strict sense. 
The important process of thermalization is a nontrivial question in a calculation 
from first principles. Thermal equilibrium keeps no memory about the time 
history except for the values of a few conserved charges. Equilibrium is time- 
translation invariant and cannot be reached from a nonequilibrium evolution on 
a fundamental level. It is striking to observes that the evolution can go very 
closely towards thermal equilibrium without ever starting again to deviate from 
it for accessible times. The observed effective loss of details about the initial 
conditions can mimic very accurately the irreversible dynamics obtained from 
effective descriptions in their range of applicability. 

For out-of-equilibrium calculations there are additional complications which 
do not appear in vacuum or thermal equilibrium. A major new aspect for 
approximate descriptions concerns secularity: The standard perturbative time 
evolution suffers from the presence of spurious, so-called secular terms, which 
grow with time and invalidate the expansion even in the presence of a weak 
coupling. Here it is important to note that the very same problem can 
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appear as well for nonperturbative approximation schemes such as standard 
1/N expansions beyond leading order, where N denotes the number of held 
components. To obtain a uniform approximation in time, where an expansion 
parameter controlling the error at early times is valid also at late times, requires 
selective summation of an infinite series of contributions. 

If a nonequilibrium initial state is evolved from early to very late times, it is 
in general crucial that energy is conserved exactly for the given approximation. 
If thermal equilibrium is approached, then the late-time result is insensitive to 
the details of the initial conditions and becomes uniquely determined by the 
energy density and further conserved charges. This property can be conveniently 
implemented, if the dynamics is obtained from an effective action by a variational 
principle. The analogue in classical mechanics is well known: if the equations of 
motion can be derived from the principle of least action, then they will not admit 
any friction term without further approximations. 

Many of these requirements can indeed be achieved using efficient functional 
integral techniques based on so-called n-particle irreducible (nPI) effective 
actions. They often provide a practical means to describe far-from-equilibrium 
dynamics as well as thermalization from first principles. Closely related are 
approximations based on identities relating different n-point functions such 
as Dyson-Schwinger equations, or identities between derivatives of correlators 
described by the functional renormalization group. 

Local symmetries as encoimtered in gauge theories pose even stronger 
restrictions on possible approximation schemes. Here, the derivation of a sequence 
of effective theories, which disentangle dynamics on different time and distance 
scales can be very efficient. In particular, for weakly-coupled theories close to 
equilibrium at high temperature a wealth of information can be obtained from 
effective descriptions such as kinetic theory. 

Some aspects of nonequilibrium systems are often successfully described 
using classical-statistical held theory methods. Important examples involve 
the dynamics of nonequilibrium instabilities such as parametric resonance, or 
universal properties of dynamic critical phenomena near second-order phase 
transitions. Classical-statistical held theory can be implemented on a space- 
time lattice in Minkowski space and simulated on computers. Classical Rayleigh- 
Jeans divergences and the lack of genuine quantum effects limit their use. It is, 
therefore, important to understand for which time and distance scales classical 
simulations can adequately reproduce the behavior of the underlying quantum 
theory of interest. 

Of course, this list of methods is not complete and these notes are restricted 
to only some of the major held-theoretic techniques and applications. In the 
following we will introduce the functional integral formulation of nonequilibrium 
quantum held theory, which provides a very efficient starting point to derive the 
different approaches that will be discussed in later chapters. 
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2.2 Nonequilibrium generating functional 

We consider first a scalar quantum field theory with Heisenberg held operator 
$(x). Though most of the developments are completely equivalent for relativistic 
as well as nonrelativistic theories, in the following we will always assume a 
relativistic setup for definiteness and point out the relevant differences where 
necessary. 

Knowing all correlation functions, which are represented as the trace over 
time-dependent Heisenberg held operators for given initial density operator go, 
completely specihes the quantum system at any time. Therefore, all information 
about nonequilibrium quantum held theory is contained in the generating 
functional for correlation functions 

Z[J,R-g 0 ] = Tr (2.1) 

for given g 0 . It is the generalization of the partition function for nonequilibrium 
systems in the presence of sources J{x) and R(x,y). Since this is a generating 
functional for the trace over Heisenberg held operators, the time argument x° 
of the held operator $(a; 0 ,x) is evaluated along a closed real-time contour C 
appearing in the integrals, where we write f c = f c dx° f d d x in d spacial 
dimensions. This will be explained in more detail in section 12.31 and should be 
distinguished from computations restricted to scattering matrix elements between 
asymptotic incoming and outgoing states, which can be formulated without a 
closed time contour C. Employing the latter is more general, in the sense that 
vacuum quantum held theory as well as thermal and nonequilibrium systems 
can be described. The time-path C starts at time t 0 running forward along the 
real-time axis and then runs back to t 0 . Graphically, C is depicted as 



C+ 

* N 

to 

^ J 

c~ 00 


where we shifted parts of the curve slightly away from the real axis purely 
for visualization purposes. Contour time ordering along this real-time path is 
denoted by Tc . It corresponds to usual time ordering along the forward piece C + 
and reversed ordering on the backward piece C~ . In particular, any time on C~ 
is considered later than any time on C + . 

We have introduced the generating functional with two source terms on the 
contour, J{x) and R(x,y), which can be extended straightforwardly to take into 
account further source terms if necessary. Standard functional differentiation 
with respect to sources is extended to include time arguments along the closed 
time path. In particular, 

= -v) = ^( x ° - y°) <K x - y), ( 2 - 2 ) 
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where the Dirac Sc (x° — y°) is dehned on the closed time contour to be zero 
everywhere except at x° = y° if either x° and y° are both on C + or both on C~, 
where it is infinite with 


dx 0 5c(x°) = 1- 


(2.3) 


Setting all source terms to zero in (12.ip we obtain the normalized partition 
sum 

Z[J,R;Qo]\j , r=0 = Tr {po} = 1 . (2.4) 


Nonequilibrium correlation functions correspond to expectation values of 
products of Heisenberg held operators. These can be obtained by functional 
differentiation of (12. 1 p . For instance, the one-point function or macroscopic field 
(b reads 


5Z[J,R\ go] 
iSJ(x) 


j,r =o 


Tr{po<F(:c)} 


(<F(a:)) = (j>{x) . 


(2.5) 


Two-point, functions are obtained from second functional derivatives and n- 
point functions involve n derivatives. Together with the time-ordered two-point 
function 


5 2 Z[J,R-q 0 } 


iSJ(x)iSJ(y ) 


j,r =o 


= Tr{p 0 T c $(a;)$( 2 /)} = (T c $(a;)<%)> 


we can introduce the connected two-point function or propagator as 


( 2 . 6 ) 


G(x, y) = <T C *(i)*(tO) - <t>(x)m . (2.7) 


Here G(x,y) is dehned on the contour C and the time-ordering of Tc$(x)<F(|/) 
has to be evaluated according to the positions of the time arguments of the helds. 
In terms of the Heaviside step function 0(x° — y°) this reads explicitly 


T c $(x)$(y) 


{ <F(a;) < I)(|/) 9(x° — y°) + $(?/)$(a;) 9(y° — x°) for x°, y° on C + 
$(x)$(y) 9(y° — x°) + &{y)&{x) 9(x° — y°) for x°, y° on C~ 
&(y)$(x) for x° on C + , y° on C~ 

$(x)$(y) for x° on C~, y° on C + 

$(x)$(y) 9 c (x° - y°) + $(y)$(x) 9 c (y° - x°) , (2.8) 


where the last equation defines the contour step function 9 c (x° — y°). 

For an efficient notation, the held <E>(a;) may be written as < h ± (x°,x) where 
the ±-index denotes on which part of the contour the time argument x° is 
located. Then the contour integration for the linear source term in (I2.ip takes 
the form 


<&(x)J(x) 


' x,C 



dx° 


d d x ($ + (a:) J + (x) — $ (x)J (x)) 


(2.9) 
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where the minus sign comes from the reversed time integration along C~. 
Similarly, the bilinear source term in (12.ip can be decomposed into four terms 
introducing R ++ (x,y), R (x,y), R + ~(x,y ) and R~ + (x,y ) according to the 
different possibilities to locate x° and y° on C ± . For instance, the one-point 
function for vanishing sources is obtained as 


8Z[J, R ; p 0 ] 
i8J + (x) 


j,r=o 


4 > 0 ). 


( 2 . 10 ) 


We emphasize that the held expectation value is the same if obtained from a 
derivative with respect to either J + or J - in the absence of sources, since time¬ 
ordering plays no role for a one-point, function. Taking the second functional 
derivative of the generating functional (12.ip with respect to the source J + and 
setting J and R to zero afterwards, we obtain 


5 2 Z[J,R-,q 0 ] 


i8J + {x) i8J + (y) 


J,R=0 


(d >(x)<&(y) 9(x° — y°) + &(y)$(x) 9(y° — x °)) 
G ++ (x,y) + <f>(x)</>(y). (2.11) 


Here we used that contour ordering corresponds to standard time ordering if all 
time arguments are on C + . We also introduced the notation G ++ (x, y) in order to 
distinguish this correlator from the other possible second functional derivatives 
with respect to the sources J + and J~ setting J, R = 0 afterwards. These can be 
written as: 


5 2 Z[J,R-q 0 } 


i8J ( x)i8J ( y ) 


j,r =o 


8 2 Z[J, R] p 0 ] 


i8J + (x ) i8J~(y ) 
8 2 Z[J, R] p 0 ] 


i8J ( x)i8J + (y ) 


j,r =o 


j,r =o 


(<F(a:)$(p) 9(y° — x°) + &(y)$(x) 9(x° — y °)) 
G~~(x,y) + ct)(x)(t)(y ), 

(<F(p)d>(x)) = G + ~(x,y) + (j>{x)(j>{y) , 

{<&(x)${y)) = G~ + (x,y) + 0(x)0(y) . (2.12) 


While any nonequilibrium two-point correlation function can be written in terms 
of the components (12.lip - (I2.12p . not all of them are independent. In particular, 
using the property 9(x° — y°) + 9(y° — x°) = 1 of the Heaviside step function one 
obtains the algebraic identity for J, R = 0: 


G ++ (x,y) + G (x,y) = G + ( x,y) + G + (x,y) . (2.13) 


Later we will observe that a further, so-called fluctuation-dissipation relation 
exists for the special case of vacuum or thermal equilibrium density matrices. 
However, for general out-of-equilibrium situations this is not the case. 
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2.3 Functional integral representation 

Above the time-ordering along the closed contour C appears as a bookkeeping 
device that allows one to conveniently describe different components of 


nonequilibrium correlation functions. To simplify the evaluation of correlation 
functions, in the following we write the generating functional (12.lj) in terms of 
a functional integral representation. In this construction the closed time-path 
appears because we want to compute correlation functions which are given as the 
trace over the density operator with time-ordered products of Heisenberg held 
operators. Representing the trace as a path integral will require a time path 
where the initial and final times are identified. 

We evaluate the trace using eigenstates of the held operator <f )± at time t 0 , 

**(*«>, x) 1^) = ipf (x) 1^) , 


(2.14) 

such that (j2.1|) may be written as 



Z[J,R;£o] = l[dtp£](ip + \ g 0 T c expi 

&(x)J(x) 

x,C 


+ \ [ $(x)R(x,y)$(y) 

^ Jx,y,C 

1 1<^ + } ■ 

(2.15) 

The integration measure is 



j [d^o i = f • 


(2.16) 

With the insertion 

J I d Fo] \V~)(V~\ = 1 , 
we mav briner (12.151) to a form 


(2.17) 

Z{J,R]Qo] = / [dyyj-][dv? 0 ] (<P + \ Qo W ) (<p , t 0 \ <f + , t 0 ) JR 

. (2.18) 


Here the transition amplitude in the presence of the sources is given by 


(<P ,to\v + ,to) JR = {v | Te exp i 


&(x)J(x) 


' x.C 


1 

2 


'x,y,C 


$(x)R(x,y)$(y) > \(f + ) . (2.19) 


This matrix element can be written as a functional integral over “classical” fields 
ip(x) on the closed time path C, whose derivation is given below. For a quantum 
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theory with classical action 5[</?] it reads 


(<d ,t 0 \<p + ,t 0 )j iR 


Here the classical action includes the time contour integral over the respective 
Lagrangian density £(x), i.e. S = f c dx° f d d xC(x). The functional integration 
goes over the field configurations <p(x) depending on space and times x° > t 0 
that satisfy the boundary conditions <p±{x Q = x ) = V^o ( x ) f° r x ° 011 The 
prime on the functional measure is indicating that the integration over the fields 
at x° = t 0 is excluded. 

Putting everything together, we may write the generating functional (12. 1 p as 


fo 


Vo 


?V ex P i ^ ‘S'M + / <p(x)J(x ) 
Jx,C 


1 

2 


(p(x)R(x,y)(p(y)\ . 


'x,y,C 


( 2 . 20 ) 


Vo 

Z[J,R,qo]= I [d9?+][dv?o](^ + |Po|^“) y'^e < { SM+ ^.c J(x)v(x)+ 3^,cV ( * )fl( *-J' ) ^ ( w ) } 

Vo 

" -V- - ---—--V-- 

( 2 . 21 ) 

initial conditions quantum dynamics 

The above expression displays two important ingredients entering nonequilibrium 
quantum field theory: the quantum fluctuations described by the functional 
integral with action S, and the statistical fluctuations encoded in the averaging 
procedure with the matrix elements of the initial density operator £o- 

Formulations of closed time path generating functionals typically employ a 
time interval starting at some to and extending to the far future as mentioned 
above. Causality implies that for any n-point function with finite time arguments 
the contributions of an infinite time path cancel for times exceeding the largest 
time argument of the n-point function. We will see this explicitly when we derive 
time evolution equations for correlation functions in section 12.61 To avoid un¬ 
necessary cancellations of infinite time path contributions we will often consider 
finite time paths. The largest time of the path can be kept as a parameter and 
is evolved in the time evolution equations. 

Details: Construction of the path integral 

To simplify the evaluation of correlation functions we wrote the above generating 
functional in terms of a functional integral representation using (12.201) . Here 
we describe the construction of the corresponding transition amplitude following 
standard presentations and refer to the literature for more information. 
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We start from a Schrodinger picture where the held operators are time- 
independent. For the scalar held theory the Hamiltonian H[U, $] is expressed in 
terms of the held <b(x) and the conjugate momentum held operator n(x). We 
work with basis sets of eigenstates for which $(x) or n(x) are multiplicative 
operators, respectively, 

$( X M = ^(x)l^), n(x)|7r) = 7r(x)|7r). (2.22) 

The completeness and orthogonality conditions read 

MVi) = 5[<Pj ~<Pi\ ( 2 - 23 ) 


(njlTTi) = S[7rj - TTi] . (2.24) 

To become familiar with the notation we recall that the <5-functional may be 
represented as a product of individual 5-functions, one for each point x in space, 

- (pi] = JJ 5 (<pj(x) - <pi(x)) . (2.25) 

X 

If we write each 5-function in the momentum representation 

5 (ifj(x) - <pi(x)) = j d7 ^ exp {iTTi(x) (^j(x) - <^(x))j (2.26) 

we obtain 

-ipi\ = J] f d?r 2 ^ exp |i7Ti(x) (^(x) - <^(x))j 

exp ji J d d x7Ti(x) (<Pj(x) - v?i(x))| . (2.27) 

Writing the same 5-functional using orthogonality and completeness conditions 
as 5[ipj — <Pi\ = j [d7Tj/(27r)] {(pj\TTi)(ni\(pi), we recover by comparison with (j2.27[) 
the familiar overlap 

{<Pj\ni) = exp ji J d d a;7r i (x)v?j(x)| . (2.28) 

If the dynamics is described by a Hamiltonian H then the state |(^ 0 ) at time 
to evolves into e~ lH( ' tf ~ t °' > \ipo) at time tf. The transition amplitude for the held 
in configuration <fio( x ) at time t 0 to evolve to at time tf is thus 

(¥>/, tf | 75), to) = {<p f \ \<p 0 ) . (2.29) 



[d(p\\ip)(<p\ = l , 


and 


d7T 

27T 


171") (7T | = 1 , 
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To conveniently evaluate this amplitude we write it as a path integral. For this 
we divide the interval tf — t 0 into N +1 equal periods. Next, we insert a complete 
set of states at each division and consider the limit N —> oo: 


N 


{<Pf,tf\<Po,t 0 ) - 


••• Me 


—iH(tN—tN- 1) 


Ivdv-i) 


2—1 
—iH(ti—to) I 


Vo)- (2.30) 

As N gets large, At = (U+i — U) 0 and we may expand the exponentials, 
(<f i+ 1 \e~ lHAt \ (pi) ~ (ipi + i\l-iHAt\(pi) 


= 6 [<p i+1 - w] - iAt (cp i+ 1| H\<pi). 


(2.31) 


The first term in the above sum is a <5-functional which can be written in terms of 
its momentum representation (12.271) . It remains to evaluate the operators in the 
matrix element (<fi+ 1 | AT[II, $] | tpi). We consider first the case of a matrix element 
for a contribution h n [ II] that only depends on the conjugate field momentum. We 
write, inserting a complete set of states, 

(<Pi+l\ Mn] |7Ti) (7T i\<Pi) 

MV exp ji J d d x7Ti(x) V+i(x) — M x ))j , (2.32) 



M+1| K[U] fa) = 


d7 n 


where (12.281) is employed for the second equality. Typical Hamiltonians can be 
written as a sum of a field and a conjugate momentum term, i.e. i/[n,$] = 
h v [T] + h n [n] . For the held contribution h v [T], we employ a symmetric operator 
ordering. In general, matrix elements including also products of $ and n 
operators can be evaluated by a symmetric operator ordering. For instance, 
the symmetrization of a term like $n —* {$, n}/2 = (<I>n + n$)/2 leads to the 
matrix element 


. . $n + n$ 

Wi+ll -o- 



Wi) = 


Hi exp 


fi +1 + <Pi 
2 





XI Tj(x) (^i + l(x) 



(2.33) 


where 


In general, we can write 




fi +1 + 74 

2 


(2.34) 


(ip i+1 \H[Il,$} | <pi) = 


d7 Tj 
2n 


H[Hi, ipi] exp < i / d d x7Ti(x) M +i ( x ) - 74 ( x )) 


(2.35) 
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Therefore, (12.311) may be represented as 


<¥>t+i|e lHAt \ipi) = 


d7 Ti 
2 vr 


(l-iAt H[Ki, ipi} + 0{At 2 )) 


x exp li J d d XTT j(x) (<p i+ i(x) - <#(x)) 

"d7Tj 


2tt 


exp < i / d d x7Ti(x) ( 9 ? i+ i(x) - tpi (x)) -iAtH^p^ \ . (2.36) 


Putting everything together gives 

. N N 

{iff, t f | (f 0 , t 0 ) = jin^ / JJ [d^] JJ 


i=l j =0 


d7Tj 

2tt 


“id 


x exp ^ ?At 2^ 1 / d< W x ) 
5=0 

V>(tf,x)=Vf(x) 


<Pj+ l(x) - ^-(x) 
At 


HfrjiVj] 

dp{x 


ip(to,x)=ipo(x) 


exp ^ d x ° (^J ^ dx7r i x ) ^^,J — H[n, </?]^ | (2.37) 


Above we have defined </ 2 at + i(x) = </?/(x) and identify (ifi(x),ti) = </?(£*, x) 
employed for the continuum notation. We emphasize that all references to 
operators are gone in (12.371) . 

If the Hamiltonian is quadratic in n{x), as for a real scalar field theory with 
H[tt, <f] = I d d x Q vr 2 + -^ (W) 2 + <f 2 + V {(p)^ (2.38) 


with mass parameter m and interaction part V{p), then one can carry out the 
functional integration over n in (12.371) . Completing the squares one finds from 
the Gaussian integral 

J ex P { i7r i(7 5 i+1 - <Pj) = (2 TriAt)-? exp j ^ + 2 A /^ | 

(2.39) 

that the net effect is to replace n by the time derivative of p in the exponential 
of (12.371) . Applying the continuum notation the matrix element then becomes 


{iff, tf | p 0 , t 0 ) 


<p(t f ,x)=<p f (x) 



<p(to,x.)=<po(x) 
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where the classical action reads 


= i' dx ° I A “ x I \ (!?) “ \ (Vv? - t ^ 2 ■ VM 


(2.41) 


Here the measure means 


r r N 

/ Sl'tp = lim / dy9(tj, x) (27riAt) _ 2 ; NAt = tf — t 0 . 

J J i=1 x 


(2.42) 


The final time tf in the above matrix elements is, of course, arbitrary. 
Applying the same construction to the closed time path of section [2T2l first to the 
forward piece C + and subsequently to the backward piece C~, is straightforward. 
Adding to the Hamiltonian appropriate source terms as employed in section 12.21 
leads to the generating functional for correlation functions. 

2.4 Initial conditions 

To understand in more detail how the initial density matrix enters calculations, 
we consider first the example of a Gaussian density matrix for a real scalar field 
theory. For simplicity, we neglect for a moment the spatial dependencies setting 
d — 0, i.e. we consider quantum mechanics. The generalization to higher d will 
typically be straightforwardly In this case the most general form of a Gaussian 
density matrix can be parametrized in terms of five real parameters 0 O > 0o> r l 
and cr: 


{<P + \QoW ) = 



+ i ^ [(^0 - ^ o ) 2 - 


-.o)-^ 1 [(y?o -0o) 2 + Oo -0o) 2 ] 
Oo - 0o) 2 ] + a 4 ^ 2 1 (^o - 0o)(^o - 0 o)}, 


(2.43) 


whose values will be further discussed below. Gaussianity refers here to the 
fact that the highest power of <p q appearing in the exponential of (12.431) is two. 
Therefore, density matrix averages of field operators at initial time only involve 
Gaussian integrals. 

In order to interpret the above parameters and to see that we indeed consider 
the most general Gaussian density matrix, we first note that (12.431) is equivalent 
to the following set of initial conditions for one- and two-point functions: 


</>o = Tr {Q 0 ${t)} lt=to , 0o = Tr {p 0 cGI>(f)} Mo , (2.44) 

1 We note that for homogeneous field expectation values taking into account spatial 

dependencies essentially amounts to adding a momentum label in Fourier space (see also below). 
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£ 2 = Tr {go^(tymt r )} lt=t , =to -0000, (2.45) 

tv = ^ Tr {^o (dt$(t)$(lf) + ^(t)d t ^(t'))} lt=t , =to - 0 O 0 O , (2.46) 

2 

V 2 + ^2 = Tr {«oc»f4l>(0a t /4>(0)}| i=t , =io - 0o0o • (2.47) 

In contrast to the symmetrized combination (12.461) . we note that the anti¬ 

symmetrized initial correlator involving the commutator of <3> and d t & at to is 
not independent because of the operator commutation relation 

[<h(f),«9ph(t)] = i. (2.48) 

The equivalence between the initial density matrix and the initial conditions for 
the correlators can be verified explicitly. For instance 

/ OO 

-oo 

1 r°° ci 

= y!i?.L drfexp { _ ^ w “ 4) } = l! (2 ' 49) 

where one notes that the RHS of (I2.43jl does not depend on 0 O , a and rj for 
l Pq ={ Pq- Similarly, 

1 r°° 

Tr {p 0 = -^== j <p+ exp | - ^(<Po “ 0o) 2 } = 0o , (2.50) 

which is simply obtained by the shift <Pq —$■ <Po + 0o, etc. Similarly, since only 
Gaussian integrations appear one finds that all initial n-point functions with 
n > 2 can be expressed in terms of products of the one- and two-point functions. 

Of course, higher initial time derivatives are not independent as can be 
observed from the field equation of motion. For instance, for the above scalar 
theory Hamiltonian (12.38j) with mass m and the interaction part 14(4*) = A<F 4 /4! 
the corresponding field equation reads: 

(a 2 $) = -m 2 (<I>) _ — ($ 3 ). (2.51) 

6 

Since ($ 3 ) is given at initial time in terms of one- and two-point functions for 
Gaussian p 0 ; also second and higher time derivatives are not independent. We 
conclude that for our case the most general Gaussian density matrix is indeed 
described by the five parameters appearing in (12.431) . In particular, all observable 
information contained in the density matrix can be conveniently expressed in 
terms of the correlation functions (12.44 j) - (j2.47p . 
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For further interpretation of the initial conditions we note that 


Tr £>o — / dy2o"d</? 0 (+ + M + ){<P M+ + ) 

d<^o dy?Q exp j - [((^+) 2 + (<^-) 2 ] + 

(2.52) 

where again we shifted the integration variables by 0 O to arrive at the second 
equality. The latter shows that for cr > 1 the density matrix describes a mixed 
state, which may be stated in terms of a non-zero occupation number / with 
a = 1 + 2 /0 For a = 1 the “mixing term” in (I2.43P is absent and one obtains a 
pure-state density matrix of the product form, 

0 o = |tf><*|, (2.53) 



with Schrodinger wave function 

^ + l*> = (2 ; 2) i/4 ex P {^o+o + - (4^2 + * |r) ^0 - 0o) 2 } • ( 2 . 54 ) 


In order to go beyond Gaussian initial density matrices and to field theory in 
d spatial dimensions, one may generalize the above example and parametrize the 
most general density matrix as 

(+ + M+-> = Afe ihc ^, (2.55) 


with normalization Af and he [+] expanded in powers of the fields: 


hc[<p\ = «o+ / «i (x)<p(x) + - / a 2 (x,y)<p(x)<p(y) 

Jx,C " Jxy,C 


1 

3! 

1 

41 


«3 (x,y,z)<p(x)(p(y)(p(z) 


' xyz,C 


' xyzw,C 


a 4 (x, y, z, w)ip(x)ip(y)ip(z)(p(w) + ... (2.56) 


Here the integrals have to be evaluated along the forward piece, C + , and the 
backward piece, C~, of the closed time path with <p + (t 0 ,x) = <Pq(x) and 
(p~{t 0 ,x) = +0 (x) as described in sections 12.21 and 12.31 Since the density matrix 
£0 is specified at time to only, all time integrals in (12.56(1 can contribute only 
at the endpoints of the closed time contour. As a consequence, the coefficients 

2 For the special case of thermal equilibrium the occupation number is given by the Bose- 
Einstein distribution. 
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ai(x), a 2 (x, V ), otz(x, y, z), ... vanish identically for times different than f 0 . More 
precisely, one has 


ati(x)<p(x) 


' x.C 



«2 (x,y)<p(x)<p(y) 


J d d x «(x) ytf(x) + a x (x) <p 0 ( x )} , 

J d d xd d |/|a^ + (x,y)^+(x)99+(y) 

+ a+“(x,y) <p£(x)po(y) 

+ « 2 " + (x,y) (fio ( x )^o(y) 

+ a^'( x ,y)^o ( x )v?o(y)} (2.57) 


and so on. We note that «o appearing in (12.56ft is an irrelevant constant that 
can be taken into account by a rescaling of the overall normalization Af in (I2.55[) . 
For a physical density matrix the other coefficients are, of course, not arbitrary. 
Hermiticity of the density matrix, g 0 = Qq, implies iaf = ( ia 7)*, iaf + = (ia^ _ )*, 
ia\~ = (iaf + )*, etc. If the initial state is invariant under symmetries, there are 
further constraints. For example, for an initial state which is invariant under 
$ —> — <E>, all a n (xi,, x n ) with odd n vanish. If the initial state is homogeneous 
in space, the a n (x \,..., x n ) are invariant under space translations x, —y x, : + a for 
% — 1,..., n and arbitrary spatial vector a. In this case, they can be conveniently 
described in spatial momentum space, 


a. 


±- 


(X!,...,x n ) = 


d d pi d d p r 

(2vr) d 


x a^'"(pi,.. 


(2vr) d 

,Pn) • 


exp 


f ( 27F ) d 5 (Pl + • • • + Pn) 
3 =1 J 

(2.58) 


Using the parametrization (12.551) and (12.56ft for the most general initial density 
matrix, one observes that the generating functional (12.21ft introduced above may 
be written as 


Z[J, R ; go] 



R(x,y)ip(x)ip(y)+^i f xyz C a 3 {x,y,z)cp(x)v(y)<f(z)+...) 


(2.59) 

Here we have neglected an irrelevant normalization constant and rescaled the 
sources in (I2.2ip as J(x) + ot\{x) —> J(x) and R(x,y ) + a 2 (x, y) —* R(x,y). The 
sources J and R may, therefore, be conveniently used to absorb the lower linear 
and quadratic contributions coming from the density matrix specifying the initial 
state. This exploits the fact that the initial density matrix is completely encoded 
in terms of initial-time sources for the functional integral. 

The generating functional (12.591) can be used to describe situations involving 
initial density matrices for arbitrarily complex physical situations. However, 
often the initial conditions of an experiment may be described by only a few 
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lowest n-point functions. For many practical purposes the initial density matrix 
is well described by a Gaussian one. For instance, the initial conditions for the 
reheating dynamics in the early universe at the end of inflation are described by 
a Gaussian density matrix to rather good accuracy. Also many effective theories 
for the description of nonequilibrium dynamics, such as Boltzmann equations, 
exploit an assumption about Gaussianity of initial conditions. 

From (I2.59P one observes that for Gaussian initial density matrices, for which 
a;* = 0 for i > 3, one has 


Z 


J, R ; 


(Gauss) 


-> Z[J,R], 


(2.60) 


As a consequence, in this case the nonequilibrium generating functional 
corresponds to the generating functional with linear and bilinear source terms 
introduced in (12.11) for a closed time path. 

In field theory non-Gaussian initial density matrices pose no problems in 
principle but require taking into account additional initial-time sources. An 
alternative to initial-time sources is to represent a non-Gaussian initial density 
matrix with the help of an additional imaginary piece of the time contour 
preceding the closed time path C. We refer to the literature at the end of this 
chapter for further details concerning this alternative representation. 


2.5 Effective actions 

The functional Z[J,R-,g 0 \ given by (12.591) is the nonequilibrium quantum field 
theoretical generalization of the thermodynamic partition function in the presence 
of source terms. In thermodynamics, Legendre transforms of the logarithm of 
the partition function lead to equivalent descriptions of the physics. Similarly, 
in nonequilibrium quantum held theory Legendre transforms with respect to the 
various source terms lead to different representations of a free energy functional T. 
The Legendre transform with respect to the linear source term, ~ J, will lead to 
the so-called one-particle irreducible (1PI) effective action r[</>] parametrized by 
the one-point, function <ft. An additional Legendre transform with respect to the 
bilinear source term, ~ i?, gives the two-particle irreducible (2PI) effective action 
T [(f), G], which is parametrized in addition by the two-point function or propagator 
G. In the presence of higher source terms this procedure can, in principle, be 
continued to n-particle irreducible (nPI) effective actions up to arbitrarily high n. 

In the absence of approximations, all these generating functionals give 
equivalent descriptions of the physics and any choice is only a matter of 
convenience or efficiency. For nonequilibrium initial-value problems it is often 
particularly helpful to employ effective actions, since they are described in 
terms of correlation functions whose initial values are typically much better 
accessible than those of source terms. Furthermore, the ability to find suitable 
approximation schemes in practice can strongly depend on the appropriate choice 
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of the functional representation. In a following chapter it will be demonstrated 
how 2PI effective actions can be applied to describe important processes such 
as thermalization in quantum field theory. Later chapters include also further 
applications, such as nonequilibrium instabilities or an efficient derivation of 
kinetic theory in its range of validity. 

In this section we discuss the construction of 1PI and 2PI effective 
actions. The latter provides a powerful starting point for approximations in 
nonequilibrium quantum field theory with Gaussian initial density matrix as 
specified in section 12.41 It should be emphasized again that the use of a 
2PI effective action for Gaussian initial density matrices represents a priori no 
approximation for the dynamics — in an interacting quantum held theory higher 
irreducible correlations build up in general corresponding to non-Gaussian density 
matrices for times t > t 0 . It only restricts the possible setup described by the 
initial conditions for correlation functions. 

To be able to explain a variety of approximation schemes in a following 
chapter, we consider here the example a quantum held theory for a real, N- 
component scalar held cp a (a — 1,..., N) with 0(7V)-symmetric classical action 


*%] = 


' x.C 


Hi 

~d> X (p a (x)d fl ip a (x) ~ —(p a (x)(p a (x) 


X 


4!iV 


(Va(x)Va(x)) 2 > • (2.61) 


Here summation over repeated indices is implied as well as the integration over 
the closed time path C as introduced in section [2721 Starting from Z[J,R ], given 
by (I2.59P with cq = 0 for i > 3, we write 

Z[J,R ] = exp (iW[J, i?]) 


= expi< + / J a (x)ip a (x) + - Rab(x,y)(p a (x)ip b (y) > (2.62) 

J l JX,C 2 J X yfi J 

with real J a (x), and R ab (x,y) = R ba (y,x). 

The macroscopic field <f a and the connected two-point function G ab in the 
presence of the source terms ~ J a (x) and ~ R a b{x, y ) can be obtained by variation 
of W[J, /?,], 


-iSlnZ[J,R] 


5W[J,R] _ _ 

5J a {x) SJ a (x) 

5W[J,R] = 1 

5R a b(x,y ) 2 


Z[J,R} i6J a (x) 

M x )My) + Gab(x,y ) 


SZ[J ' R] = <Pa(x) , (2.63) 


(2.64) 


These definitions are in accordance with (12. 5 p and (12.71) for vanishing sources, 
i.e. J,R = 0. However, we keep J and R nonvanishing in (I2.64p and to ease 
the notation the dependence of <f = 4>(J,R) and G = G(J,R) on the sources is 
suppressed. 
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Before constructing the 2PI effective action, we consider first the 1PI effective 
action. It is constructed via a Legendre transform with respect to the source term 
linear in the field, 


= w[,f R] - 


' x,C 


5W[J,R] 

5J a {x) 


J a (x) = W[J,R}- / J a (x)<l> a {x). (2.65) 


' x,C 


To complete this construction one has to note that the relation between <f> and J 
is independent, i.e. inverting f = 8W[J, R]/5J yields J = J R (4>) on the RHS of 
(I2.65p . Taking the functional derivative with respect to the held one finds 


Sfa(x) 


'vfi 


5W[J,R] 5J b {y) 
5J b (y) S(f>a(x ) 





-Ja{x). ( 2 . 66 ) 


and the source derivative yields 

(5T R [0] _ 5W[J,R] | f 5W[J,R] 8J c (z ) 

5R ab (x,y) 5R ab (x,y ) J zfi 5J c (z ) 8R ab (x,y) 
5W[J,R] 

5R ab (x,y ) ' 



M z ) 


SJ c {z) 
SR ab (x,y ) 


(2.67) 


The second equalities both in (12.66(1 and (12.67(1 arise from cancellations using 

(HMD. 

The standard 1PI effective action, T[</>], is obtained from (12.65(1 for R — 0, i.e. 
T[0] = p r<=0 [((>]. For non-zero R the functional r fl [0] may be viewed as the 1PI 
effective action for a theory governed by a modified classical action, [</?], with 


S R [(f] = S[<p] + ]- f Rab(x,y)(p a (x)<p b (y). (2.68) 

Z J xy,C 

As a consequence, it is straightforward to recover for r R [</>] “textbook” relations 
for the 1PI effective action taking into account R. For instance, evaluating 
T B [(j)} in a saddle-point approximation (“one-loop”) around the macroscopic held 
configuration 0 a (x) one hnds up to irrelevant constants 


r R(lloop) [0] = S R [(f>] + ^ Tr c ln[Go x (0) - iR ] , (2.69) 

where the trace involves the integration over space-time coordinates on the closed 
time path C as well as summation over held indices. We derive (12.69(1 below, 
which is indeed the standard one-loop result for the 1PI effective action with 
S[(j)\ —> 5 R [0] and Gf 1 (</>) —> Gf 1 ^) — iR. Here the classical inverse propagator 
= 5 2 S[cf)]/8(j) a (x)5(j) b (y) for the action (12.61(1 reads 

iG oib( x >y'’<!>) = ** (o x + m 2 +-^ (f c (x)(j) c (x)j 5 ab 8c(x^y) 

~3^a( x )M x )^d x ~ y) ■ (2.70) 
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It is instructive to compare this to what is obtained from the approximation 
(I2.69P for the propagator by taking the functional derivative with respect to 
r. sr R ( n °°rt[(/)]/5R = (cjxj) + [Gq 1 (0) — iR]~ 1 )/2. Since this has to agree to 
(4>4> + G)/2 using (12.671) and the definition (12.651) . we obtain G~ l = Gq 1 (0) — iR 
in this approximation. 

We now perform a further Legendre transform of T R [</>] with respect to the 
source R in order to arrive at the 2 PI effective action , T[</>, G], which is a 
functional of the field (p a {x ) and the propagator G ab (x,y ): 


r[0,G] = r 


R r 


8T 


r r 


' xy,C 8Rg b { x i 'll) 

5W[J,R} 


R a b{x,y ) 


/ x = o [M x )Mv) + G ah (x,: 

dRab{x,y) 2 


- 2 


'xy,C 


[(fa(x)My) + G ab (x, y)} Rab(x , y) . (2.71) 


Here we have used first (I2.67P and then the definition (12.641) . Of course, the two 
subsequent Legendre transforms, which have been used to arrive at (I2.7ip . agree 
with a simultaneous Legendre transform of W[J,R] with respect to both source 
terms: 


T[<j>,G] = W[J,R} - 


' x.C 


6W[J,R] 

5J a (x) 


Ja{x) 


6W[J,R] 


Rab(x,y) 


= W[J, R] 


f a (x)J a (x) 


’ x,C 


1 

2 


'xy,C 


Ixy,C ^Rabix^y ) 

[4> a (x)(j) b (y) + G ab (x, y)\ R ab (x, y) .(2.72) 


Taking the field and propagator dependence of the sources into account, one 
obtains from (12.721) after cancellations the stationarity conditions 


5(f)a{x) 
5G ab (x,y ) 


' Ja (x~) 


Rab(x,y)</> b (y ), 


'y,c 


- -R ab (x,y ), 


(2.73) 

(2.74) 


which are the quantum equations of motion for <f> and G. 

To get familiar with the 2PI effective action, we note that the following 
approximate expression 


T (lloop) [</>,G] = ^[0] + l T r c lnG- 1 + ^Tr c {(G o - 1 (0)-G- 1 ) G} 

= T^ 11001 ^] - ]- [ [<f> a (x)(l> b (y) + G ab (x,y)}R ab (x,y) (2.75) 

" Jxy,C 
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indeed corresponds to the Legendre transform of the one-loop result (12.69ft for 
the 1PI effective action as written in the second line of (12.75j) . To see this, 
we first observe that the stationarity equation (12.741) yields in this case 


Ap(lloop) \(k O] j j 1 

JgJ^) ' = -2°-^ + = - 2 R^,V) (2.76) 

or G” 1 = Gq x (0) — iR. Using this to replace (Gq 1 (0) — G _1 ) in the first line of 
(I2.75|) yields with the definitions (I2.68|) and (I2.71|) the result. 

To go beyond this approximation it is convenient to write the exact T[0, G] 
as the one-loop type expression in the first line of (12.75[) and a “ rest ”, T 2 [</>, G]: 


T[0, G] — S[<t>] + -Tr c hr G 1 + — Tr^ G 0 1 (0) G + r 2 [0, G] + const 


(2.77) 


Here we have written Tr C jG^ 1 G as an irrelevant constant which can be adjusted 
for normalization. To get an understanding of T 2 [0, G] we vary this expression 
with respect to G, which yields 


G a<b(. x >y) = G oi(u2/;0) -iRab(x, y) - £ afe (x, ?/;</>, G) , (2.78) 

where we have defined: 


j: ab (x,y](j),G) = 2i 


^r 2 [0,G] 

5G ab (x,y ) 


(2.79) 


This gives the proper self-energy, T, ab (x,y, f,G), in terms of a functional 
derivative of T 2 [4>, G] with respect to G ab (x,y). 

Inverting (12.781) . we may express for further interpretation the full propagator 
G as an infinite series using a compact matrix notation: 

G = (Go 1 — iR)~ l + (Go 1 — iR)~ l E (Gq 1 — iR)~ l 

+ (Go 1 - iR r 1 S (Go 1 - iR) _1 S (Go 1 - iR)- 1 + ... (2.80) 

A term is called one-particle irreducible (1PI) if it remains connected when an 
arbitrary internal propagator (Gq 1 — iR)- 1 is removed. For instance, removing in 
the expression (G^ 1 — iR)- 1 £ (Gq 1 — iR)~ l £ (Gq 1 — iR)- 1 the inner propagator, 
we end up with two disconnected pieces (G^ 1 — iR)- 1 £ and £ (Gq 1 — iR) 1 . 
Therefore, that contribution is one-particle reducible. If we would start classifying 
all one-particle reducible and irreducible structures contributing to G accordingly, 
we would end up with the series (I2.80P with the self-energy £ containing all 1PI 
contributions. If this would not be the case, then we could always write G as 
the series in (I2.80j) plus a correction AG. However, since the sum of the series is 
known to be given by (I2.78[) this correction must be zero, i.e. AG = 0. 
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Most importantly, from the fact that E(0, G) contains only 1PI contributions 
one can conclude the important property of r 2 [0, G\ that it only contains two- 
particle irreducible (2PI) contributions with respect to the full propagator G. 
A contribution is said to be two-particle irreducible if it does not become 
disconnected by removing two inner propagators. Suppose Y 2 [0, G\ had a two- 
particle reducible contribution. The latter could be written as rCGT 7 , where 
GG denotes in a matrix notation two propagators connecting two parts P and 
T 7 of a contribution. Then E(0, G) would contain a contribution of the form 
rGT 7 since it is given by a derivative of P 2 with respect to G. Such a structure 
is one-particle reducible and cannot occur for the proper self-energy. Therefore, 
two-particle reducible contributions to r 2 [0, G\ have to be absent. 

This can be conveniently illustrated using a diagrammatic language, where 
propagators are associated to lines that can meet in the presence of interactions. 
For instance, for a quartic interaction this is exemplified for a two- and a three- 
loop graph contributing to T 2 below. Diagrammatically, the graphs contributing 
to £(0, G) are obtained by opening one propagator line in graphs contributing to 
r 2 [0, G], where the corresponding contributions are schematically given on the 
right: 



In these diagrams all lines are associated to the propagator G. Using the 
expansion (12.801) we see that each diagram corresponds to an infinite series of 
diagrams with lines associated to the “classical” propagator (G0 1 — iR ) _1 , and 
the above two- and three-loop diagrams contain, e.g., so-called “daisies” and 
“ladder” resummations: 




The diagrammatic loop expansion of the 2PI effective action will be discussed in 
detail in a subsequent chapter. 
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Details: One-loop effective action 


We derive the approximate result (12.69ft for the effective action r R [0] using a 
saddle-point approximation. Starting from the Legendre transform (12.65ft . we 
have 

r R [0] = -i\nZ[J,R} - [ J a {x)(f> a {x). (2.81) 

J X 

Using the defining functional integral (j2.62[) this can be written as 


e irR[(p] = J Shp exp i + J J a (x)(<p a (x) - 4> a (x)) 

= I S>Lp exp % \ S R [(j) + <p] + / J a {x)ip a (x) 


= e isR M / exp i <| S R [<f> + <p] - S R [0] + / J a (x)<p a (x) \ , (2.82) 


where for the second equality above we employed a field shift 

<Pa{x) ip a {x) + <f> a (x) . 


(2.83) 


Separating the classical part, we may write 

= s fl M + rf[0] (2.84) 

where the fluctuation part reads 


rf[0] = -iln I @ip expz <( S K [<j) + (f>\ - S K [(p] + / J a (x)ip a (x) 




r r 


(2.85) 


Furthermore, using the definitions (I2.68p and (12. 70j) one can verify that 


S R [cp + (p\- S R [cj)} = ]- [ (p a {x) [iG 0 ' ab (x, y■ </>) + R ab (x, y)\ <p b (y) 

^ J xy 

+ £' int [</?,0] + f ( p a {x ) ^ . (2.86) 

Jx 0(pa{x) 

Here we introduced the interaction term 

Sint[<P,<i>] = j Mx)Va(x)M X )Mx) - j {.Va(x)(p a (x)) 2 . (2.87) 

Starting from a classical action (I2.6ip with a quartic self-interaction, one observes 
from (12.871) that an effective cubic interaction term appears in the functional 
integral for the effective action in the presence of a non-zero field </>. 
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From plugging (12.86ft into (12.85ft we obtain with the help of the field equation 
(12.661) , which reads 

««(*]. ATM = (2.88, 


+ 


S<j>a{x) 8<t> a {x ) 

an exact functional integro-differential equation for Tf [0]: 


Tf[0] = —fin / @<p expi 


'xy 


*Sint \}P 5 0] 


' ( 

, Lf>a X S<t> a (x) J ’ 


<Pa(x) [iG 0 ^ ab (x, y\ 0) + R ab {x, y)] tp b (y) 

(2.89) 


The exponential in the integrand contains a term that is quadratic in the fields, 
a contribution 5j nt with cubic and quartic field powers as well as a linear source¬ 
like term ~ <5r^[0]/50. The linear term represents a so-called tadpole which 
guarantees a vanishing expectation value of the fluctuating field, {(p)j,R = 0, 
after the shift (12.831) . Its role is taken into account by taking 0 to describe the 
macroscopic field. If we neglect the cubic and quartic field powers contained in 
Sint, we arrive at the approximate saddle-point or one-loop expression 


rf ( hoop) [0] = -i In/^exp|-- J<p a (x) [G 0 ^(x, y, 0) - iR a b{x, y)] <p b {y) 

= —fin (det [G0 1 (0) — fit]) 2 + const 

(2.90) 

which leads to (12.691) up to an irrelevant constant. 


= - Tr In [G 0 x (0) — iR] + const, 


2.6 Nonequilibrium evolution equations 

2.6.1 Propagator evolution equations 

The functional representation of the nonequilibrium 2PI effective action, T[0, G], 
employs a one-point, 0, and a two-point function, G, whose physical values may 
be computed for all times of interest starting from the initial time t 0 . The 
equations of motion are given by the stationarity conditions (j2.73[) and (|2.74]h 
We will consider first the scalar field theory in the symmetric regime, where 
T[0 = 0, G] = T[G] is sufficient. We come back to the case of a nonvanishing field 
expectation value and consider also additional field degrees of freedom including 
fermionic and gauge fields below. 

The stationarity condition (j2.74[) for the propagator equation leads to (12.78p . 
which reads G~ b (x,y) = G^ ab (x,y) - E ab (x,y) - iR ab (x,y). Alternatively to 
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(12.801) . we may invert this equation by writing 


G ab (x,y) = G 0 , ab (x,y) + / G 0 , ac (x, z) [E cd (z,w) + iR cd (z,w)]G db (w,y). 

J zw,C 

(2.91) 

This can be directly verified with the help of (I2.78P by acting with G 1 on the 
left of (12.911) to get GG -1 = 1 and with Gf 1 — S — iR on the right of (I2.9ip to 
get in matrix notation: 

(Go + G 0 [£ + iR] G) (Gg 1 - S - iR) = G 0 Gg 1 - G 0 [S + iR] 

+ G 0 [S + iR] G (Gg 1 - E — iR) 

" -V-' 

= 1 G” 1 (2.92) 

Following section 12.41 the initial Gaussian density matrix leads to a non¬ 
zero source term R(x,y) at initial times x°,y° = to only- Therefore, (12.911) may 
be used to relate source terms to initial values of the propagator and its first 
derivatives for the relativistic theory (I2.6ip . However, typically the initial values 
for correlation functions are known and it would be favorable to prescribe them 
directly without recourse to sources. In particular, we may rewrite the equation 
of motion as a partial differential equation. For instance, starting from (I2.78P 
convolution with G using J_ c G“ c * (x, z)G cb (z, y) = S ab Sc(x — y) leads to 

G o%( x i z )G c b{z,y) ~ / [E ac (x, z) +iR ac (x,z)]G cb (z,y) = 5 ab d c {x - y). 

,C ’ Jz,C 

(2.93) 

In the absence of a macroscopic field we have 

V) = ~ [□* + ™ 2 ] fiabfic{x ~ V) (2-94) 

and one obtains the evolution equation for the propagator for 0 = 0: 




(2.95) 


which is an exact equation for known self-energy £. Since a source term 
representing the initial density matrix has support at initial time only, it gives 
a vanishing contribution to f vC R ac (x,z)G cb (z,y) for x° > t 0 . To solve this 
equation for some given approximation of E ab (x, y), the values for G ab (x,y ) and 
first derivatives have to be prescribed at x°, y° = t 0 . Then (I2.95P is used to evolve 
G ab (x,y ) for times larger than t 0 . We use this to describe in the following the 
evolution equation for physical propagators setting R = 0 in (I2.95P . 
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It is convenient to introduce a decomposition of the two-point function 
G into spectral and statistical components. The corresponding evolution 
equations for the spectral function and statistical propagator are fully equivalent 
to the evolution equation for G, but have a simple physical interpretation. 
While the spectral function encodes the spectrum of the theory, the statistical 
propagator gives information about occupation numbers. Loosely speaking, the 
decomposition makes explicit what states are available and how often they are 
occupied. This interpretation will become more clear as we proceed. 

For the neutral scalar held theory there are two linearly independent real¬ 
valued two-point functions, which may be associated to the real and the 
imaginary part of G. More precisely, we consider the expectation value of the 
commutator and the anti-commutator of two fields, 


commutator: p ab (x,y) = i([$ a (x),$ b (y)}), (2.96) 

anti-commutator (0 = 0): F ab (x, y) = -({<L a (z), $&(?/)}). (2.97) 

Here p ab (x,y ) denotes the spectral function and F ab (x,y ) the statistical two-point 
function. For neutral scalar fields the real functions obey F ab (x,y ) = F ba (y,x ) 
and p a b{x,y ) = —p ba (y,x). We note from (I2.96P that the spectral function p 
encodes the equal-time commutation relations 


Pab^x^ y^\ x o~y 0 0 , d x op ab fx^ y^\ x o = yO 5 ab 6(pt y) . (2.98) 


The spectral function is also directly related to the retarded or advanced 
propagator, G% b {x,y) = p a b{x,y)Q(x° - y °) = G^(y,x), respectively. 

The decomposition identity for spectral and statistical components of the 
propagator in the absence of sources reads: 


G a b{x, y) = F ab (x, y) - - p ab {x , y) sgn c (x° - y°) 


(2.99) 


Here sgn c (x° — y°) = 9c(x° — y°) — 9 c (y° — x°) using the contour step function 
introduced in (12.81) . With this the above decomposition is directly understood as 


Gat(x,y) ($a(x)$t{y))0 c (x° -y°) + (4j(»)$„(:r))0c(!/ 0 - x°) 

= f <{*«(*), * 6 (l0» («c(x° - V°) + »c(y° ~ x 0 )) 
!<[$„(*), *k(y)]> (»c(x° - y°) - 8c(.y° - *“)) 

2 V---v--' 

sgn c (x° -y°). 


Making the contour ordering explicit, we distinguish the four propagators 
G tb y) = F ab(x, y) — - pab(x, y) sgn(x° - y°) , 
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'I 

G ab~( x >y) = F ab(x, y) + -p ab (x,y) sgn(x° - y °), 

'l 

G tb (z, y) = F ab (x, y) + - p ab (x, y ) , 

G ab(x : y) = F ab (x,y) - ^p ab (x,y ) , (2.100) 

which correspond to the definitions (12.lip - (I2.12p with the standard sign 
function, sgn(a; 0 — y°) = 9(x° — y°) — 6(y° — x°). Since F ab (x, y) and p a b{x, y) are 
not time-ordered correlation functions, no distinction between C + and C~ for the 
location of their time arguments has to be taken into account here. 

To obtain a similar decomposition for the self-energy, we separate £ in a 
“local” and “nonlocal” part according to 

^ab(x,y) = -iT, ( $(x)8(x-y) + Y, ab (x,y). (2.101) 

Since E® just corresponds to a space-time dependent mass-shift it is convenient 
for the following to introduce the notation 

W S° ) m 2 S ab + S ( £(x ). (2.102) 

To make the time-ordering for the non-local part of the self-energy, Yl ab (x,y), 
explicit we can use the same identity as for the propagator (12.991) to decompose: 

Zab(x, y) = £ ab(x , y) - ^££ b (z, y) sgn c (x° - y°) . (2.103) 

Out of equilibrium we have to follow the time-evolution both for the statistical 
propagator, F, as well as for the spectral function, p. The evolution equations 
are obtained from (12.951) with the help of the identities (I2.99jl and (j2.103p . Most 
importantly, once expressed in terms of F and p, the time-ordering is explicit 
and the respective sign-functions appearing in the time-ordered propagator can 
be conveniently evaluated along the closed real-time contour C. 

With the notation (12.1011) the time evolution equation for the time-ordered 
propagator (12.951) reads 

[□Ac + Ml b (x)\ G cb (x,y)+ i / T, ac (x,z)G cb (z,y) = -i8 ab 5 c {x-y ), (2.104) 

J z,C 

where we have set R = 0. For the evaluation along the time contour C involved 
in the integration J[ = f c dx: 0 f d d z we employ (I2.99j) and (I2.103jl : 

i / Z ac {x,z)G cb {z,y ) = i / ^ c (x, z)F cb (z,y) 

J z )C J z 

~^ac(x, z)p cb {z , y) sgn c {z° - y°)-^ p ac (x, z)F cb (z , y) sgn c (x° - z°) 
~\Kc(x,z)p cb (z,y) sgn c (x° - ^°)sgn c (^ 0 - p°)| . (2.105) 
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The first term on the RHS vanishes because of integration along the closed time 
contour C. To proceed for the second term one splits the contour integral such 
that the sign-functions have a definite value, for instance 


dx; 0 sgn c (^° 


dz°(-l) + 


'to 




(2.106) 


for the closed contour with initial time to- We emphasize that the contributions 
from times later than y° simply cancel. To evaluate the last term on the RHS of 
(I2.105P it is convenient to distinguish the cases 


(a) 9 c (x° - y°) = 1: 


rto 


dz" sgn c (ar — z°)sgn c (z° — y°) = / dz°(-l)+ dz" + dz u (-l), 


'*0 


(2.107) 


(b) d c {y° — x°) = 1: 


d^ u sgn c (a; 0 - z°)sgn c (z u - y") = I dz u (-l) + / d z" + I dz u (-l) . 


v o „.o 


rto 


'to 


(2.108) 

One observes that (a) and (b) differ only by an overall sign factor sgn c (x°-y°). 
Combining the integrals therefore gives: 


i J T, ac (x,z)G cb (z,y) = j d d z 


d -° s ^c(x,z)F db (z,y ) 


'to 


dz°Za C (x,z)p cb (z,y)--sgn c (x u - y") I dz" Z p ac (x, z)p cb (z, y) [>(2.109) 


'to 


One finally employs 


C\ x G ab (x, y) = \3 x F ab (x, y) - - sgn c (a; 0 - y°)D x p ab (x, y) - i5 ab 5 c (x - y) (2.110) 

such that the (5-term cancels with the respective one on the RHS of the evolution 
equation (12.1041) . Here we have used 


- [Pab{x,y) sgn c (x° - y 0 )] = -^sgn c (x° - y°)d 2 x0 pab(x, y) 

—i5 c (x° - y°)d x op ab (x,y ) 

S -V-' 

- i5 ab 5c{x-y ) , 

where (I2.98P is employed for the last line and the term Pab{x, y)5 c (x° - y°) 
is observed to vanish since ~ p ab (x,y)\ x o =y o = 0. Comparing coefficients, which 
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here corresponds to separating real and imaginary parts, one finds from (12. 109[) 
and (I2.110P the equations for F ab (x,y ) and p ab (x,y). 

Using the abbreviated notation J'* 2 dz = J^ 2 d z° d d z we arrive at the 

coupled evolution equations for the statistical propagator and the spectral function: 



( 2 . 111 ) 


For the considered case with 0 = 0, the self-energies dependent only on F and 
p, i.e. M 2 = M 2 (F), £* = E F (p, F) and = Tj p (p,F). We note that the local 
self-energy correction (j2.102j) encoded in M 2 does not depend on the spectral 
function because the latter vanishes for equal-time arguments. 

One observes that (12.1111) are causal equations with characteristic “memory” 
integrals, which integrate over the time history of the evolution. We emphasize 
that the presence of memory integrals is a property of the exact theory and in 
accordance with all symmetries. In particular, the equations describe a unitary 
time evolution without further approximations. The equations themselves do 
not single out a direction of time and they should be clearly distinguished from 
phenomenological nonequilibrium equations, where irreversibility is typically put 
in by hand. Since the nonequilibrium evolution equations are exact for known 
self-energies, they are fully equivalent to any kind of identity for the two-point 
functions and are often called Kadanoff-Baym or Schwinger-Dyson equations. 

Note that the initial-time properties of the spectral function have to comply 
with the equal-time commutation relations (I2.98[) . In contrast, for F ab (x,y ) as 
well as its first derivatives the initial conditions at t 0 need to be supplied in order 
to solve these equations. To make contact with the discussion of initial conditions 
in section 12.41 we consider for a moment the spatially homogeneous case for 
which F ab (x,y) = F ab (x°, y°; x - y) = / d d p/(2ir) d exp[ip(x - y)]F ab (x°, y°] p) 
and equivalently for p ab (x,y). In terms of the Fourier components F ab (t,t'] p) 
the solution of the integro-differential equations (12 . 111 j) can be solved with the 
following initial conditions respecting the 0(N) symmetry of the theory: 

Fabif, t , p)| t=t'=to £p d a b j 

— (dfF ab (t, t , p) + dt'F ab {t, t , p))|£_£/_£ 0 £p??p dab 1 

d t d t 'F ab (t,t'- p)\ t=t/=to = (vl + ^jd ab . (2.112) 
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Here we have used that the required correlators at initial time are identical to 
those given in (12.461) for the considered case 0 = 0 if the momentum labels 
are attached for d > 0. Accordingly, these are the very same parameters 
that have to be specified for the corresponding Gaussian initial density matrix 
(12.431) . We emphasize that the initial conditions for the spectral function equation 
are completely fixed by the properties of the theory itself: the equal-time 
commutation relations (12.981) specify p a b(t,t p)\ t=t , =to = 0, d t p ab (t,t'] p)\ t=t , =to = 
S ab and d t d t 'p a b(t, p)\t=t'=t 0 — 0 for the anti-symmetric spectral function. 


2.6.2 Non-vanishing field expectation value 


In the presence of a non-zero field expectation value, 0^0, the same 
decompositions (12.99f) and (I2.103P apply. Equivalently, one can always view 
(I2.99|) as defining F and p from the connected propagator G. This means that 
the symmetric part, which is described by statistical two-point function, becomes 
in the presence of a macroscopic field 

F a b(x,y ) = ^({$ a (a;),$&( 2 /)}> - M x )Mv) , (2.113) 

whereas the anti-symmetric spectral function is still given by p ab {x,y ) = 
i([<I> a (a;), d*^?/)]). The general form of the scalar evolution equations for the 
spectral and statistical function (12.1 lift remains the same for 0 ^ 0. The only 
change compared to the symmetric regime is that the functional dependence 
now includes field dependent terms M 2 = M 2 (0, F), Ep = Ep(0, p, F) and 
£ p = S P (0,p, F ). 

For the N -component scalar field theory (12.61 jl with the field dependent 
inverse classical propagator (I2.70|) acting on the l.h.s. of the evolution equations 
for F ab (x,y) and p ab (x : y ), one has (0 2 = 0 a 0 a ): 

M ab(x) = (m 2 + 0 2 (x)^ 8 ab + ^ 0 a (x)0 6 (a;) + iff (a) • (2.114) 


In this case the evolution equations for the spectral function and statistical two- 
point function (j2. 111 j) are supplemented by a differential equation for 0 given 
by the stationarity condition (I2.73[) . For this we have to compute the functional 
derivative of (12.77 ]) . i.e. 


ST 


SS i <5{Tr G0 i (0)G} 


54>a(x) 8<j) a (x ) 


S<J>a(x) 


+ 


8 To 


84>a{x ) 


= ~Ja(x) 


(2.115) 


according to (I2.73j) for R = 0. Again, since a source term representing the initial 
density matrix has support only at time t 0 , we have J a (x) = 0 for x° > t 0 in the 
absence of external sources. With 


Tr Go X (0)G = -- 


a x + m 2 + ^ 0 2 ( X ) 


8 ab + 7^M x )M x ] 


F ba {x,X ) 

(2.116) 
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for the iV-component scalar field theory this yields the field evolution equation: 


| (tlx + m 2 + ^ [fi 2 (x) + F cc (x, x)] ^ S ab + ^ F ab (x, x) J <f> b (x) = ) 

(2.117) 

The solution of this equation requires specifying the field and its first derivative 
at initial time. In the context of the above discussion for spatially homogeneous 
fields this just corresponds to specifying fi a [t 0 ) and d t (p a (x 0 )\ x o =to . 

2.6.3 Evolution equations for fermions and gauge fields 

The different dynamical roles of spectral and statistical components is a generic 
property of nonequilibrium field theory and not specific to scalar field degrees 
of freedom. In terms of spectral and statistical components the equations for 
fermionic fields or gauge fields have very similar structures as well. To be specific 
we consider Dirac fermions T(x) and T = dCy 0 , in matrix notation with Dirac 
matrices 7 ^ for p — 0,..., 3^1 The time-ordered fermion propagator is A(x, y) = 
(T(x)T(y)) 6 l c (a : 0 — y°) — (ty (y)ty (x)) 9c(y° — x°), where the minus sign is a 
consequence of the anti-commuting property of fermionic fields. Correspondingly, 
in contrast to bosons, for fermions the field anti-commutator is associated to the 
spectral function^ 

anti-commutator: p^\x, y) — i({\P(x), ^(t/)}), (2.118) 

commutator: F^\x, y) = i([\I/(x), 0(7/)]) . (2.119) 

In terms of spectral and statistical components the propagator reads 

A (x,y) = F if \x,y) - l -p U \x, y) sgn c (x° - y°) . (2.120) 

The equal-time anti-commutation relations for the fields are again encoded in 
p(f\x,y). For instance, for Dirac fermions one has 

7 °p {f) (x,y)\ x o =y o =iS(x- y), ( 2 . 121 ) 

which will uniquely specify the initial conditions for the evolution equation for 
pd) similar to what is observed for bosons. 

For Dirac fermions with mass tttv' the free inverse propagator reads 

iA^ 1 (x,y) = \ip x - m U) ] 5{x - y ), (2.122) 

where j) = Similarly to the bosonic case, we may write in the absence of 

sources 

A ~ l (x,y) = A^ l (x,y) - E {/) (x,?/) , (2.123) 

3 Dirac matrices obey (yCy 1 "} = 2 g^ v with g^ v = diag(l, —1, —1, —1). 

4 We use (!') = (it) = 0 in the absence of sources. 
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which defines the proper fermion self-energy T,^\x,y). The corresponding 
decomposition for the fermion self-energy reads 

£ (/) 0 ,y) = ^ ( P(x,y) - *-Z ( p f \x,y)sgn c (x° - y°) . (2.124) 

If there is a local contribution to the proper self-energy, this is to be separated in 
complete analogy to the scalar equation (12.1011) and the decomposition (12.1241) is 
taken for the non-local part of the self-energy. 

From the equation of motion for the time-ordered fermion propagator (12.1231) 
we obtain a suitable time evolution equation by convoluting with A. For a free 
inverse propagator as in (I2.122P for Dirac fermions this yields 

[iftx ~ rn,( f) ] A(x,y) - i f S (/) (x, z)A(z, y) = i5 c (x - y) . (2.125) 

Jz,C 


Following the lines of the above discussion for scalars one finds for the fermion 
case the coupled evolution equations 


[ip x — m^] ( x , y) 

\}$x ~ p u \x,y) 


dzT, ( ' f \x,z)F {f \z,y) - I d z^\x,z)p if \z,y), 


'to 


I to 


dz T, ( ' f \x,z)p u \z,y) 


(2.126) 


Similarly, the nonequilibrium evolution equations for gauge fields can be obtained 
as well. For instance, for a gauge field A^(x) with Lorentz indices p = 0,..., 3 
the full time-ordered propagator D^(x,y) = (TcA fJ, (x)A l/ (y)) may be written as 


D^(x, y) = Fd{x, y) - 1 -Pd 0, y ) sgn c (x° - y°) . (2.127) 

For a theory with free inverse gauge field propagator given by 

iD o,lv( x >y)= - (l-C 1 )d tl d v \ x 5 (x-y) (2.128) 

for covariant gauges with gauge-fixing parameter £ and vanishing macroscopic 
field, (A^(x)) = 0, one finds the respective equations from (12 . 111 j) by 

[a, + M 1 }p(x,y) —> - [jyn- (1 - ryarl.Pot*,#) (2.129) 

and equivalently for F^ u (x,y). Of course, the respective Lorentz and internal 
indices have to be attached to the corresponding self-energies on the RHS of the 
equations. 
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Details: Functional integral for fermions 


In section 12.31 we derived a bosonic functional integral in terms of eigenvalues 
of the bosonic field operators. Fermionic operators anti-commute, which is 
intimately connected with the Pauli exclusion principle. In order to represent 
fermionic integrals, we will have to deal with anti-commuting eigenvalues of 
fermionic field operators. Anti-commuting numbers and functions are called 
Grassmann variables. Her we recall some basic aspects which we need for our 
purposes and refer to the literature for more extensive discussions. 

For a single anti-commuting variable 77 we have 


{v,v} = 0 , 


(2.130) 


such that 7777 = — 7777 , thus rf = 0. The derivative operator d/d 77 can be defined 
in an analogous way as for ordinary numbers and we write 


^ 3 _ 

d77 ’ 



= 1 . 


(2.131) 


As a consequence of (12. 13011 . the power series expansion of any suitable function 
f(r]) may be written as 

/ ( 77 ) = a + br] (2.132) 

with ordinary commuting numbers a and b. Therefore, d 2 f ()]) / dr] 2 = 0 and we 
have {d/d 77 , d/d 77 } = 0. 

In order to define integration for Grassmann variables, we start by asking 
what could be the analogue of the bosonic Gaussian integral, J dxexp{— ax 2 } = 
\J 7 T /a ? An important property is certainly the invariance under the shift of the 
integration variable x —> x + b. To carry over such a property to the fermionic 
integral, we define 

J dr] (a + brj) = b. (2.133) 

In fact, since r ] 2 = 0 the integral must be linear in 77 and the only linear function 
which respects the property of shift invariance is a constant. Consequently, 
integration for Grassmann variables acts just like differentiation. To write a 
Gaussian integral, we consider a complex anti-commuting variable 77 = (77# + 
i Vl )/y/2 where 77# and 77/ are real-valued anti-commuting variables. With 
V* = ( Vr ~ * ? ?7')/ v / 2 we find, for instance, 7777* = (77^ + 777/77# — 777#?7/ + r ] j )/2 = 
= The Gaussian integral is 


J dr]*dr] e~ b71 * v 


J dr]*dr] (1 — 677 * 77 ) 


= b 


dr]*dr] 7777 * 


b. 


(2.134) 
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where the sign change in the second equality conies from changing the order of 
the variables to be able to use first / dr/ rj = 1 and then J dy* y* = 1 . 

Next we want to consider the space of anti-commuting functions, 77(01). 
Generalizing (12.1311) . functional differentiation is defined as 

{^T^(y)} = S(x~y). (2.135) 

The corresponding Gaussian functional is 


S>y*2>y exp J y*(x)B(x,y)y(y ) j = det B , (2.136) 

where f < 3)y* l 3ly = U x dy*(x)dy(x). 

The construction of the 2PI effective action for fermionic holds proceeds along 
very similar lines as for bosons. However, one has to take into account the anti¬ 
commuting behavior of the fermion fields. A main differences compared to the 
bosonic case can be observed from the one-loop part (T2 = 0) of the 2PI effective 
action. Using again the Dirac fermions introduced in the previous section and, 
for comparison, a real scalar held <p with vanishing held expectation values it 
involves the integrals: 



/ — • (/) 

e lS ° = —i In (det Aq x ) = — iTrlnAg 1 , 
Bosons: —i In j QJip e lS ° — — i In (det Gq 1 ) 2 = ^Tr lnd// 1 . (2.137) 

Here S q = f d 4 xd 4 ytp(x)iAQ 1 (x, y)ip(y) denotes a fermion action that is bilinear 
in the Grassmann helds. For Dirac fermions the free inverse propagator A ( f 1 (x, y ) 
is given by (j2.122jl and ^ = 7)^7°• For the bosons S 0 is given by the quadratic 
part of (12.611) . Comparing the two integrals one observes that the factor 1/2 for 
the bosonic case is replaced by —1 for the fermion helds because of their anti- 
commuting property. With this difference, following along the lines as for the 
bosonic case, one hnds that the 2PI effective action for fermions can be written 
in complete analogy to (12.77f) . Accordingly, for the case of vanishing fermion held 
expectation values, (T) = (T) = 0, one has: 


T[A] = —i Trln A i Tr A 0 X A + r 2 [A] + const 


(2.138) 


Here r 2 [A] contains all 2PI diagrams with lines associated to the contour 
time ordered propagator A (x,y) = (T c T(a:)'k(7/)). The trace “Tr” includes 
integration over time and spatial coordinates, as well as summation over held 
indices. 
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As for the bosonic case of equation (12.741) . the equation of motion for A in 
absence of external sources is obtained by extremizing the effective action: 


^AL = 0 . 

5A(x : y) 

Using (12.1381) this stationarity condition can be written as 

= Aq x (x, y) - S (/) (x, y\ A) , 
with the proper fermion self-energy: 



(2.139) 


(2.140) 


(2.141) 


2.7 The special case of thermal equilibrium 

For the neutral scalar field theory the canonical thermal equilibrium density 
matrix is 

g 0 =^e-^ H (2.142) 

with the partition sum Z 0 = For some time t 0 we take g(t 0 ) = g 0 as 

the density matrix in the generating functional (12.ip . The construction of a path 
integral representation follows then closely section 12.21 The main difference is 
that e~P H = ma y be interpreted as an evolution operator in imaginary 

time. This simplifies taking into account a thermal density matrix which is 
not Gaussian in general. Following the steps of section 12.21 we evaluate the 
trace using eigenstates of the Heisenberg field operator <E>(x) at time x° = t 0 , 
<f>(£o,x)|<£>) = <£>o(x)|<£>). The partition function may then be written as 

Z 0 = j [dy?o] M Qp W) ■ (2.143) 

The interpretation of the canonical equilibrium density matrix as an evolution 
operator in imaginary time allows us to consider 

{(f\e~' m \ip) = (cp,t 0 -i/3\(p,t 0 ) (2.144) 

as a transition amplitude. As a consequence, we can use (I2.40[) applied to 
imaginary times in order to represent it as a path integral. Since we are 
calculating a trace, we have to identify 

ip(t 0 -i/3,x) = cp(t 0 ,x) = <p 0 (x). (2.145) 

This periodicity of the field in thermal equilibrium represents an important 
constraint which is, in general, not present out of equilibrium. 
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In order to write down a generating functional also for real-time correlation 
functions in thermal equilibrium, we extend with T c the time ordering to include 
the imaginary time branch from to to to — i/3. Here any time on the imaginary 
branch is considered to be later than any time on the previously discussed real¬ 
time contour C. Graphically, Cp is given below: 



We also include sources J{x) and R(x,y ) and use the notation f xC = 
J Cb dx° f d d x. The generating functional for contour time ordered correlation 
functions is then given by 

Zp[J,R] = f [d<^ 0 ] (tp\ gp T Cp expi j f $(x)J(x) 

J \Jx,Cp 

+ \f ^{x)R{x,y)^{y) i |^> (2.146) 

^ Jxyfip J 

Inserting complete sets of states along the contour following section [2721 we obtain 
Zp[J,R] = I @<pexp i<Sp[(f\+ J(x)<p(x) + ]- f R(x,y)<p(x)<p(y) > , 

J periodic I J x,Cp ^ Jxy,Cp 

(2.147) 

where the term periodic refers to the periodicity condition for the fields (j2.145j) . 
The thermal equilibrium contour action, Sp, for a theory with quartic self¬ 
interaction reads 

Sp[<p\ = j* ^ \]^^if{x)d^{x) - ^(p(x)(p(x) - ^ (cp(x)cp(x)) 2 1 . (2.148) 

For the specific case of calculating static quantities, the real-time extent 
of Cp can be taken to zero. Static quantities can then be obtained from the 
corresponding generating functional with the imaginary-time interval [0, — i(3] 
only, which is often formulated in “Euclidean” space-time with the introduction 
of a Euclidean time r = it and Euclidean action Se = —iS. Because of the 
periodicity condition (12.1 45 jh the Euclidean time dimension is compact. In 
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Fourier space this leads to a discrete set of so-called Matsubara frequencies, which 
we will not consider here. Quantities like the spectral function are, however, 
difficult to obtain from such a formulation since they would require analytic 
continuation. The above time contour including the real-time axis gives direct 
access to these quantities. We emphasize that in either formulation - including 
real times or not - all correlation functions are computed here with the thermal 
density matrix (12.142)1 such that the system is always in equilibrium. 

Contour time ordered correlation functions in thermal equilibrium can be 
obtained from the generating functional (12.147 ) ) by functional differentiation. For 
instance, the two-point function is given by 


1 8 2 Zp[J,R] 


Zp 8J(x)5J(y) 


j,r =o 


(T c ^(x)^(y)) = G^\x -y) + U- 


(2.149) 


For the last equality we have used that thermal equilibrium is translation 
invariant, such that the one-point function, </>, is homogeneous and, accordingly, 
the two-point function depends only on relative coordinates, G^ eq \x,y) = 
G (eq ) (x — y). Similar to f!2.99p . we may write 


G [eq \x - y) = F {eq \x -y)- 1 - p {eq \x - y) sgn C/3 (x° - y°) (2.150) 

with the contour sign function sgn C/3 (x° — y°) = d Cp (x° — y°) — 0 Cf3 (y° — x°). 

We now derive an important consequence of the periodicity condition (j2.145)1 
for correlation functions. At the initial time to of the contour Cp, we have 


(T Cj 4(i)4M)U.„ 0 = <J>(y)tM)U„„ 


and at the end point t Q — i/3, 


(2.151) 


(Tc /3 $(^) ®(y))\ x °= to -ip = {®(x)$(y)) \ x o =to -ip. (2.152) 


The periodicity condition for the fields in 


thermal equilibrium then implied 


($( 2 /)$(x))| a: o =to = ($(x)$( 2 /))| x o =to _ i/3 . (2.153) 

Using the decomposition (12.150)1 this reads 


( F ^ q \x-y)+ l -p^ q \x 



F^ eq \x-y) 



) x°=to—i/3 

(2.154) 


Translation invariance also makes it convenient to consider the Fourier transform 
with real four-momentum, (cn,p), 



d(jJ(¥ l p iu{x O_ y O )+ip(x _ y) ^(eq)/ 

(2tt)^ 1 ,P 


(2.155) 


5 


The periodicity condition is also called Kubo-Martin-Scliwinger (“KMS”) condition. 
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and equivalently for the spectral part. Taking for a moment for granted that these 
integrals can be properly regularized and defined for the considered quantum held 
theory, (12.1541) reads in Fourier space 


P (eq) (^P)) . (2.156) 

Solving for F ( eq ) ( uj , p) and using that (1 + e~ /3aj )/(I — e^ w ) = {eP w + 1)/(e^ w — 1) = 
1 + 2/{eP w — 1), the periodicity condition for bosons can be written as@ 


F (eq) (u;, p) + -p (eq) (w, p) = e~^ ( F (eq) (o;, p) - - 


= -iF + Mu)] P ) 


(2.157) 


with fpioj) = (& 6uj — 1 ) _1 denoting the Bose-Einstein distribution function. 
Equation (12.1571) relates the spectral function to the statistical propagator. This 
is also called fluctuation-dissipation relation, which will be discussed further in 
later chapters. While // eq ) encodes the information about the spectrum of the 
theory, one observes from (12.157jl that the function encodes the statistical 
aspects in terms of the occupation number distribution In the same way 

one obtains for the Fourier transforms of the spectral and statistical components 
of the self-energy the thermal equilibrium relation 

E^ q) (u;,p) = + s ? q) (wp)- (2.158) 

From (I2.157P and (12.158jl one obtains the identity 

E( eq )(o;, p)F( eq )(tn, p) - E { p\co, p)p (eq) (cn, p) = 0 . (2.159) 

Later we will see that this condition can be rewritten as the difference of a ’’gain” 
and ’’loss” term in kinetic descriptions in their range of applicability. 

For fermions the contour-ordered propagator reads in thermal equilibrium 
A( eq )(x — y) = (T(x)'L(|/)) 6 i c / 3 (x 0 — y °) — (4 fly)^>(x))6c p (y 0 — x° ) with the minus 
sign from the anti-commuting property of fermions. If we proceed along the lines 
as for bosons above, i.e. evaluating the contour-ordered correlation function once 
at the initial time and once at the final point of the contour, this minus sign leads 
to the anti-periodicity condition 

(V(y)'Sf(x)) | x o =to = -(tf(x)’® r (y))| I o =to _ i 0 . (2.160) 

Using the decomposition (12.120)1 for fermions, one can translate this again into a 
relation between the spectral and statistical two-point function. The difference 

6 In our conventions the Fourier transform of the real-valued anti-symmetric function p[x , y) 
is purely imaginary, while that of the symmetric function F{x,y ) is real. 
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is that the Bose-Einstein distribution in (I2.157P is replaced by the Fermi-Dirac 
distribution f^\u) = {eP^ + 1) _1 according to 1/2 + fp{oS) —> 1/2 — in 

the respective relation. 

It is important to realize that for a general out of equilibrium situation the 
spectral and statistical components of correlation functions are not related by 
a fluctuation-dissipation relation. Such a relation is a manifestation of the 
tremendous simplification that happens if the system can be characterized in 
terms of thermal equilibrium correlation functions. An even more stringent 
reduction occurs for the vacuum, where fp{uj) = 0 such that the spectral and 
statistical function agree up to a normalization. In this respect, nonequilibrium 
quantum held theory is more complicated, since it admits the description of more 
general situations and encompasses the thermal equilibrium or vacuum theory as 
special cases. 


2.8 Nonrelativistic quantum field theory limit 

For the nonrelativistic limit of a relativistic quantum held theory, one considers 
processes with typical momenta of particles that are small compared to their 
mass m. Consequently, we can approximately write \/p 2 + m 2 ~ m + p 2 /(2m). 
A related important aspect of the nonrelativistic limit is the absence of anti¬ 
particles, such that the typical chemical potential associated to the difference 
between the particle and anti-particle numbers is about /i ~ m. Therefore, in a 
nonrelativistic context it is often convenient to introduce the different chemical 
potential 

/i nr — n — m. (2.161) 

To be specihc, we consider a relativistic quantum theory for a complex scalar held 
ip(t, x). This theory can be directly related to the held theory for a = 1,..., N 
real held components ip a (t,x.) of section [2751 bv taking N = 2 and identifying 


<P 


1 

71 


(<Pi + i<P2) ■ 


(2.162) 


In the nonrelativistic regime, this complex scalar held theory can be effectively 
described in terms of the Gross-Pitaevskii held theory employed already in 
section PI for the physics of ultracold Bose gases. 

In order to relate the relativistic and nonrelativistic descriptions, it is 
convenient to start from the corresponding functional integral expression (I2.37|) . 
which represents the partition function before the conjugate momentum held is 
integrated out. Since the derivation in thermal equilibrium and out of equilibrium 
follows mostly along the same lines using the results of the previous sections - 
with the main difference that the time integration runs along the closed time 
path for the nonequilibrium system and includes the interval [0, —i/3] in thermal 
equilibrium - for the following we do not distinguish them in the notation. 
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This will be sufficient to derive, in particular, the nonrelativistic version of the 
relativistic evolution equations (12.1 lip . 

In terms of the complex field ip(t, x) and its conjugate momentum field 7 r(t, x), 
the functional integral for the partition function reads 


Z\p] = J StipStTr* £>7r exp <j i / dt ( / d d x [ 7 wd t <p + 7r*d t tp*} — [H — pQ] 

with the Hamiltonian 


H — I d 3 x ( 71*71 + (Vv?*)(V<y 2 ) + m 2 ip*(p + -^(ip*ip) 2 


(2.163) 

(2.164) 


and charge 

Q = i J d 3 x (tp*7 r* — Trip ). (2.165) 

Here we write the partition function without sources for brevity, which can always 
be added along the lines of section 12.21 

The term in the exponential of (12. 163j> is essentially the classical action for the 
interacting theory. We first consider the quadratic part of that action by setting 
A = 0 in (12 .164ft , which can be written in a compact matrix notation in spatial 
Fourier space as 


i dt 


d d p ( 


—1 


(2tt) 


p),-=Tr{t,p) 




u 


—l 0 id t + p 

id t + p — oj 


p) 

- 7 = 7 r*(t, p) 


(2.166) 

with c 0 = \J p 2 + m 2 . To proceed, we look for a canonical transformation such 
that the above quadratic part is diagonal in the new fields. This is achieved by 
introducing 


X = \/g V' 


1 * _ _ /w * 1 

sphf ' x 


(2.167) 


together with their complex conjugates. Inverted, these read ip — (x + X*)/v / 2cu 
and Ti = iyjuj/2(x* ~ x)- I 11 terms of the new fields, the quadratic part (12.166ft 
becomes a sum of a particle (+p) and an anti-particle (—p) contribution: 


i J dt J i^y 1 ^ x P ' 1 ^ ~ u + ^ + x *^’ p ) ^ idt ~ u ~ ^ p)} • 

(2.168) 

Using that cu ~ m + p 2 /(2m) and (I2.16ip . we can write for the particle 
contribution in configuration space 


i 


dt 


d d 


x 




V 2 

id t + ~- 1 " Put 

2m 



(2.169) 
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(2.170) 


Of course, taking into account the interaction in (12. 164|> for A 7^ 0, with 

^OV ) 2 = -^2 (x*x + ax + x*x* + x*x ) 2 

one observes that particles and anti-particles do not decouple in general. However, 
since in the nonrelativistic limit all characteristic scales are taken to be much 
smaller than the mass m with /1 — m <C m ~ /1, the term ~ x*[ m + lAx — ^ m X*X 
dominates the quadratic anti-particle part of the action. Loop corrections 
involving the “heavy” anti-particle modes are suppressed and the coupling terms 
involving anti-particle fields may be approximately neglected. In section [6] we 
show that this is an excellent approximation for the infrared scaling regime near 
nonthermal fixed points. 

With these approximations, the final result for the nonrelativistic limit of 
Z\/j] ~ Z nr [fi nr ] reads 


Zm[p> nr] = / W^Xexp < i dt d d x [ X 


V 2 

idt + X-1“ /^nr 

2m 


x - §(x*x ) 2 


Here the nonrelativistic coupling g is not dimensionless and given by 

A 


9 = 


24m 2 


(2.171) 

(2.172) 


To study the nonequilibrium evolution we can follow the developments of 
previous chapters, now starting from the nonrelativistic action for the complex 
field x( x ) with t = x° on the closed time path C, which reads for /x nr = 0: 




V 2 

iO x 0 T — 
2 m 



(2.173) 


Varying this classical action with respect to the field gives the Gross-Pitaevskii 
equation employed in (11.30[) . The second derivative gives the inverse classical 
propagator. For a complex field there are four combinations of the two derivatives 
with respect to x or X*- These can be efficiently described using the following 
index notation, 

XT Or) = x( x ) , X 2 (*r) = X*( x ) , (2.174) 

such that the inverse classical propagator matrix reads in terms of Xa f° r a = 1,2: 


^S[x,x*] 
Sxa(z)'5xt(y) 



-^+fxc 


+ 9X»Wx Sc(x-y). 

(2.175) 


Here = diag(l, —1) denotes the third Pauli matrix. 

Following the corresponding steps of section 12.61 we can write down the 
quantum evolution equations. For instance, the equations for the nonrelativistic 
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statistical two-point function F^ r \x, y) and the spectral function p^\x,y) read 
in the absence of a condensate!)] 

[ia 3 ac d x o - ft ac (x)] F c ( ° r) (x, y) = [ d^ E^ nr) (x, z)F^ r) (z, y) 

Jto 

- [ dz^ ni) (x,z)p < ff ) {z,y), (2.176) 

Jto 

/*X° 

[ia 3 ac d x o-n ac {x)\ p { ™\x,y) = / dz ^ p a { f r \x, z)p < ff\z, y), (2.177) 


where 




ab 


fAfe'r’b,*) + gF£\x,x) (2.178) 


and Efi” 1 (x, y) and T p( f' r> (x, y) denote the corresponding nonrelativistic 
statistical and spectral parts of the self-energies, respectively. 


2.9 Bibliography 

• 2PI or so-called ^-derivable approximation schemes have been put forward 
in J. M. Luttinger and J. C. Ward, Ground state energy of a many fermion 
system. 2, Phys. Rev. 118 (1960) 1417 and in G. Baym, Self consistent 
approximation in many body systems, Phys. Rev. 127 (1962) 1391, or 
J. M. Cornwall, R. Jackiw and E. Tomboulis, Effective Action For 
Composite Operators , Phys. Rev. D 10 (1974) 2428. 

• General discussions on 2PI or higher effective actions include H. Kleinert, 
Higher effective actions for bose systems, Fortschritte der Physik 30 (1982) 
187. A.N. Vasiliev, Functional Methods in Quantum Field Theory and 
Statistical Physics , Gordon and Breach Science Pub. (1998) and J. Berges, 
nPI effective action techniques for gauge theories, Phys. Rev. D 70 (2004) 
105010. 

• In the context of nonequilibrium dynamics, apart from L.P. Kadanoff, 
G. Baym, Quantum Statistical Mechanics, Benjamin, New York (1962), 
some reviews are K. c. Chou, Z. b. Su, B. 1. Hao and L. Yu, Equilibrium 
And Nonequilibrium Formalisms Made Unified, Phys. Rept. 118 (1985) 1, 
E. Calzetta and B. L. Hu, Nonequilibrium Quantum Fields: Closed Time 
Path Effective Action, Wigner Function and Boltzmann Equation, Phys. 
Rev. D 37 (1988) 2878, P. Danielewicz, Quantum Theory of Nonequilibrium 
Processes, /, Ann. Phys. 152 (1984) 239 and J. Berges, Introduction 

7 In full analogy to (12.961) and (12.971) . we can write in terms of the Heisenberg field operators 
Xa the two-point functions F^ T \x,y) = ({Xa(x) , xliv)}) / 2 and = i{[x a {x) ,x\{y)])- 


58 









to nonequilibrium quantum field theory, AIP Conf. Proc. 739 (2005) 3; 
arXiv: hep-ph/0409233. 

• An introduction to the close relation of these techniques with the 
nonequilibrium functional renormalization group is given in J. Berges 
and D. Mesterhazy, Introduction to the nonequilibrium functional 
renormalization group, Schladming lectures, arXiv: 1204.1489 [hep-ph], 

• Aspects of renormalization are discussed in H. van Hees and J. Knoll, 
Renormalization in self-consistent approximations schemes at finite 
temperature. I: Theory, Phys. Rev. D 65, 025010 (2002). J. P. Blaizot, 
E. Iancu and U. Reinosa, Renormalization of phi-derivable approximations 
in scalar field theories, Nucl. Phys. A 736, 149 (2004). J. Berges, 
S. Borsanyi, U. Reinosa and J. Serreau, Nonperturbative renormalization for 
2PI effective action techniques, Annals Phys. 320 (2005) 344. U. Reinosa 
and J. Serreau, 2PI effective action for gauge theories: Renormalization, 
JHEP 0607 (2006) 028. The renormalization of initial-value problems 
is specifically addressed in M. Garny and M. M. Muller, Kadanoff-Baym 
Equations with Non-Gaussian Initial Conditions: The Equilibrium Limit, 
Phys. Rev. D 80 (2009) 085011. 

• For a discussion of the relation between relativistic and nonrelativistic 
theories, which is available on the arXiv, see e.g. T.S. Evans, The condensed 
matter limit of relativistic QFT , in: Proc. 4th Workshop on Thermal Field 
Theories and their Applications, Dalian, China, 1995, Y.X. Guiet. et. al. 
(eds.), World Scientific, 1996; http://arxiv.org/abs/hep-ph/9510298. 


59 


3 Thermalization 


3.1 Two-particle irreducible loop or coupling expansion 

3.1.1 Expansion of the effective action 

Loop or coupling expansions of the 2PI effective action proceed along the same 
lines as the corresponding expansions for the standard 1PI effective action, with 
the only difference that 

• contributions are parametrized in terms of “dressed” propagators, which 
are obtained from the stationarity condition as in (I2.74j) . instead of classical 
propagators 

• and only 2PI contributions are kept. 

Using the decomposition (12.771) . the 2PI effective action contains the contributions 
from the classical action supplemented by a saddle-point (“one-loop”) correction 
and a 2PI contribution r 9 . For the example of the O (AQ-symmetric classical 
action (12.61 j) for the scalar field <f a (x), the contributions to T 2 [ 0 , G] are 
constructed from the effective interaction 

iS int [(/),ip\ = - f i^<j> a (x)<p a {x)(p b (x)v b (x) - f i^(<Pa{x)v a (x)) .(3.1) 

As discussed in section 12.51 this is obtained from the classical action by shifting 
(f a (x) —* <t> a (x)+(f a (x ) and collecting all terms cubic and quartic in the fluctuating 
held f a (x). As for the case of the 1PI effective action, in addition to the quartic 
interaction there is thus an effective cubic interaction for non-vanishing held 
expectation value. To obtain loop or coupling expansions, one Taylor expands 
(13. ip using e x = ^2 n x n /n\ such that only Gaussian functional integrals have to 
be performed to compute the corrections. 

To be specihc, we consider hrst the scalar held theory for a one-component 
held, i.e. N = 1 in ( 13. lj) . To lowest order one has T 2 [0, G\ = 0 and we recover 
the one-loop result given in (I2.75p . Further corrections can be very efficiently 
classihed with a standard graphical notation. The propagator G is associated to 
a line and the interactions are represented as points where four (three) lines meet 
for a four- (three-)vertex: 


G(x,y) = 1 




(3.2) 


Since T 2 [0, G] is a functional, which associates a number to the helds 0 and G, 
only diagrams with closed loops of propagator lines and vertices can contribute. 
We may classify the contributions to T 2 [0, G\ according to their number of closed 
loops, i.e. 

r 2 [ 0 , g\ = rf oop) [ 0 , g\ + rf oop) [ 0 , g\ + ... (3.3) 
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At two-loop order there are two contributions, r^ 21oop ) [<^ G\ — r^fG] + 
r^ 2b) [0,G]. The field independent contribution is 



(3.4) 


where we have made explicit the different factors coining from the overall —i in 
the defining functional integral for T[0, G] (see (12.621) h the factor 3 from the three 
different ways to form closed loops with a single four-vertex and the factor from 
the vertex itself. The field-dependent two-loop contribution is given by 


r^ 2b) [0,G] = f Q, 

Jxy,C 

(3.5) 

where the number of different ways to form closed loops from two three-vertices 
is 6 , and the factor 1/2 comes from the Taylor expansion of the exponential to 
second order. At three-loop order there are three contributions, r ^ 310013 " 1 [0, G\ = 

rf a) [G] + r? b) [0, g\ + [0, g], with 

r^ 3a) [G] = i^- [ g 4 

48 Jxyfi 

rfV.G] = f a 2 

* Jxyz,C 



-(3c) r 


>.G] = -i s 


G(x, y)G(y, z)G(z, x)G(x, w)G(w, y)G(w, z) 


' xyzw,C 





(3.6) 


Of course, care has to be taken when estimating the order of the coupling 
for a given contribution. For instance, one observes already at the classical level 
from 5S/S(j)(x ) = 0 for the simple case of a space-time independent field 0 that 
m 2 (f) + A0 3 /6 = 0 . As a consequence, for the case with spontaneous symmetry 
breaking with 0^0 and m 2 < 0 one finds 0 = ±-^/— 6 m 2 /A. Since parametrically 
0 ~ 1/VA, both of the above two-loop contributions should be taken into account 
to first order in the coupling A. Similarly, all three-loop contributions should be 
included at the next order. For general out-of-equilibrium situations the power 
counting can be much more involved even in the presence of a small coupling. We 
will consider explicit examples of time-dependent or “dynamical” power counting 
schemes for far-from-equilibrium dynamics in later chapters. 

The 2PI loop expansion exhibits of course much less topologically distinct 
diagrams than the respective 1PI expansion. For instance, in the symmetric phase 
(0 = 0 ) one finds that up to four loops only a single diagram contributes at each 
order. At fifth order there are two distinct diagrams which are shown along with 
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Figure 4: Topologically distinct diagrams in the 2PI loop expansion shown to 
five-loop order for 0 = 0. The suppressed prefactors are given in (13.8p . 


the lower-loop graphs in Fig. [4] suppressing prefactors. For later discussion, we 
give here the results up to five loops for the O (IV)-syrrimetric field theory with 
classical action (I2.6ip for general N. For 0 = 0 by virtue of O(N) rotations the 
propagator can be taken to be diagonal: 


G oi (i,y) = G(x,y)S ab . 


(3.7) 


This represents the most general form of the propagator for the symmetric regime 
if the initial conditions respect the O(N) symmetry. One finds to five-loop order 

with r 2 [G] = r 2 [0 = 0,G]: 


r f oop) [G]L _ G(5 

Z 1 i ' ( ~*ab= ( ~s d ab 

rf oop) [G]U _ G5 

Z L J I Gr a b — GrO a 5 


r (41oop) 

1 2 


[ G ]l G ab =GS ab 


rf oop) [G]L 

Z L -I I Gr a b — G rddl) 


A (IV + 2) 


i\ 2 (N + 2 ) 


G 2 (x, x ), 


' x.C 


G 4 (x,y ), 
48 3 N J xyfi 1 ,yh 

A 3 (N + 2)(iV + 8 ) 


48 


27N 2 


G (x,y)G ( x,z)G ( z,y ), 


' xyz,C 


iX 4 (N + 2)(N 2 + 6N + 20) 


128 811V 3 

x G 2 (z,w)G 2 (w,x) 


G 2 (x,y)G 2 (y,z ) 


J xyzw,C 

iX 4 (N + 2) (51V + 22) 


32 


811V 3 


G 2 (x,y ) 


' xyzw,C 


x G(x, z)G(x, w)G 2 (z, w)G(y, z)G(y, w ). 


(3.8) 


3.1.2 Self-energies 

In section 12.61 we have derived the coupled evolution equations (12.1111) for the 
statistical propagator F and the spectral function p as well as the field equation 
(12. 11 7j) . A systematic approximation to the exact equations can be obtained from 
the above loop or coupling expansion of the 2PI effective action. This determines 
all the required self-energies to a given order in the expansion from the variation 
of the 2PI effective action using (12.741) . We emphasize that all classifications of 
contributions are done for the effective action. Once an approximation order is 
specified on the level of the effective action, one obtains a closed set of equations 
and there are no further approximations on the level of the evolution equations 


62 



















required. This ensures the “conserving” properties of 2PI expansions such as 
energy conservation, since all approximate equations of motion are obtained from 
a variational principle. 

For the computation of self-energies with 

X(x,y) = -iW(x)5 c (x^y) + Z(x,y) = 2*^^ (3.9) 

oG[x,y) 

we again start with the one-component scalar field theory using the self-energy 
decomposition (I2.101|) . The field-independent two-loop contribution (j3.4j) gives 

r-p( 2 a)iv-Yj \ 

= -i^G(x,x)S c (x-y). (3.10) 

With the above self-energy decomposition we obtain 

E (0) = ^F{x,x), (3.11) 

where we have used that p(x, x) = 0. This contribution corresponds to a space- 

time dependent mass shift in the evolution equations. Similarly, for the field- 
dependent contribution one finds 

= “T ° 2{x ' ■ < 3 ' 12 > 

To write this in terms of statistical and spectral two-point functions, we use (I2.99j) 
and consider 

G 2 (x,y) = (^F(x, y) - ^sgn c (x° - y°)p(x,y)^J 

= F 2 (x,y) - ^p 2 {x,y) -isga c (x° - y°)F(x,y)p(x,y) , (3.13) 

where sgn^(x° — y°) = 1 has been employed. With the self-energy decomposition 
(I2.103P this leads to 

S F (x, y) = -y (V 2 (z, y) - ip 2 (x, y)^j (j>{x)(j>{y ), 

Z p (x,y) = -A 2 F(x,y)p(x,y)(p(x)(j)(y). (3.14) 


It remains to compute the loop contributions to the held evolution equation 


(I2TT7D . i.e. 


aifWG] 

54>(x) 




G 3 (x,y)cj)(y ). 


(3.15) 
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Similar to the above steps, we consider 


G 3 (x,y) = \F(x,y) - -sgn c (x° - y°)p(x,y) 

3i 

= F 3 (x, y) - — sgn c (x° - y°)F 2 {x , y)p(x, y) 

-\ F i x i y)p 2 ( x , y) + l s £ n c( x ° - y°)p 3 ( x , y) 


(3.16) 


using, in particular, sgn))(x° — y°) = sgn c (x° — y°). Following the discussion in 
section 12.61 we split the contour integration as 


d y°G 3 (x,y) = I dy u ( ~ F 2 (x, y)p(x, y) + -p 3 {x, y) 


+ 


'to 


' xo 


3 i 


d 2/° ( ^F 2 {x 1 y)p(x 1 y) - l -p 3 {x,y) 


= -3 i 


'to 


dy°p( x ,y) [F 2 (x,y)~— P 2 (x,y)\ . (3.17) 


12 ' 


Putting everything together gives 


«f b V,G] a 2 r" f 2 i 2 

= tt / d yp( x >y) F ( x ^y) - 7 fp fay) w > 


5(j){x) 


'to 


12 ' 


(3.18) 


where we use again the abbreviated notation J f 2 d z = f t 2 d^° f d d z. The 
evolution equations (12. 11 1)1 and (12. 1 17)1 together with the loop contributions 
(13.lip . (13 . 1 4)1 and (13. 18)1 represent a closed set of equations for F(x, y), p(x, y) and 
(f>{x). A corresponding discussion can be done starting from the above three-loop 
expressions for r 2 [0, G\ and, similarly, at higher loop orders. 

For later discussions, we present here also the case of the scalar O(N) 
symmetric held theory with a vanishing held expectation value starting from 
the three-loop F 2 . The propagator in the symmetric regime can be written as 
(13.7p . From F 2 [G]\G ab =G& ab , we obtain the self-energy as 

Z ab {x,y) = 2i 8V j^ ab=G5ab = 2& ^ 2 ^ Gab= ° Sab , (3.19) 

oG a b(x,y) N oG(x, y) 

where we have employed G(x,y ) = G a b(x,y)S a b/N for the last equality. Writing 
' F ab{x 1 y) = T l (x,y)5 a b and using the decomposition (12.10111 we thus have 

2 i fiT 2 [G}\G ab =G6 ab 


N 5G(x,y) 

From the two-loop contribution to T 2 [<j] given in (13. 8 p one hnds 

i(2)r 


(3.20) 


»™\c^a S . t = _ lX (l+A G(X)X)Sc(x _ y y (3. 21) 


N 5G(x,y ) 


6N 
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Using the propagator decomposition (12.991) in spectral and statistical components, 
this leads to the one-loop self-energy 


E ( 0 ) (x) 


A 


iV + 2 
6 N 


F(x,x ), 


(3.22) 


which corresponds to (13.lip for N — 1. Similarly, from Y 2 [G] at three-loop, one 
finds 


=-(21oop) , , >2^ 

s (X,y) = -A 


18 N 2 


G 3 (x,y). 


Using the decomposition (12. 103ft for the self-energy, we have 


Xf(x,v) 

Vp(x,y) 


-A 2 


IV + 2 
18 N 2 


F(x,y) 


F 2 (x,y) 


-A 2 


IV + 2 
6N 2 


p(x,y) 


F 2 (x, y) 


\p 2 (x,y) , 

T^p 2 (x,y) , 


(3.23) 


(3.24) 

(3.25) 


which enter (12.1021) and (12.1 lip . 


3.1.3 Solving nonequilibrium evolution equations 

The nonequilibrium evolution equations obtained from the loop expansion of 
the 2PI effective action can be used as a starting point for a wealth of further 
approximations that give analytical insight into the dynamics. For instance, we 
will later show that already two-loop self-energies as in (I3.24p and (I3.25P include 
and go beyond standard kinetic or Boltzmann equations describing two-to-two 
scattering of particles. However, it is very instructive to consider their numerical 
solution without further approximations. In particular, this will allow us to 
observe the striking process of thermalization in quantum field theory from first 
principles. 

Despite the apparent complexity of the evolution equations (12.1 lip and (12.117P 
they can be very efficiently implemented and solved on a computer. Here it is 
important to note that all equations are explicit in time, i.e. all quantities at 
some time t can be obtained by integration over the explicitly known functions 
for earlier times t past < t for given initial conditions. The time-evolution 
equations (12.1111) for p(x,y)\ x o =ti y o =t2 and F(x,y)\ x o =tljV o =t2 do not depend on 
the RHS on p(x,y ) and F(x,y) for times x° > t\ and y° > t 2 and similarly for 
the field evolution equation (I2.117p . To see this, we note that the integrands 
vanish identically for the upper time-limits of the memory integrals because 
of the anti-symmetry of the spectral components, with p(x,y)\ x o =y o = 0 and 
Ti p (x,y)\ x o =y o =0. As a consequence, only explicitly known quantities at earlier 
times determine the time evolution of the unknowns at later times. The numerical 
implementation, therefore, only involves sums over known functions. 

Solving the evolution equations requires the specification of (Gaussian) initial 
conditions as discussed in section 12.41 We consider first the one-component 
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scalar field theory for spatially homogeneous systems and Fourier transform with 
respect to the spatial coordinates. The equal-time commutation relations fix the 
Fourier modes of the spectral function p(x°, x°, p) = 0, d x op(x°, y°] p) = 1 and 
d x od y op(x°,x°] p) = 0 also at initial time. Here we take 

F(x°,y°-,p)\ x o =yO=0 = ~(/p(°) + 0’ 

d x od y oF(x°,y°;p)\ x o =y o =0 = uj p ^/ p (0) + 0 , 

d x oF(x°,y°;p)\ x o =yO=0 = 0, 0(^ O )U°=o = 0, d x o(j)(x°)\ x o =0 (3.26) 


with ca p = y/p 2 + m 2 as initial conditions. Since the initial macroscopic field and 
its derivative are taken to be zero, this remains so at all times for the considered Z 2 
symmetric theory with the corresponding invariance of the action under (j) —> — <fi. 
Here / p (0) plays the role of an initial occupation number distribution which, for 
instance, may be taken as 


/p(0) = A/"exp (-t^GpI - IPtsD 2 ^) • (3.27) 

The real parameter cr controls the width of the initial distribution and Af is 
a normalization constant. The initial condition is clearly far from thermal 
equilibrium and reminiscent of two colliding high energy wave packets. More 
precisely, it describes a spatially homogeneous collection of particles which move 
with approximately the same momentum peaked around p ts and — p ts . Of course, 
other initial distributions can be considered which we will also do in the following. 

In addition, we have to state renormalization conditions in order to define the 
quantum field theory as for any renormalizable theory in vacuum or thermal 
equilibrium. From this perspective, it is rather simple to consider first a 
field theory in one spatial dimension. In this case, the one-loop correction 
(13. lip evaluated with a momentum regulator A shows a logarithmic regulator 
dependence. This can be analyzed at initial time, where for the spatially 
homogeneous initial conditions (I3.26[) we find for / p (0) = 0 that the large 
momentum behavior for d = 1 is given by 


E ( 0 ) (f 0 ) 


x r A dp 1 

2 J 2^2^p 


~ In A. 


(3.28) 


We absorb this regulator dependence in a bare mass parameter p with the 
replacement m 2 t —y p 2 = m 2 — dm 2 in (12. 102)1 . The counterterm 5m 2 cancels 
the divergent vacuum contribution coming from the one-loop graph. The finite 
part of 5m 2 is fixed by requiring that the renormalized one-loop mass parameter 
in vacuum (/ p (0) = 0) equals m and we express all dimensionful scales in units 
of m. Since higher loop corrections are finite for one spatial dimension as A is 
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Figure 5: Evolution of the equal-time two-point function F(t,t] p) with Fourier 
modes |p| = 0,3,5 from the 2PI three-loop effective action. The evolution is 
shown for three very different nonequilibrium initial conditions with the same 
energy density. (All in units of m.) 


sent to infinity, the renormalization procedure is straightforward in this casejf] A 
crucial question of nonequilibrium dynamics is how quickly the system effectively 
looses the details about the initial conditions. Thermal equilibrium keeps no 
memory about the time history except for the values of a few conserved charges. 
As a consequence, if the real scalar field theory approaches thermal equilibrium 
then the late-time result would be uniquely determined by energy density. 

As a first example, we show in Fig. [5] the time dependence of the equal¬ 
time propagator F(t, f; p) for three Fourier modes |p|/m = 0,3,5 and three 
very different initial conditions with the same energy density. The coupling 
is chosen for all runs to be A/m 2 = 10. The large value suggests important 
corrections beyond the two-loop self-energy approximation, however, for the 
moment it allows us to discuss the relevant qualitative properties in a pronounced 
way and we will later consider the question of convergence in more detail. For 
the solid line in Fig. ED the initial conditions are close to a mean held thermal 
solution with inverse temperature /3 = 0.1/m, the initial mode distribution for 
the dashed and the dashed-dotted lines deviate more and more substantially 
from a thermal equilibrium distribution. One observes that propagator modes 
with very different initial values but with the same momentum |p| approach 
the same large-time value. The asymptotic behavior of the two-point function 
modes is uniquely determined by the initial energy density. Concerning the 
underlying nonequilibrium processes, we will see that this is also a very interesting 
observation in view of strong kinematic restrictions for scattering processes in one 

8 This is different for higher spatial dimensions, where a consistent renormalization 
will require taking into account non-Gaussian corrections. For further reading on the 
renormalization of 2PI effective actions, see the literature in section [2751 
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spacial dimension. 

The shown results are obtained from a rather simple discretization of the 
equations for the Fourier components F(t,t'; p) and p(t,t'] p). They are based 
on a time discretization t = not , t' = mat with stepsize at such that F(t,t') t—>■ 
F(n, m), and 


d^F(t,t') i ^ (f(ti + 1 , m) + F(n — 1 , m) — 2 F(n, m)V 
a t \ / 

• n —1 

dtF(t,t') n- q t (V(0, m)/2 + ^F(Z, m) + F(n,m)/2) , 


(3.29) 

(3.30) 


where we have suppressed the momentum labels in the notation. Above, the 
second derivative is replaced by a finite-difference expression which is symmetric 
in at -H- —at- It is obtained from employing subsequently the lattice “forward 
derivative” [F(n + l,m) — F(n,m)\/at and “backward derivative” [F(n,m) — 
F(n — 1 ,m)\/a t . The integral is approximated using the trapezoidal rule with 
the average function value [F(n, m) + F(n + 1, m )]/2 in an interval of length a t . 
The time-discretized version of (12.1111) reads then for the one-component theory 


F{n + 1, m; p) = 2 F(n, m; p) — F[n — 1, m; p) 

— a^ |p 2 + m 2 + ^ J F(n,n ; k)| f(ti, m; p) 

-a 3 t |s p (n, 0; p) F( 0, m; p)/2 - E F (n, 0; p) p(0, m; p)/2 (3.31) 

m—1 

+ ^ P) "b P) “ ^ P) /°(*» P)) 

z=i 

n—1 

+ 5^ S p (n,Z;p)F(Z,m;p) 

l=m 

and similarly for the spectral function. These equations are explicit in time: 
starting with n — 1 , for the time step n + 1 one computes successively all 
entries with m = 0,..., n + 1 from known functions at earlier times. The 
above discretization leads already to stable numerics for small enough stepsize a t , 
but the convergence properties may be easily improved with more sophisticated 
standard estimators if required. 

As for the continuum, the propagators obey the symmetry properties 
F(n,m ) = F(rn,n) and p(n,m ) = —p(m,n). Consequently, only “half” of 
the (n, m)-matrices have to be computed and p(n,n) = 0. Similarly, since the 
self-energy E p is antisymmetric in time one can exploit that E p (n,n) vanishes 
identically. As initial conditions one has to specify F(0,0 ;p), F(l,0;p) and 




F(l,l;p), while p(0,0;p), p(l,0;p) and p(l, l;p) are fixed by the equal-time 
commutation relations (12.98j) . 

It is crucial for an efficient numerical implementation that each step forward 
in time does not involve the solution of a self-consistent or gap equation. This 
is manifest in the above discretization. The main numerical limitation of the 
approach is set by the time integrals (“memory integrals”) which grow with 
time and therefore slow down the numerical evaluation. Typically, the influence 
of early times on the late time behavior is suppressed and can be neglected 
numerically in a controlled way. In this case, it is often sufficient to only take 
into account the contributions from the memory integrals for times much larger 
then a characteristic inverse damping rate. An error estimate then involves a 
series of runs with increasing memory time. 

For scalars one may use for the spacial dependence a standard lattice 
discretization with periodic boundary conditions. For a spatial volume V = 
(N s a s ) d in d dimensions with lattice spacing a s one finds for the momenta 



(3.32) 


where n* = 0,..., N s — 1. This can be easily understood from acting with the 
corresponding finite-difference expression (13.29jl for space-components: d 2 e~ ipx K > 
e -ipx^ e i P a s _|_ e -ipa s _ 2 y a 2 s = —e~ ipx As\\\ 2 {pa s /2)/a 2 . Exploiting the lattice 
symmetries reduces the number of independent lattice sites. Taken this explicitly 
into account becomes the more important the larger the space dimension. For 
instance, on the lattice there is only a subgroup of the rotation symmetry 
generated by the permutations of p x ,P y ,Pz and the reflections p x -B- — p x etc. for 
d = 3. Exploiting these lattice symmetries reduces the number of independent 
lattice sites to ( N s + 1)(N S + 3)(IV S + 5)/48. The self-energies may be calculated 
in coordinate space, where they are given by products of coordinate-space 
correlation functions, and then transformed back to momentum space. The 
coordinate-space correlation functions are available by fast Fourier transformation 
routines. 

The lattice introduces a momentum cutoff n/a s , however, the renormalized 
quantities should be insensitive to cutoff variations for sufficiently large n/a s . 
In order to study the infinite volume limit one has to remove finite size effects. 
This may be done by increasing the volume until convergence of the results is 
observed. For time evolution problems the volume which is necessary to reach the 
infinite volume limit to a given accuracy can depend on the time scale. This is, in 
particular, due to the fact that finite systems can show characteristic recurrence 
times after which an initial effective damping of oscillations can be reversed. 
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3.1.4 Nonequilibrium evolution of the spectral function 


In order to understand the nonequilibrium dynamics of the spectral function p 
vs. the statistical two-point function F, it is helpful to consider for a moment the 
free field theory. With E = 0 in the nonequilibrium evolution equations (12. 11 ID . 
the above initial conditions (13.26[) lead to the plane wave solutions 

F <f, “V,i/ 0 ; p) = 2- (/ p (0) + i) cos K(x“ - /)] , 

= — s in [w p (x° - y 0 )] (3.33) 

(Up 

with frequency cu p . We may bring this into a more suggestive form by introducing 
the center coordinate A 0 = (x° + y 0 )/ 2 and the relative coordinate s° = x° — y°, 
with respect to which a Fourier transformation to momentum space is performed 
(Wigner transformation). To analyze the spectral function we may perform a 
Wigner transformation and write 

p2X° 

ip(X°;u, p) = / ds°e iu}S °p{X° + s°/2,X° - s°/2;p). (3.34) 

J—2X° 


The i is introduced such that p(X °; cu, p) is real. Since we consider an initial-value 
problem with x°,y° > 0, the time integral over s° is bounded by its maximum 
value for y° — 0 where s° = 2A° and its minimum value for x° = 0 where 
s° = —2A°. After performing a Wigner transformation of the free spectral 
function ('13.331) . we find 


P 


(free) (A 0 ;(u,P) = 


sin [(cu — (Up)2W°] sin [(a; + cu p ) 2X°] 


(Up((U (Up) 


(Up ((U + (Up) 


(3.35) 


For hnite A" 0 this spectral function shows a rapidly oscillating behavior, while its 
envelope is peaked at cu = ±cu p . In the limit A 0 —* oo it reduces to 

p (free) ((u,p) = 27rsgn((u) <5 (cu 2 — (Up) , (3.36) 


which describes the familiar form of the free spectral function in Fourier space. 
In general, the positivity condition sgn(cu)p(A°; cu, p) > 0 can only be shown 
to hold in the special case that the initial density matrix commutes with the 
full Hamiltonian, such as in thermal equilibrium. In this case the system is of 
course stationary and independent of A 0 . As a consequence, the interpretation 
of the Wigner transformed nonequilibrium spectral function as the density 
of states should be taken with care. Nevertheless, as a consequence of the 
equal-time commutation relation the Wigner transform obeys the sum rule 
J dcu/(27r) cup(A°; cu, p) = 1 for the free as well as the interacting theory. 
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Figure 6: Wigner transform of the spectral function as a function of u at 
A" 0 = 35.1 for p = 0 in units of m. Also shown is a fit to a Breit-Wigner 
function (dotted) with (u; 0 ,r 0 ) = (1.46,0.37). The inset shows a blow-up around 
the three-particle threshold 3cuo- The expected bump is small but visible. Here 
the coupling is A/m 2 = 4. 


Going beyond the free theory, we know that the equilibrium spectral function 
would acquire a “width 1 ' in the presence of a non-zero imaginary part of the 
self-energy, i.e. £ p ^ o|H Such a width in Fourier space is related to a damping 
of correlation functions in time. To obtain a qualitative understanding of this 
statement, we may consider the case of some spectral function which is assumed 
to be strictly exponentially damped, i.e. p(x°, y° ; p) = e _7p l x ~ y ^E ~ 1 sin[F , p (a: 0 — 
?/ 0 )]. The effective frequency E p and rate q p are allowed to depend on the time 
A" 0 . The corresponding Wigner transform reads p(A°;o;,p) = pBw(Af°; u;, p) + 
5p(X°] uj, p), where psw denotes the Breit-Wigner function 


Pbw(-A°;o;, p) 


2cur p (A°) 

[cu 2 - E 2 (X 0 )} 2 + u; 2 T 2 p (X 0 ) 


(3.37) 


with a width r p (A°) = 27 P (A°). The additional contribution 5p(A°;o;,p) 
vanishes exponentially as exp(—r p A°), such that the finite-time effect is indeed 
suppressed for A 0 A l/r p . 

For the numerical solution we solve the evolution equations (12.1 lip using the 
above two-loop self-energy contributions. The initial conditions are given by 
(I3.26P with the initial distribution fl3.27p . For the plots we have used a space 
lattice with spacing ma = 0.3 and a time lattice a^/a = 0.25. The system size 
is mL = 24. In Fig. Owe display the Wigner transform p(A°;cu,p) for the zero 
momentum mode for mX° = 35.1. One clearly observes that the interacting 

9 For a vanishing w-dependence the function l4 eq )(w, p) = E/ eq ^(w, p)/2 w plays the role of a 
decay rate for one-particle excited states with momentum p. 
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theory has a continuous spectrum described by a peaked spectral function with 
a nonzero width. The inset shows a blow-up of the zero mode around the three- 
particle threshold 3uo/m = 4.38. The expected enhancement in the spectral 
function is small but visible. In the figure we also present a fit to a Breit- 
Wigner spectral function. While the position of the peak can be fitted easily, 
the overall shape and width are only qualitatively captured. In particular, the 
slope of p(X°]co,p) for small u is quantitatively different. We also see that the 
Breit-Wigner fits give a narrower spectral function (smaller width) and therefore 
would predict a somewhat slower exponential relaxation in real time. 

3.2 Two-particle irreducible 1/N expansion 

3.2.1 Power counting 

The above loop or coupling expansion is restricted, a priori, to weakly coupled 
systems. However, even in the presence of a weak coupling the dynamics can be 
strongly correlated. Important examples concern systems with high occupation 
numbers, where the statistical correlation function F grows large. In this section 
we discuss a nonperturbative approximation scheme for the 2PI effective action. 
It classifies the contributions to the 2PI effective action according to their scaling 
with powers of 1/N, where N denotes the number of held components: 

r 2 [0,G] = r^°[0,G] + rf>,G] + t™ lo [0,g] + ... 

~ iV 1 ~ N° ~ TV -1 

Each subsequent contribution T^ 0 , r^ LO , r^ NLO etc. is down by an additional 
factor of 1/N. The importance of an expansion in powers of 1/iV stems from 
the fact that it provides a controlled expansion parameter that is not based 
on weak couplings. It can be applied to describe physics characterized by 
nonperturbatively large fluctuations, such as encountered near second-order phase 
transitions in thermal equilibrium, or for extreme nonequilibrium phenomena 
such as parametric resonance and the subsequent emergence of strong turbulence 
just to mention two examples which are discussed in later chapters. For the latter 
cases a 2PI coupling or loop expansion is not applicable. The method can be 
applied to bosonic or fermionic theories alike if a suitable held number parameter 
is available, and we exemplify it here for the case of the scalar O(N )-symmetric 
theory with classical action (I2.61|) . 

We present a classihcation scheme based on O (iV)-invariants which 
parametrize the 2PI diagrams contributing to T[0, G\. The interaction term of 
the classical action in (12.6ip is written such that S[(p] scales proportional to N. 
From the helds (j) a alone one can construct only one independent invariant under 
O(N) rotations, which can be taken as tr (jxf> = 0 2 = 0 a 0 a ~ N. The minimum 
0o of the classical effective potential for this theory is given by (/>/ = N(—6 m 2 /A) 
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Figure 7: Graphical representation of the ^-dependent contributions for T 2 = 0. 
The crosses denote held insertions ~ 0 a 0 a for the left figure, which contributes 
at leading order, and ~ 0 a 0;, for the right figure contributing at next-to-leading 
order. 

for negative mass-squared m 2 and scales proportional to N. Similarly, the trace 
with respect to the held indices of the classical propagator Go is of order N. 

The 2PI effective action is a singlet under O(N) rotations and parametrized 
by the two helds 0 a and G a &. To write down the possible O(N) invariants, which 
can be constructed from these helds, we note that the number of 0-helds has to 
be even in order to construct an O(N) -singlet. For a compact notation we use 
(4>4>)ab = 0a06- All functions of 0 and G, which are singlets under O(N), can be 
built from the irreducible (i.e. uonfactorizable in field-index space) invariants 

0 2 , tr (G n ) and tr(00G n ). (3.38) 

We note that for given N only the invariants with n < N are irreducible — there 
cannot be more independent invariants than helds. We will see below that for 
lower orders in the 1/N expansion and for sufficiently large N one has n < N. 
In particular, for the next-to-leading order approximation one finds that only 
invariants with n < 2 appear, which makes the expansion scheme appealing from 
a practical point of view. 

Since each single graph contributing to T[0, G] is an 0(iV)-singlet, we can 
express them with the help of the set of invariants in (13.381) . The factors of N in 
a given graph have two origins: 

• each irreducible invariant is taken to scale proportional to N since it 
contains exactly one trace over the held indices, 

• while each vertex provides a factor of 1/iV. 

The expression (12.771) for the 2PI effective action contains, besides the classical 
action, the one-loop contribution proportional to Tr InG -1 + Tr G0 x (0)G and a 
nonvanishing T2[0, G] if higher loops are taken into account. The one-loop term 
contains both leading order (LO) and next-to-leading order (NLO) contributions. 
The logarithmic term corresponds, in absence of other terms, simply to the free 
held effective action and scales proportional to the number of held components 
N. To separate the LO and NLO contributions at the one-loop level consider 
the second term Tr G0 (0)G. From the form of the classical propagator (12.701) 
one observes that it can be decomposed into a term proportional to tr (G) ~ N 
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and terms ~ (A/6 N) [tr (00) tr (G) + 2 tr (00G)]. This can be seen as the sum of 
two “2PI one-loop graphs” with held insertion ~ 0 a 0a and ~ 0 a 06, respectively, 
as shown in Fig. [71 Counting the factors of N coming from the traces and the 
prefactor, one sees that only the first contributes at LO, while the second one is 
NLO. 


According to the above rules one draws all topologically distinct 2P1 diagrams 
and counts the number of closed lines as well as the number of lines connecting two 
held insertions in a diagram following the indices. For instance, the left diagram 
of Fig. [7] admits one line connecting two held insertions and one closed line. In 
contrast, the right hgure admits only one line connecting two held insertions by 
following the indices. Therefore, the right graph exhibits one factor of N less and 
becomes subleading. Similarly, for the two-loop graph below one finds: 

a a a b a b 




The same can be applied to higher orders. We consider hrst the contributions to 
r 2 [0 = 0,G] = r 2 [G], i.e. for a vanishing held expectation value and discuss 0 ^ 0 
below. The LO contribution to F 2 [G] consists of only one two-loop graph, whereas 
to NLO there is an infinite series of contributions which can be analytically 
summed: 


r 2 °[G] 

= 4 !jv/ G aa(x,x)G bb (x,x) , 

(3.39) 

r^ LO [o] 

= ^Tr ln[ B(G) ], 

(3.40) 

B{x,y;G) 

= 5(x~y) + i-^G ab (x,y)G ab (x,y). 

(3.41) 


In order to see that (I3.40p with (I3.4ip sums the following infinite series of diagrams 



one can expand: 


Tr ln[ B(G)) 


i-^Gah( x , x )G a b( x , x ) 


\\ x {^^ Gab ^ X ' V ^ Gab ^ X ' y ^ 


i-^ G a'b'( y , x ) G a'b'( y , x ) 


(3.42) 
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The first term on the RHS of (13.421) corresponds to the two-loop graph with 
the index structure exhibiting one trace such that the contribution scales as 
tr G 2 /N ~ N°. One observes that each additional contribution scales as well 
proportional to (tr G 2 /N) n ~ iV° for all n > 2. Thus all terms contribute at the 
same order. 

The terms appearing in the presence of a nonvanishing field expectation value 
are obtained from the effectively cubic interaction term in (13. ip for 0 ^ 0. One 
first notes that there is no ^-dependent graph contributing at LO. To NLO there 
is again an infinite series of diagrams ~ N° which can be summed: 

r!r>,G] = r“[G], (3.43) 

r* LO [0,G] = r^ LO [0 = o,G] + ^y I(x, y] G)Mx)G ab (x,y)My)X3-M) 

X X f 

I{x,V,G) = 7 rTjG ai {x,y)G^(x,y)~i — / I(x, z;G)G ai (z,y)G ab (z,y). 

(3.45) 

The series of terms contained in (I3.44p with (I3.45P corresponds to the diagrams: 



The functions I(x,y ; G) and the inverse of B(x,y, G ) are closely related by 

B~\x, y ; G) = S(x - y) - il{x, y, G) , (3.46) 

which follows from convoluting (13.41 j) with B and using (I3.45p . We note that 
B and / do not depend on 0, and T 2 [0, G\ is only quadratic in 0 at NLO. 

It is straightforward to apply the above description of an expansion of the 
2P1 effective action in the number of field components to boson or fermion 
field theories with “vector-like” TV-component fields. It is very helpful that for 
the TV-component theory discussed above one can analytically sum the infinite 
series of NLO contributions analytically. The situation is different for theories 
with “matrix” fields, such SU(N) gauge theories relevant for QCD, where no 
corresponding closed expression is known. 

3.2.2 Self-energies in the symmetric regime 

We consider first the case of the scalar O(N) symmetric field theory with a 
vanishing field expectation value such that F a b(x,y ) = F(x,y)6 a b and p a b(x,y ) = 
p(x,y)5 a b. The case 0 7^ 0 is treated below. In the 1/iV-expansion of the 
2P1 effective action to next-to-leading order the effective mass term M 2 (x]G) 
appearing in the evolution equations (12.1111) is given by 

M 2 (x-F ) = m 2 + F(x,x). (3.47) 

6 N 
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One observes that this local self-energy part receives LO and NLO contributions. 
In contrast, the non-local part of the self-energy (12.10111 is nonvanishing only 
at NLO: T,(x,y,G) = —\/(3N)G(x,y)I(x,y) and using the decomposition 
identities (j2.99j) and (I2.103P one finds 


Xf(x,v) = 


Z p (x,y) = 


( F (x, y) - y)I p {x, y ) ) , 


A 


3N 


F(x, y)I p (x, y) + p(x, y)I F (x, y) 


(3.48) 


(3.49) 


Here the summation function (13.45ft reads in terms of its statistical and spectral 
components!^ 


I F (x,y) = n F (x,y)~ / dz I p (x, z)U F (z,y) + / dz I F (x, z)U p (z,y), 


to 


to 


I p {x,y) = U p (x,y)~ dz Ip(x,z)Up(z,y), 


(3.50) 


where 


H F (x,y) 
n p (x,y) 


(Fab(x,y)F ab (x,y) - ^pab(x,y)p ab (x,yj^ , 
j]yF ab (x,y)p ab (x,y), (3.51) 


using the abbreviated notation J)* 2 dz = f** d z° J ^ d d z. Here we kept in 
(13.51(1 the general notation without summing over held indices, since the same 
expressions will be used for the case with a nonvanishing held expectation value 
below. We note that F(x, y) and p(x, y) along with the other statistical and 
spectral components of the self-energies are real functions. 


3.2.3 Nonvanishing field expectation value 

In the presence of a nonzero held expectation value (j) a the most general 
propagator can no longer be evaluated for the diagonal conhguration (13.71) . For 
the TV-component scalar held theory ((2.61(1 one has with <fi 2 = 0 a 0 a : 

M 2 ab {x- 0, F) = (m 2 + [F cc (x, x) + cp 2 (x)] ^ 5 ab 

+ [Fab{xpx) + 4) a {x)(t) b (x)} . (3.52) 

10 This follows from using the decomposition identity for the propagator (12.99[) and I(x,y ) = 
I F (x, y) - |/ p (x, y) sgn c (a; 0 - y°). 
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The self-energies Sf b (x y) = E^ b (x,y; (f), p, F) and Z p ab (x,y) = E p ab (x,y, (p, p, F) 
are obtained from (|3.44|) as 

£ ab ( x , 2/) = y ) bo(®)0ft(j/) + ^(x, 2/)] - ^ip(x, y)pab(x, y) 

+P F (x,y)F ab (x,y) - jP p (x,y)p ab (x,y) }, (3.53) 

s ab( x >2/) = “3^{ / p( x ’2/)[0a(a:)06(2/) + -P 1 a6(®,2/)] + ^F(a:,J/)Po6(®,2/) 

+P p (x,2/)Fa6(^,2/) + ^V(ic,2/)Paft(a;,2/)}- ( 3 . 54 ) 


The functions I F (x,y) = I F (x,y,p,F) and I p (x,y ) = I p (x,y,p,F ) satisfy 
the corresponding equations as for the case of a vanishing macroscopic field 
given above. The respective ^-dependent summation functions Pf(x, y) = 
Pp(x, ?/; 0, p, F) and P p (x, p) = P p (z, p; </>, p, P) are given by 


P F {x,y) = ~\H F (x,y)~ / dz [F p (x, z)/ F (z, y) + I p (x, z)H F (z, y)\ 


I to 


+ / dz [H f (x, z)I p (z , p) + /f(®, 2)-Hp(z, 2/)] 

J to 

nx° ny° 

- dz dv I p (x,z)Hp(z,v)I p (v,y) 

Jto Jto 

pz° 

+ dz dvI p (x,z)H p (z,v)I F (v,y) 

Jto Jto 

ry° ry° 

+ / dz d v Ip(x,z)H p (z,v)I p (v,y)>, 

Jto Jz° J 


(3.55) 


P P (x,y) = 


A 


3N 


Hp{x,y) - / dz [P p (x, z)I p (z,y) + I p (x, z)H p (z,y)j 


+ / dz dvI p (x,z)H p (z,v)I p (v,y)\, 


(3.56) 


with 


H F (x,y) = -4> a (x)F ab (x,y)4> b (y) , H p (x,y) = -<j> a (x)pab(x,y)</> b (y) ■ (3.57) 


The time evolution equation for the field (12. 11 7j ) for the 2PI effective action to 
NLO (I3.44p is given by 


Dl + W^ x) 


5 ab +Ml h {x-,(j) 



(j) b (x) 
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(3.58) 


_A_ 
3 N 



d y [Ip(x, y)F ab (x , y) + I F (x, y)pab(x , y)} <j) b (y) 



0 ,F,p) <j> b (y ) . 


3.2.4 Dephasing at leading order 


For simplicity we consider spatially homogeneous field expectation values 0 a (t) = 
(<F a (t, x)), such that we can use the modes F ab (t, t'\ p) and p a b(t, t'] p) after spatial 
Fourier transformation to describe the dynamics. 

The LO contribution to the 2PI effective action ([3739]) adds a time-dependent 
mass shift to the free field evolution equation. The resulting effective mass term, 
given by (13.52)1 for N —» oo, is the same for all Fourier modes and consequently 
each mode propagates “collisionlessly”. There are no further corrections since 
according to (13.40(1 and ((3.44(1 the self-energies T,p and T, p are 0(1/N) and vanish 
in this limit. Therefore, the evolution equations (I2.11ip for this approximation 
read: 


[dt + P 2 + M 2 (t] 0, F)] Fab(t, t’] p) = 0 , 
[dt + p 2 + M 2 (t] 0, F)] Pab (t, t'\ p) = 0, 


(3.59) 


9? + F^( () + m 2 («; 0 ,F) 


4> b (t) = o 


with 


M 2 (t; 0, F) = "i 2 +^ 


Fecit11] p) +0 2 (^) 


L-'p 


(3.60) 


where f p = f d d p/ (27r) d . In this case one observes that the evolution of F and 
<j) is decoupled from p. Similar to the free field theory limit, at LO the spectral 
function does not influence the time evolution of the statistical propagator. The 
reason is that in this approximation the spectrum consists only of “quasiparticle” 
modes of energy oj p (t) = y/p 2 + M 2 (t) with an infinite life-time. The associated 
mode particle numbers are conserved for each momentum separately. In contrast, 
in the interacting quantum field theory beyond LO, direct scattering processes 
occur such that mode occupancies can change. We note that there are also 
no memory integrals appearing on the RHS of (13.59(1 as a consequence of this 
approximation. 

As an example, we consider here the dynamics for d — 1 with the initial 
conditions at t 0 = 0: 

/p(0) + 1/2 


F(0,0; p) = 


dtF(t, 0; p) |t= 0 = 0 , 


.^p 2 + M 2 (0) 

F(0, 0; p)d t d t iF(t, t 1 ; p)| t= ^ =0 = [/ P (0) + 1/2] 2 


(3.61) 


0(0) = <9 t 0(f)| t=o = 0 


(3.62) 
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Figure 8: LEFT: Comparison of the LO and NLO time dependence of the 
equal-time correlation modes F(t,t] p) for the initial condition (13.611) . The 
importance of scattering included in the NLO approximation is apparent: the 
evolution deviates from the nonthermal LO solution and the initially peaked 
distribution decays, approaching thermal equilibrium at late times. RIGHT: 
Effective particle number distribution for a peaked initial distribution in the 
presence of a thermal background. The solid line shows the initial distribution 
which for low and for high momenta follows a Bose-Einstein distribution, i.e. 
ln[l + l// p (0)] ~ e p (0)/T 0 . At late times the nonthermal distribution equilibrates 
and approaches a straight line with inverse slope T eq > T 0 . 


and refer to the suggested literature in section l3~4l for further studies in dimensions 
d > 1 and for different initial conditions. Here F a b(t, p) = F(t, t'; p )S ab , which 
is valid for all times with these initial conditions. The mass term M 2 (0) is given 
by the gap equation (13.601) in the presence of the initial nonthermal particle 
number distribution (13.271) . 

As the renormalization condition we choose the initial renormalized mass in 
vacuum, tur = M(0)|/(o)=o = 1, as our dimensionful scale. In these units the 
particle number is peaked around |p| = p ts = 5m r with a width determined 
by a = 0.5 m R and amplitude A f = 10. We consider the effective coupling 

V(6 m? R ) = 1. 

On the left of Fig. |S]we present the time evolution of the equal-time correlation 
modes F(t, t\ p) for different momenta: zero momentum, a momentum close to 
the maximally populated momentum p ts and about twice p ts . One observes that 
the equal-time correlations at LO are strictly constant in time. This behavior can 
be understood from the fact that for the employed initial condition the evolution 
starts at a time-translation invariant nonthermal solution of the LO equations. 
There is an infinite number of LO solutions which are constant in time, depending 
on the chosen initial condition details. This is in sharp contrast to the well- 
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founded expectation that the late-time behavior may be well described by thermal 
equilibrium physics. Scattering effects included in the NLO approximation indeed 
drive the evolution towards thermal equilibrium, which is discussed in detail in 
section 13.2.51 below. 

A remaining question is what happens at LO if the time evolution does not 
start from a time-translation invariant nonthermal solution of the LO equations, 
as was the case for the initial condition employed above. In the left graph of 
Fig. [9] we plot the evolution of M 2 (t ) in the LO approximation as a function 
of time f, following a “quench” described by an instant drop in the effective 
mass term from 2 M 2 to M 2 at initial time. As a consequence, the mass term 
appearing in the initial conditions (I3.62p is not the same as the one appearing 
in the evolution equations (I3.59p . The initial particle number distribution is 
/ p (0) = l/(exp[-y/p 2 + 2M 2 (0)/T 0 ] — 1) with T 0 = 2M(0) and 0(0) = 0(0) = 0. 
The sudden change in the effective mass term drives the system out of equilibrium 
and one can study its relaxation. 

In Fig. [9] we present M 2 {t ) for three different couplings A = Ao = 0.5M 2 (0) 
(bottom), A = 10Ao (middle) and A = 40Ao (top). All quantities are given in 
units of appropriate powers of M(0). Therefore, all curves in the left graph 
of Fig. [9] start at one. The time-dependent mass squared M 2 {t) shoots up 
in response to the “quench” and stays below the value 2M(0) 2 of the initial 
distribution. The amplitude of initial oscillations is quickly reduced and the 
evolution is rapidly approaching constant asymptotic values. Since there are no 
direct scatterings included in this approximation, the damping is due to a simple 
dephasing phenomenon of harmonic oscillations with a time-dependent frequency 
shift. Correspondingly, the correlator F(t, f'; p), which is a function of the relative 
times t — t' and of t + t' at early times, loses its dependence on t + t' asymptotically 
such that F(t,t ; p) becomes a constant. 

To understand this in more detail, we compare the asymptotic values with the 
self-consistent solution of the LO mass equation (13.601) for constant mass squared 
M 2 p and given particle number distribution / p (0)l)j 


Mgap m ’ + 6 ./ 2tt l /p(0) + 2 J ^/p 2 + M 2 


(3.63) 


gap 


The result from this gap equation is M 2 ap = {1.01,1.10,1.29}M 2 (0) for the three 
values of A, respectively. For this wide range of couplings the values are in good 
numerical agreement with the corresponding dynamical large-time results, which 
can be inferred from Fig. [9] as {1.01,1.11,1.31}M 2 (0). One concludes that the 
asymptotic behavior at LO is well described in terms of the initial particle number 

11 Here the logarithmic divergence of the one-dimensional integral is absorbed into the bare 
mass parameter m 2 using the same renormalization condition as for the dynamical evolution 
in the LO approximation, i.e. m 2 + ^ f ^ (/ p ( 0) + ^) (p 2 + M 2 ( O))” 1 / 2 = M 2 ( 0). 
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Figure 9: LEFT: Shown is the time-dependent mass term M 2 (t ) in the LO 
approximation for three different couplings following a “quench”. All quantities 
are given in units of appropriate powers of M( 0). RIGHT: Time dependence 
of the equal-time zero-mode F(t, £; p = 0) after a “quench”. The inset shows 
the mass term M 2 (t), which includes a sum over all modes. The dotted lines 
represent the Hartree approximation (LO + ), while the solid lines give the NLO 
results. The coupling is X/6N = 0.17M 2 (0) for N = 4. 

distribution / p (0). We emphasize that the latter is not a thermal distribution for 
the late-time mass terms with values smaller than 2M 2 (0). 

The question of how strongly the LO late-time results deviate from thermal 
equilibrium depends of course on the details of the initial conditions. Typically, 
time- and/or momentum-averaged quantities are better determined by the LO 
approximation than quantities characterizing a specific momentum mode. This 
is related to the prethermalization of characteristic bulk quantities, and we refer 
to the further discussions of this topic to the literature listed in section 13.41 

Fig. [9] shows the equal-time zero-mode F(t, £; p = 0) along with M 2 (t), which 
includes the sum over all modes for comparison. Here we employ a “quench” 
with a larger drop in the effective mass term from 2.9M 2 (0) to M 2 (0). The 
initial particle number distribution is / p (0) = l/(exp[A/p 2 + M 2 (0)/T 0 ] — 1) 
with T 0 = 8.5M(0). In the figure the dotted curves show the dynamics obtained 
from an “improved” LO (Hartree) approximation, LO + , that takes into account 
the local part of the NLO self-energy contribution and is often employed in the 
literature. The resulting equations have the very same structure as the LO ones, 
however, with the LO and NLO contribution to the mass term M 2 (t ) included 
as given by (I3.47P below. The large-time limit of the mass term in the LO + 
approximation is determined by the LO + stationary solution in complete analogy 
to the above discussion. We also give in Fig. [9] the NLO results, which are 
discussed below. 

The effective loss of details about the initial conditions is a prerequisite 
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Figure 10: LEFT: Shown is the evolution of the unequal-time correlation 
F(t, 0; p = 0) after a “quench”. Unequal-time correlation functions approach 
zero in the NLO approximation and correlations with early times are effectively 
suppressed ( X/6N = (5/6 N ~ 0.083) M 2 (0) for N = 10). In contrast, there is 
no decay of correlations with earlier times for the LO approximation. RIGHT: 
The logarithmic plot of | p(t, 0; p = 0) | and | F(t, 0; p = 0)| as a function of time t 
shows an oscillation envelope which quickly approaches a straight line. At NLO 
the correlation modes therefore approach an exponentially damped behavior. (All 
in units of M(0).) 


for thermalization. At LO this is obstructed by an infinite number of 
spurious conserved quantities (mode particle numbers), which keep initial-time 
information. An important quantity in this context is the unequal-time two- 
point. function F(t, 0; p), which characterizes the correlations with the initial time. 
Clearly, if thermal equilibrium is approached then these correlations should be 
damped. On the left of Fig. [TU] the dotted line shows the unequal-time zero-mode 
F(t, 0; p = 0) following the same “quench” at LO as for Fig. [9] left. One observes 
no decay of correlations with earlier times for the LO approximation. Scattering 
effects entering at NLO are crucial for a sufficient effective loss of memory about 
the initial conditions, which is discussed next. 

3.2.5 Thermalization at NLO 

In contrast to the LO approximation, at NLO the self-energies and S p ~ 
0(1/IV) do not vanish. For the initial conditions (I3.6ip all correlators are diagonal 
in field index space and 0 = 0 for all times. In this case the evolution equations 
derived from the NLO 2PI effective action (13.401) are given by (12.1 lip with the 
self-energies (I3.47p - fl3.50p . As discussed above, the NLO evolution equations are 
causal equations with characteristic “memory” integrals, which integrate over the 
time history of the evolution taken to start at time to = 0. 
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We consider first the same initial condition (13.61 jl as for the LO case discussed 
above. The result is shown on the left of Fig. [ 8 ] for N = 10. One observes that 
the NLO corrections quickly lead to a decay of the initially high population of 
modes around p ts . On the other hand, low momentum modes get populated such 
that thermal equilibrium is approached at late times. 

lu order to make this apparent one can plot the results in a different way. For 
this we note that according to (I3.6ip the statistical propagator corresponds to 
the ratio of the following particle number distribution 

fp( t ) + \= [ F (t,t-,p)K(t,t-,p) -Q 2 (t,t-p)] 1/2 , (3.64) 


at initial time t = 0 and the corresponding mode energy. Here we have defined: 


K(t,t']p) = d t d? F(t,If ip) , Q{t,if\ p) = ^[d t F(t,t'-,p) + d t ,F(t,t']p)] . 

(3.65) 

Since we employed Q(0, 0; p) = 0 for the initial conditions (I3.6ip the mode energy 
is given by 


e P (f) 


f K(t,t; p) \ 1/2 

\F(t,t]p)J 


(3.66) 


at initial time, such that F(t, f; p) = [/ p (f) + l/2]/e p (f). For illustration of the 
results we may use (13.641) and (I3.66P for times t > 0 in order to define an effective 
mode particle number and energy , where we set Q(t,t'\ p) = 0 in (13.641) for the 
moment. 

The behavior of the effective particle number is illustrated on the right of 
Fig. [SJ where we plot ln(l + l// p (f)) as a function of e p (f). For the corresponding 
plot of Fig. Owe have employed an initial peaked distribution at Pt s /m R = 2.5, 
where we added an initial “thermal background” distribution with temperature 
T 0 /m R = 4@ Correspondingly, from the solid line ( t = 0) in the right figure 
one observes the initial “thermal background” as a straight line distorted by 
the nonthermal peak. The curves represent snapshots at equidistant time steps 
A tm R = 10. After rapid changes in f p (t ) at early times the subsequent curves 
converge to a straight line to high accuracy, with inverse slope T eq /To = 1.175. 
The initial high occupation number in a small momentum range decays quickly 
with time. More and more low momentum modes get populated and the particle 
distribution approaches a thermal shape. 

The crucial importance of the NLO corrections for the nonequilibrium 
dynamics can also be observed for other initial conditions. The right graph of 
Fig. [9] shows the equal-time zero-mode F(t,t ; p = 0), along with M 2 {t) including 
the sum over all modes, following a “quench” as described in section 13.2.41 While 
the dynamics for vanishing self-energies and is quickly dominated by the 


12 The initial mass term is M(0)/m R = 2.24 and A/6 N = 0.5m|j for TV = 4. 
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spurious LO stationary solutions, this is no longer the case once the NLO self¬ 
energy corrections are included. 

In particular, one observes a very efficient damping of oscillations at NLO. 
This becomes even more pronounced for unequal-time correlators as shown for 
F(t, 0 ; p = 0 ) in the left graph of Fig. [TQ1 On the right of Fig. [TO] the approach 
to an approximately exponential behavior is demonstrated for \p(t, 0 ; p = 0 )| and 
\F(t, 0; p = 0)| with the same parameters as for the left figure. The logarithmic 
plot shows that after a non-exponential period at early times the envelope of 
oscillations can be well approximated by a straight line. 

The strong qualitative difference between LO and NLO appears because an 
infinite number of spurious conserved quantities is removed once scattering is 
taken into account. It should be emphasized that the step going from LO to 
NLO is qualitatively very different than the one going from NLO to NNLO or 
further. In order to understand better what happens going from LO to NLO 
we consider again the effective particle number (13.641) . It is straightforward by 
taking the time derivative on both sides of (13.64(1 to obtain an evolution equation 
for f p (t) with the help of the relations (12.1 lift : 



(3.67) 

Here t 0 denotes the initial time which was set to zero in (12.1 lip without loss 
of generality. Since Ejr ~ 0(1/N) as well as £ p , one directly observes that at 
LO, i.e. for N —y oo, the particle number for each momentum mode is strictly 
conserved: d t f p (t) = 0 at LO. Stated differently, (13.64(1 just specifies the infinite 
number of additional constants of motion which appear at LO. In contrast, once 
corrections beyond LO are taken into account then (13.64(1 no longer represent 
conserved quantities. 

3.3 Transport equations 

We have seen in the previous sections that after a transient rapid evolution, 
which is characterized by damped oscillations of correlation functions with 
frequency determined by the renormalized mass for zero spatial momentum, 
the subsequent evolution of equal-time correlation functions such as F(t,t) p) 
becomes comparably smooth. Therefore, after the oscillations become effectively 
damped out, an approximate description taking into account only low orders 
in an expansion in derivatives with respect to the center coordinates (t + t') 


84 















is expected to be suitable. This will be used in the following to obtain a set 
of equations, which form the basis of kinetic theory and which will be used to 
describe the transport of conserved quantities in section | 6 ] 


3.3.1 Gradient expansion 

We first note that the spectral function (12.961) is directly related to the retarded 
propagator, Gr, or the advanced one, Ga , by 

G R (x,y) = @(x°-y°)p(x,y), G A (x,y) = -@(y°-x°)p(x,y). (3.68) 

Similarly, the retarded and advanced self-energies are 

Sr(x,?/) = 0(x° - ?/ 0 )£ p (x,;?/), T A (x,y) = -0(y°^ x°)T p (x,y ). (3.69) 

Here we consider the case of a vanishing held expectation value, </>(x) = 0. With 
the help of the above notation, interchanging x and y in the evolution equation 
(12. 11 ip for F(x, y) and subtraction one obtains 

(□* - + M 2 (x) - M 2 (y) ) F (.x , y) 

= j d d+1 z 9 (z°) (V (x, z ) £ a {z, y) + G R (x, z) T, F (z, y) 

- Sr (x, z) F (z, y) - (x, z) G a (z, y) j . (3.70) 

The same procedure yields for the spectral function 

(n x -n y + M 2 (x) — M 2 (y))p(x,y) 

= I d d+1 z (g r (x, z) S p (z, y) + p (x, z) £4 (.z , y) 

- Tip (x, z) G a (z, y) - Sr (x, z) p (z, y) ) . (3.71) 


So far, the equations (13.701) and (13.711) are fully equivalent to the exact equations 

mm- 

Transport equations are obtained by prescribing F, p and derivatives at a 
finite time using the equations with t 0 -)■ -(X) as an approximate description. 
Furthermore, one employs a gradient expansion to (j3.70[) and (13.711) . In practice, 
this expansion is carried out to low order in the number of derivatives with respect 
to the center coordinates 



(3.72) 


and powers of the relative coordinates 


s M = x M — y^ . 


(3.73) 
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Figure 11: For given finite 2X° = x° + y° the relative coordinate s° = — y° has a 

finite range from —2X° to 2X°. 


Even for finite X° one assumes that the relative-time coordinate s° ranges from 
—oo to oo in order to achieve a convenient description in Wigner space, i.e. in 
Fourier space with respect to the relative coordinates (13.73ft . This requires a loss 
of information about the details of the initial state, which enters in the derivation 
as an assumption, and limits the use of the approximate equations to not too early 
times. 

For the description in Wigner space we introduce the Fourier transforms with 
respect to the relative coordinates, such as 


F{X,u, p) 


r 2X° 


ds° e iuJS ° 


' —2X° 
r 2X° 


[ d d se~ ips F (x + Lx- 

J — OO ' " 

ds° cos(o;s 0 ) [ d d s e~ ips F (x+-,X-~), 

J-oo ' 2 2/ 


(3.74) 


using the symmetry property F(x, y ) = F(y, x ) for the second line. We emphasize 
that the time integral over s° is bounded by ±2X°. The time evolution equations 
are initialized at time x° = y° = 0 such that x° > 0 and y° > 0. According 
to (13.731) . the minimum value of the relative coordinate s° is then given by 
— y° = —2A"° for x° = 0 while its maximum value is x° = 2X° for y° = 0 as 
illustrated in Fig. [HI Similarly, we define 


p 2X° r oo 

p{X, u, p) = —i ds° e tU}S ° 


I- 2X° 
r 2X° 


= 2 


J A* ae .-»‘p(x + i,x-^ 

ds° sin(cus°) [ d d s e~ ips p (x + -, X - -) , (3.75) 

J — ^ 2 2 / 


where a factor of i is included in the dehnition to have p(X,oj, p) real and we 
have used p(x,y) = —p(y,x). The equivalent transformations are done to obtain 
the self-energies E^(X,o;,p) and E p (X,cu,p). 

In order to exploit the convenient properties of a Fourier transform it is a 
standard procedure to extend the limits of the integrals over the relative time 
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(3.76) 


coordinate in (13.74ft and (j3.75[) to ±oo. For instance, using the chain rule 

d d 


□a; - □„ = 2- 

y ds u dX» 


and the gradient expansion of the mass terms in (13.701) td 13 l NLO 


M [X + - 




(3.77) 


then the LHS of (13.70j) becomes in Wigner space 

d 8 . t dM 2 (X)' 


d d+i se i P ^ 


2_I_u s fl - 

ds^ dXv dX' 1 


F(X+ S -,X- S - 
1 2 2 




„ d 3M 2 (X) 3 

2ff -1- v ' 

P dX" 


F{X,p ) 


(3.78) 


dXv dp y 

with p° = co. 

Similarly, one transforms the RHS of (j3.70[) and (j3.71[) using the above 
prescription. To derive this expression, we consider for functions f(x,y ) and 
g(x,y ) the integral f d d+1 zf(x, z)g(z, y). For 


= f(X XI + S -f,X TZ - S -f 


(3.79) 


we introduce the Fourier transform with respect to the relative coordinate: 


/ (X xz ,p) = J d d+1 s xz exp {ippS^y} f (x xz + X xz - . (3.80) 

This quantity may also be written as 

f(X xz ,p) = f (x xy + ^-,p^) = e,J^^ jr \f(X X y,p), (3.81) 

where we Taylor expanded to obtain the last equality. Similar steps for g(z,y) 
lead to 


g (X zy , p) = j d d+1 s zy exp {ipys^} g (x zy + X zy - ^ 

= g (X xy ~ ^,p) = expj-^Q^j^^xmP)- (3-82) 
Inverting flihSOl) gives with (137811) 

f( x ,z) = f (2 Vi exp {-Wxz} exp 7r} f( X wP )» 

9 ( z ’ y ) = f ^y+i ex P {-W'M ex P {~ S -fd^x} 9{X xy ,p') (3.83) 

13 The expression is actually correct to NNLO since the first correction is O ((s M ch>) 3 ). 
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using also the equivalent expressions for g(z, y ). With this we can write employing 
a partial integration: 


d d+1 s xy exp{ip p s% y } J d d+1 zf{x : z) g(z,y) 
d d d d 


= exp - 


dp p dX^y dpdXHy 


f (■ X xy ,p)g(X' p') | x xy =x'p~ 


— f{X X y, p)g(X X y, p) + — {f(X X y, p)] g^Xxy, p)} PB , (3.84) 


where the Poisson bracket reads 


{f(X,p);g(X,p)} pB 


df(X,p) dg(X,p) _ df(X,p) dg(X,p) 
dpy dX' J dXv dpy 


(3.85) 


The lowest-order equations are obtained by neglecting O (dx^d P/1 ) and higher 
contributions in the gradient expansion. To this order the transport equations 
then read: 


dF{X,p) 
p dxv 


2p^ 


dp(X,p) 

dX^ 


t p (X, p) F (X, p) - Z F (X, p ) p (X, p) , (3.86) 

0. (3.87) 


We note that the compact form of the gradient expanded equation to lowest 
order is very similar to the exact equation for the statistical function (12.1 1 1)1 . 
The main technical difference is that there are no integrals over the time 
history. Furthermore, to this order in the expansion the evolution equation for 
the spectral function p(X,p) becomes trivial. Accordingly, this approximation 
describes changes in the occupation number while neglecting its impact on the 
spectrum and vice versa. Moreover, one observes that in thermal equilibrium the 
fluctuation-dissipation relations (12.157jl and (I2.158jl yield dxF(X,p) = 0. 


3.3.2 Vertex-resummed kinetic equation 

In the following we consider the lowest-order expressions (13.861) and (13.871) of 
the gradient expansion, where the self-energies are obtained from the 2PI 1/X 
expansion to NLO. It is convenient to introduce a suitable “occupation number 
distribution” f(X,p). For the real scalar field theory one has 

F(X,-p) = F(X,p) , p(X,-p) = -p(X,p). (3.88) 

Without loss of generality we can write 


F(X,p) 


f{x,p) + ^ J p(X,p ), 


(3.89) 


















which defines the function f(X,p ) for any given F(X,p) and p(X,p). We 
emphasize that without additional assumptions (13.89(1 does not represent a 
fluctuation-dissipation relation (12.1571) . which holds only if f(X,p ) is replaced 
by a thermal distribution function. We will not assume this in the following and 
keep f(X,p ) general at this stage. In particular, (I3.88j) then implies the identity 


f(X,-p) = ~(f(X,p) + 1) 


(3.90) 


For spatially homogeneous ensembles ((3.87)1 implies a constant p(p) that does 
not depend on time. In contrast, the statistical function F(t,p ) to this order can 
depend on time t = X° and using (13.89)1 we can write 

Here C[f] describes the effects of interactions to lowest order in the gradient 
expansion. 

It is instructive to consider for a moment a free spectral function given by 


p {0 \p) = 2n sgn(p°) <5((p 0 ) 2 — Wp) 


(3.92) 


for a relativistic scalar field theory with particle energy uj p . Choosing a free field 
theory type spectral function with zero “width” leads to the characterization of 
the dynamics in terms of a “gas” of particles, where 


f{t, P) = f{t,P° = w p ,p) = 


^-2 p°p {0) (p) f(t,p) (3.93) 


denotes the “on-shell” number distribution. In the following we will consider the 
resummed 1/IV expansion of the 2PI effective action to next-to-leading order. If 
the spectral function p(jp) for such an approximation is taken, it is not of the 
free field form (13.9211 but acquires a non-zero “width” as pointed out in section [3] 
We may nevertheless use the same procedure to define for spatially homogeneous 
systems an effective number distribution 


r°° dn° 

f(t,p ) = / —2 p°p(p)f(t,p), 


(3.94) 


which depends on time and spatial momenta only. This will allow us to formulate 
the quantum field theory in a way reminiscent of a Boltzmann equation. The time 
evolution of that distribution is, according to (13.86)1 . given by 


df(t, p) 
dt 


= C[f](t, p) = 


dp 0 

L 


%p(t,p)F{t,p)—X F (t,p)p(p) . (3.95) 




















For the following it is useful to note that the last expression can be rewritten 
with 


^ P {t,p)F(t,p) - E F (t,p)p{p) = 

F F {t,p) + ^F p (t,p) 

1 ~ , 

F(t,p) ~ ^p(p) 

— 

F F (t,p ) - - F p (t,p) 

F(t,p ) + ^p(p) 


(3.96) 

With (13.89(1 one also has 

F{t,p) -\p(p) = f{t,p)p(p), 

F{t,p) + ^p(p) = [f(t,p) + l]p(p), (3.97) 

which will allow us to conveniently express everything in terms of factors of f(t,p ) 
or, including “Bose enhancement”, [f{t,p) + 1] and p(p). 

In the following, we consider the collision term C[f] at next-to-leading order 
in the 1/N expansion of the 2PI effective action as described in section 13.21 
To this end, we introduce retarded and advanced quantities as I R (x,y ) = 
0 (x° - y°)I p {x,y), I A (x,y) = - x°)I p (x,y) and equivalently for fl R (x,y) 

and n A (x, ?/) using the expressions (13.50P and (13.51(1 with t 0 —oo. Then we 
Fourier transform with respect to the relative coordinates to obtain the lowest- 
order gradient expansion result. To ease the notation, we suppress the dependence 
on the global central coordinate and only write the momentum dependencies, with 
f q = f d d+1 q / {2n) d+1 . For the retarded and advanced summation functions one 
finds from (13.50(1 : 


In(t,p) 


n R(t,p) 

i + n R (t,p) 


lA(t,P ) 


n.4(t,p) 

i + n A (t,p)' 


(3.98) 


Proceeding in the same way also for Ip, one obtains with I p = Ir — I A - 


I pit, p) = v eS (t, p) II pit, p) , I F (t, p) = Ves(t, P) n F (t, p ). (3.99) 


The vertex correction v e g(t,p) sums the infinite chain of “ring” diagrams 
appearing at NLO in the 1/N expansion and is given by 


Veg(t,p) 


1 

l + Il R (t,p)\ 2 ’ 


(3.100) 


with |1 + n R | 2 = (1 + n R )(l + n A ) and the one-loop retarded self-energy 

n R {t,p) = ^ JF(t,p- q)G R (q) . (3.101) 
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At next-to-leading order in the 2PI 1/N expansion, the above linear combinations 
of self-energies are then given by 

1 y2 r 

± I v eS (t,p-q ) 


X 


’ ql 


F(t,p- q~l)± -p{p -q~l) 


F it,q) ± 2 M 




,(3.102) 


which has the structure of a two-loop self-energy, however, with a time and 
momentum dependent “effective coupling” entering via v e s(t,p). In terms of 
f(t,p ) the above self-energy combinations read 

1 y2r 

ZF(t,p) + -'Z p (t,p) = J Ves(t,p-q) 


x [f(t,p 

S F (t,p) - 
x f{t,p- 


- q - 0 +!] pip - q -1) [fit, q ) +1] piq) [fit, l) + 1] pit), (3.103) 
-q- 0 pip ~q~l) f{t, q)p(q) fit, l)p(l) . (3.104) 


Putting everything together one obtains 

(t p F - {t,p) = J (2Tr) d+1 5(p - q-l-r ) 

X Ves it, p~q) | [fit, q) + 1] [fit, l ) + 1] [fit, r ) + 1] f(t, p) 

~ fit. q)f{t, l)f(t, r ) [fit, p) + 1] } p{q)p{l)p{r)p{p) . (3.105) 

To bring this expression into a form which can be directly compared to kinetic 
or Boltzmann descriptions, we map onto positive frequencies. For this we split 
the frequency integrals dq° . . . = dq° ... + J 0 °° dg° ... with q° —y —q° for 
the negative frequency part and employ f(t, —q°, q) = — [f(t, q°, q) + 1]. 

Collecting the 2 3 = 8 contributions from the different orthants in frequency 
space one finds a “collision integral” (I3.95P that is accurate to next-to-leading 
order in the large -N expansion, including processes to all orders in the coupling 
constant. We obtain with f p = f(t,p ) suppressing the global time dependence: 


C' NL °[/](p) = f ^ 2 ^ 2 [f]ipJ,q,r) [ifp + l)ifi + l)f q f r -fpflifq + l)if r + l)] 

+ J dti(a) 3 [f](p, q. r ) [{fp + l)(/i + !)(/<? + l)/r — fpflfqifr + 1)] 

+ J d£l}ii) 3 [f}{p, l,q, r) [if p + Fjflfqfr — fpifl + yifq + l)(/r + 1 )] 

+ JdFL°^ 4 [f]{p,l,q,r) [ifp + l)ifl + l)(/g + l)(/r + 1) — fpflfqfr\ ■ 

(3.106) 
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Here 


d n 2 ^ 2 [f}(p,l,q,r) 




d ^lb) 3 [f]{p,l,q,r) 


d Q 0 ^[f](pJ,q,r) 


X 


/ lqr 


5(p + / — g — r) 


A 2 f 00 dp 0 d/ 0 dg°dr° 

18ZV J 0 (2vr) 4 "( rf + 1 ) 

PpPlPqPr [Ves(P + 0 + Uefffr - (?) + W e ff(P ~ f)] , 

A 2 Z ” 00 dp°dZ 0 dg 0 dr° 


18ZV J 0 (2vr) 4 -( d + 1 ) 


t lqr 


5(p + / + q - r) 


X PpPlPqPr [V e g(p + l) + V eS (p + q) + UeffO ~ r)] , 
A 2 dp°dl 0 dq°dr 0 " 


18ZV 7 o ( 2 vr) 4 -(^ 4 ) Jlqr 

X PpPlPqPr V e g(p ~ l) , 

A 2 r°° dp°dl 0 dq°dr 0 

18N J 0 NHWT Jlqr 

X Pr,PlPnPrV P ff(v + l) 


5(p — l — q — r) 


5(p + l + q + r) 


with J q = J d d g/(27r) d . We emphasize that the above expressions still contain 
the integrations over frequencies and spectral functions. No quasi-particle 
assumptions using a free-held form of the spectral function has been employed 
yet and the only approximations are the 1/iV expansion to NLO and the gradient 
expansion underlying (13.911) . 

One observes that for sufficiently large p, for which n ri (p) 1, the vertex 
function (13. lOOjl approaches one. In this case the 2 p 2 contribution of the first 
line in (13. 106j) is reminiscent of the two-to-two scattering process in a Boltzmann 
equation. The main difference is that the latter assumes a 5-like spectral function 
such that all momenta are on shell. Therefore, in the perturbative expression 
(I6.225P off-shell processes involving the decay of one into three particles or 
corresponding 3 —> 1 annihilation processes or even 0tt4 processes are absent. 
They can occur in principle at NLO in the 2PI 1/iV expansion, which leads to 
the different terms contributing to the RHS of (13. 106j) . but they are typically 
small. At sufficiently high momenta these off-shell contributions should be 
suppressed along with all quantum-statistical corrections such that the spectral 
function approaches a 5-like behavior. In this case one recovers the standard 
Boltzmann equation for elastic two-to-two scattering. In the infrared v e g (p) may 
have a nontrivial momentum dependence, which incorporates important vertex 
corrections for the 2 H 2 scattering term that are discussed, in particular, in 
section 1 
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4 Classical aspects of nonequilibrium quantum 
fields 

It. is an important question to what extent nonequilibrium quantum field theory 
can be approximated by classical-statistical field theory. It is a frequently 
employed strategy in the literature to consider nonequilibrium classical dynamics 
instead of quantum dynamics since the former can be simulated numerically 
up to controlled statistical errors. Classical-statistical field theory indeed 
gives important insights when the number of bosonic field quanta per mode 
is sufficiently large such that quantum fluctuations are suppressed compared 
to statistical fluctuations. We will derive below a sufficient condition for the 
validity of classical approximations to nonequilibrium quantum dynamics. The 
description in terms of spectral and statistical correlation functions as introduced 
in section 12 .6.11 is particularly suitable for comparisons since these correlation 
functions possess well-defined classical counterparts. 

Classical Rayleigh-Jeans divergences and the lack of genuine quantum effects 
- such as the approach to quantum thermal equilibrium characterized by Bose- 
Einstein statistics — limit the use of classical-statistical field theory. To find out 
its use and its limitations we perform below direct comparisons of nonequilibrium 
classical and quantum evolutions for same initial conditions. One finds that 
classical methods can give an accurate description of quantum dynamics for the 
case of large enough characteristic occupation numbers. This typically limits 
their application to not too late times, before the approach to quantum thermal 
equilibrium sets in. 

Classical methods have been extensively used in the past to rule out 
“candidates” for approximation schemes applied to nonequilibrium quantum field 
theory. Approximations that fail to describe classical nonequilibrium dynamics 
should be typically discarded also for the quantum case. If the dynamics is 
formulated in terms of correlation functions then approximation schemes for 
the quantum evolution can be straightforwardly implemented as well for the 
respective classical theory. For instance, the 2PI 1/A-expansion introduced in 
section 13.21 can be equally well implemented in the classical as in the quantum 
theory. Therefore, in the classical-statistical approach one can compare NLO 
results with results including all orders in 1/N. This gives a rigorous answer to 
the question of what happens at NNLO or beyond in this case. In particular, for 
increasing occupation numbers per mode the classical and the quantum evolution 
can be shown to approach each other if the same initial conditions are applied 
and for not too late times. For sufficiently high particle number densities one 
can therefore strictly verify how rapidly the 1/N series converges for far-from- 
equilibrium dynamics. 

Though we restrict ourselves here to scalar fields, it is important to note 
that crucial applications arise also by coupling classical-statistical bosonic fields, 
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including lattice gauge fields, to fermions. For instance, this is done to 
describe the nonequilibrium real-time evolution in quantum electrodynamics 
for ultrastrong laser fields or quantum chromodynamics related to collision 
experiments of heavy nuclei. Fermions are never classical in the sense that they 
cannot be strongly occupied because of the Pauli exclusion principle. However, 
since the fermions occur quadratically in the respective actions, their functional 
integral can be treated on the lattice without loosing their genuine quantum 
nature. We refer to the reference list in section 14.51 for further reading on these 
extensions of the classical-statistical lattice simulation approach. 

4.1 Nonequilibrium quantum field theory revisited 

4.1.1 Functional integral 

In order to discuss the different origins of quantum and of classical-statistical 
fluctuations, it is convenient to rewrite the nonequilibrium generating functional 
of section 12.31 We start by employing the notation introduced in section 12.21 
where the superscripts ‘+’ and ’ indicate that the fields are taken on the 
forward branch (C + ) starting at to and backward ( C ~) along the closed time 
path, respectively. To be specific, we consider again the TV-component scalar 
field theory for which the classical action (12.61(1 can be written as 

S[<p + ,ip-} 

= j t - y^(^a W - ^ W) 2 

- \d^a(x)d^-(x) + ?L(p-(x)<p-(x) + ^ (y-(x)y-(a;)) 2 J . (4.108) 

Here the minus sign in front of the ’ terms accounts for the reversed 
time integration of the closed time contour and f x tg = f /g dx° f d d x. The 
corresponding generating functional for correlation functions then reads in this 
notation 


Z[J+,J-,R++R+-,R-+,R-;e 0 ] = J [d<pt][d<Po] eo [<pt, <Po] 



(4.109) 

Again the superscripts ‘+’ and ’ indicate that the respective time arguments 
of the sources are taken on the forward or backward branch of the closed time 
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path, and the prime on the measure means that the integration over the fields at 
initial time t 0 is excluded. 

In order to simplify the comparison with the classical-statistical field theory, 
a standard linear transformation A of the fields is introduced as 


where 



(4.110) 


(4.111) 


such that ip = (ip + + ip~)/2 and <p = ip + — ip ~, or if ± = tp ± <^/2, respectively. 
To avoid a proliferation of symbols we have used here ip , which agrees with 
the defining field in (12.61 j) only for ip + = ip~. Since this will be the case for 
expectation values in the absence of sources, where physical observables are 
obtained, and since there is no danger of confusion in the following we keep 
this notation. 

Correspondingly, we write for the source terms 



A 



1 


R F 

R a 



A 


R ++ R + - 
R~ + R— 


A T . 


(4.112) 

(4.113) 


Inserting these definitions into the functional integral (14.1091) and using that 
(14.112p and (14.1131) can be equivalently written as 



( R F R a \ 

V r k r f J 



(4.114) 

(4.115) 


one finds: 

Z[J, J, R f , R r : R a : R f ; eo] = [ [d^ 0 ] [d<£ 0 ] Qo Wo + <A)/ 2 > ipo - <Po/2} 


x 


yip^'ip exp i\s[(p,0\ + [ (<p a (x),<pa(x)) ^ 


^ 0,^0 

A I 


' X,to 


xy,t 0 


(IK igSi) (S3 


(4.116) 


Here S[ip, ip] = So[<^, (p\ + <Si nt [</?, (p\ consists of the action for the free field theory 


So[(p,0\ 



[d^ip a (x)d IJ ,ip a (x) - m 2 <p a (x)ip a (x)] 


(4.117) 
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Figure 12: Classical (left) and quantum vertex (right) in the scalar held theory. 


and the interaction part 


'hint \_Fi F] 


A 


<, Pa(x)Pa{x)<P b (x)(p b (x )- 


A 


6 NJ Xtto r ^ ,r “' ' 24A" J xto 


(p a {x)(p a (x)(p b (x)^ b (x). 


(4.118) 

The two types of vertices appearing in the interaction part (j4. 1 18f) are illustrated 
in Fig. [T2J To understand their role for the dynamics, it is important to 
note that one can also write down a functional integral for the corresponding 
nonequilibrium classical-statistical held theory. This is done in section l4~2l where 
we will observe that the generating functionals for correlation functions are very 
similar in the quantum and the classical-statistical theory. A crucial difference is 
that the quantum theory is characterized by an additional vertex: The interaction 
term ~ tp 3 , appearing with (14. 1 18j) in the functional integral of the quantum 
theory, does not occur for the classical theory. 


4.1.2 Connected one- and two-point functions 

For later use, we define in the following correlation functions in the above basis. 
From the generating functional for connected correlation functions, W — — i In Z, 
we define the macroscopic held (ft a and (j) a by 


SW 

SJ a {x) 


= M x ), 


SW 

5J a (x) 


= <j>a{x) ■ 


(4.119) 


The connected statistical correlation function F ab (x,y), the retarded/advanced 
propagators G^ A (x,y), and the “anomalous” propagator F ab (x,y) are defined 
by 


SW 


SR ab( X ^) 

sw 

t R ab( X ,y) 

sw 

SR ab( X ,V) 

sw 

SR ab( X ,y) 


(<l> a (x)<i> b (y) + F ab (x,y)) 

M x )My) - iG ab( x ,y ) 

7i [4>a(x)(j) b (y) - iG%,(x,y) 
4>a(x)4> b (y) + F ab (x, y) 


(4.120) 
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Equivalently, we can define the same connected two-point correlation functions 
by the second functional derivatives 


8 2 W 

8J a (x)8J b (y ) 
8 2 W 

8J a (x)8J b (y) 


Gab(x,y) 

iF ab (x,y) 


S 2 W 

8J a (x)8J b (y) 

8 2 W 

5J a (x)5J b (y) 


G^,(x,y) , 
iF ab (x,y) . 


(4.121) 

(4.122) 


We note that the propagators satisfy the symmetry properties G A b (x,y) = 
Gf a (y,x), F ab (x,y ) = F ba (y,x ) and F ab (x,y) = F ba (y,x). These properties 
follow directly from the definition of the propagators in terms of the second 
functional derivatives with respect to J and J. The spectral function p is 
given by the difference of the retarded and advanced propagators p a b{x,y) = 

Gab( x > y) - G tb( x , y), and G ab( x ’ y) = Pab(x , y)9(x° -y°). 

It is important to note that the anomalous propagator F vanishes in the limit 
where the external sources are set to zero. This is a consequence of the algebraic 
identity (12.13[) . since F = G ++ + G — G~ + — G + ~ , as one may readily check 
by changing the basis: 


= A 


G ++ 

G~ + 


G + ~ 

G- 


A 7 


F -iG R 

iG A F 

[G++ + G + G+~ + G~+} /4 [G++ - G + 
[G++ - G - G~+ + G+~] /2 G ++ + G— - G~ 


- G+~] /2 
-G + - 


,(4.123) 


More generally, in the absence of sources, we have 


8W 


8Ja(x ) 


8 2 W 


8J a (x) 8J b (y ) 


J , J , R R ’ A ’ F ’ P '=0 
J , J , R R ’ A ’ F ’ P =0 


4>a(x) = 0 , 

iF ab (x,y ) = 0 , 


(4.124) 

(4.125) 


and, correspondingly, arbitrary functional derivatives of the generating functional 
with respect to J vanish in the absence of sources. 

We finally mention that the inverse of the two-point function matrix (I4.123P 
in the absence of sources, where F — 0, reads 


0 

i(G R ) _1 


i(G A )- x 

(GR)" 1 • F ■ (G A )- X 


(4.126) 


Here we employed the compact matrix notation 


[(G r ) 1 -F-(G a ) l ] ab (x,y) = f (G R )J(x,z)F cd (z,w)(G A ) d 7 (w,y). 

J zw 

(4.127) 
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For vanishing sources, the exact inverse two-point function (14.126j> can then be 
written as 


( 0 i(G A )~ 1 \ = / 0 Go 1 \ _ / 0 -*S A \ 

0 i(G R ) _1 (G R ) _1 • F ■ (G A )~ l J { Gq 1 0 ) \ -i£ R E F )' 

(4.128) 


where the retarded, advanced and statistical self-energies are obtained as 

0 -*£ a \ , 1]T / £+ + -E+- 

-iZ K E F J v A M -E-+ E— 


A 


-i 


(4.129) 


from the self-energies written in the ‘±’ basis. To the retarded/advanced self¬ 
energy E r / a and the statistical self-energy E^ contribute only graphs with 
propagator lines associated to G R - A and F, which can be obtained from closed 
two-particle irreducible graphs by opening one propagator line. 


4.1.3 2PI effective action 

The 2PI effective action written in terms of the rotated variables becomes a 
functional of the field expectation values q i>, 0, and the propagators G R , G A , F, F, 


T = W - J (^> a {x)J a (x) + 4> a (x) J a (x)^j 

- \ J {ftab( X i y) + F ab(x , J/)) 

+ F ab( x ,y) ( M x )My) - iGa b ( x ,y)) 

+ Rab( x ,y) (<i>a{x)<l>b{y) - iG A b {x,y)^j 

+ R ab( X ,y ) ($a(x)$b(y) + Fab( x ,y) S ) Y (4.130) 

From (I4.130P one observes the equations of motion for the fields 


6T 

5cp a (x) 

6T 

64> a {x) 


-Ja(x)- / ( R^b{x,y)<j> b {y ) 

Jy v 

+ ^Rabi x >y)My) + ^My)R^ a (y^ x )), 

-J a (x)~ j (R%b(x,y)<j> h (y) 

Jy v 

+^Rab( x ,y)My) + \My)Rt(y,x)y 


(4.131) 


(4.132) 


as well as for the two-point functions 
5T Ip.. 

l 6G*(x, y) = ~2 R ^ X ' V) ’ 


. ST 
l 6G% > (x,y) 

ST 

SF ab (x,y ) 


~^Rab(x,y), 

R ab(x,y )• (4.133) 
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In the presence of a non-vanishing held value, 0^0 but 0 = 0, the interaction 
vertices are obtained from (j4. 1 18j) by shifting in S[(p, (p\ the held (p —$■ 0 + 93 , and 
collecting all cubic and quartic terms in the fluctuating helds <p and 0 , i.e. 




(4.134) 


The quadratic terms in the fluctuating helds are taken into account in the classical 
inverse propagator corresponding to a held dependent iG^ x (x, y\ 0) given by 
(I2.70p . For the classical-statistical held theory only those terms of (14.13411 appear, 
which are linear in (p. This is shown in the following. 

4.2 Functional integral for the classical-statistical theory 

The classical held equation of motion for the ^-component scalar held ip a (x) with 
action (12.6111 is given by 



(4.135) 


Its solution, ip C a( x ), requires the specification of the initial conditions <^(£ 0 , x ) = 
ipo, a ( x ) and 7r^(t 0 ,x) = 7T 0ia (x), with ir£(x) = d x oip c ^{x) for the considered scalar 
held theory. We define the macroscopic or average classical field by 



(4.136) 


Here W cl [<po,no\ denotes the normalized probability functional at initial time. 
The measure indicates integration over classical phase-space, 


AT 



(4.137) 


a= 1 x 


and the theory may be defined on a spatial lattice. Similarly, the connected 
classical-statistical propagator Fal(x,y) is defined by 
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The classical equivalent of the quantum spectral function is obtained by replacing 
—i times the commutator with the Poisson bracket 0 

Pab( X ,y) = -({<Pa(x),<Pb(y)} Pb)c 1 - (4.140) 

As a consequence, one finds the equal-time relations for the classical spectral 
function 


Pab( x iy)\x°=y° = 0, d x op* b (x,y)\ x o =y o = SabS(x-y). (4.141) 

Though their origin is different, we note that they are in complete correspondence 
with the respective quantum relations (12.9811 . 

In general, classical-statistical correlation functions are obtained as phase- 
space averages over trajectories given by solutions of the classical field equation 
(I4.135p . Such averages, for an arbitrary functional of the field f[(p a ], are defined 
as 

(f[<Pa])d = j [d^o][dvr 0 ] kh cl [<A), vr 0 ] f[<pf]. (4.142) 

This expression will be the starting point for constructing a functional integral 
for the classical-statistical field theory similar to the expression (I4.116P for the 
quantum theory. 

We define 

S cl [<p, 0\ = S 0 [ip, (p\ + Sfrjp, 01 • (4.143) 

The free part, S 0 [(p,f>\, is given by (I4.117p . Integrating by parts, we have 
to take care of the finite initial-time boundary values <j5 0 ,a( x ) = x) and 
7r 0ja (x) = d x oip(x)\ x o =to and we can write 


So[<p,0\ = ~ 7r 0 ,a(x)<p 0 ,a(x) + / f> a {x) (~D X - m, 2 ) . (4.144) 

J X J x,to 

The interaction part, S'djy?, <)5], reads 


5 intb,$ = 


A 


6N 


<Pa(x)(p a {x)cp b (x)(p b (x) . 


(4.145) 


This interaction part differs from (14.118ft in that it contains fewer vertices. The 
absence of vertices beyond those which are linear in <p turns out to be a crucial 
difference between a classical-statistical and a quantum field theory as is shown 

14 Recall that the Poisson bracket with respect to the initial fields is 


{A(x),B(y)} PB 



5A{x) SB(y) 

6ip 0 ,a(z) <57r 0i a(z) 


8A{x) SB(y) 
Sn o,a(z) 8<po,a(z) 


(4.139) 
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in the following. From the definition of S^ 1 ^, </5], the classical equation of motion 
(I4.135P for the field ip a (x) can be obtained as 


^^{x^ = - m 2 - -jTfi<P b (x)<Pb(x)^ <Pa(x) = 0. (4.146) 

We now rewrite this equation of motion as a ^-constraint in a functional integral 
using the Fourier transform representation 


5 


'SS d [(p,0\ 


6<p 

= / 3>{p exp <( i 


'x,t 0 
■jcl 


<fa(x) ( -Cb - m 2 - ^if b (x)(f b (x) ) lp a (x ) 


= / ®<f> exp iS cl [(fi,ip] +i / 7T 0ja (x)^0,a(x) > , 


(4.147) 


where we have used (14.144)1 for the last equality. 

In order to complete the construction of the functional integral for classical- 
statistical correlation functions, we note that any functional f[(p cl ] of the classical 
field solution can be written as 


/b cl ] = / [p~v d ] 

’<5S c1 [<a <p] 


fo 


^V/b] 8 


5ip 


j>] 


^0 


^0 


@ , (p@<pf[(p] exp \iS cl [ip,{p\ +i / vr 0ia (x)^ 0iQ (x) \ J7[^].(4.148) 


The prime on the measure again signifies that no integration over the initial 
time is implied, since this is fixed by the initial condition. What we have 
employed in the second equality of (I4.148P is the generalization of S (g(x)) = 
'^2 i 5(x — Xi)/\g'(xi)\ for zeros Xi of an ordinary function g(x) to functionals. The 
Jacobian reads 




det 


V 8<p8<p ) 


(4.149) 


Here it turns out that the Jacobian plays the role of an irrelevant normalization 
constant. Putting things together, we can write 


(f[<fia\) ci = I [dip 0 \[d7i 0 } lF cl [^ 0 ,7r 0 ] 


(4.150) 


x 


J @'(p @(p f[<p] exp jiS^b,^] +i J vr 0ia (x) 9 ? 0 ,a(x)| J[<p\. 
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The initial conditions may also be specified by a classical density matrix 
Qq bo + <,5o/2, v^o — <Po/ 2], which is characterized by the Fourier transform of the 
phase-space probability distribution VF cl boj tto] ; 


Qo bo + <A)/2,<A) - <A)/ 2 ] 


J [d7r 0 ]fF cl [^o,7r 0 ]exp 


F 0 ,a 


(x) ^0,a 



(4.151) 

Adding also sources, we obtain the generating functional for classical-statistical 
correlation functions: 


Z cl [J, J, R*', R k , R a , R*'- qoI = I [d^ 0 ][d^o] 6o bo + <Po/2, ip 0 - <p 0 /2] 


r?R r? A rW 


X 


fo^o 

-1 / 


b exp 2 <j A b) + / (^ab)bab)) 

J X, to 


Ja(x) 

Ja{pc) 


xy,t 0 




Comparing with the quantum generating functional in (14. 1 16j) . and using that 
the Jacobian J\^p\ plays the role of an irrelevant normalization constant, we find 
that the generating functionals for correlation functions are very similar in the 
quantum and the classical-statistical theory. The main difference is that the 
quantum theory is characterized by more vertices. 

In particular, all definitions for correlation functions given in section 14.1.21 for 
the quantum theory apply as well for the respective classical correlators. For 
instance, from (14.1221 ) we infer for the classical spectral function in the absence 
of sources: 


Prifay) = i J [d^ 0 ][d^o] Qo bo + £o/2, <£o - <fo/ 2] 

x J ®'ip S>'(p [(fi a (x)(p b (y) - <p a (x)ip b (y)] . (4.153) 

^0,00 

In order to understand the equivalence with the definition as a phase-space 
average of the Poisson bracket in (14.140j) . we may consider the retarded case 
x° < y° first and take x () = t 0 . With the help of (14.151j) we can then write 


-i J [d(^ 0 ] [cbo] Qo bo + V^o/2, ipo - <p 0 /2] J9'<p &<p <p a (x)<p h (y) e lScl[ifiM 

<Po,<fo 

= -i y[d(^o][d7ro]W cl bo,vro] J&<p ®(p(p a {x)ip b {y) e iSd [^ ] +>^»^°,»Jb] 

^0,^0 

- /[d^o][d7r 0 ]IF cl bo,7ro] x \ . [&<p®<p<p b {y) e iS ° l ^ +i J[tp\ 
J Jvr 0 , a (x) J 
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= - j [dipo][d7lo]W[ip 0 ,7T 0 ] 
= - / [d^o][d7r 0 ]W^o,7r 0 ] 


8<Pb(y) 

'0,o (x) 
cl 


<57T| 


(tp.x) (y) 

<fy?o,c(z) <57 To, c (z) 


where for the last line we have employed (j4.148j) and summation over repeated 
held indices is implied. Doing the corresponding steps for the retarded case 
x° > y° taking y° = t 0 shows the equivalence. 


4.3 Classicality condition 


Since the generating functionals for correlation functions are very similar in the 
classical and the quantum theory, the same techniques can be used to derive 
time evolution equations of classical correlation functions. In particular, the 
classical dynamic equations have the very same form (12.1111) and (I2.117|) as their 
quantum analogues with the replacements 0 a (x) —> 0^, F ab (x,y ) —> F£l(x,y), 
and p a b{x,y ) —$■ pf b (x,y). The corresponding classical-statistical self-energies 
S^ cl (x, y) and D^ cl (x, y) have in general the same diagrammatic contributions 
but are lacking certain terms due to the reduced number of vertices. 

In order to compare classical and quantum corrections to self-energies, we note 
that the classical-statistical generating functional (14.152ft exhibits an important 
reparametrization property: If the fluctuating fields are rescaled according to 

ip a {x) ->• (p' a (x) = V\ip a {x) , Cp a (x) -4 <p' a (x) = ~^= 0 a (x) ( 4 . 154 ) 


then the coupling A drops out of S' cl [(^, (p\ = 5' 0 [^, (p\ + S'fnJvb 0\ defined in (14.1171) 
and (14.1451) . The free part S'ofcp, c/3] remains unchanged and the interaction part 
becomes 

J <Pa{ x ) ( Pa{x)(p , b (x)lf , b (x). (4.155) 

Moreover, the functional measure in (14.1521) is invariant under the rescaling 
(14.1541) , and the sources can be redefined accordingly. Therefore, the classical 
statistical generating functional becomes independent of A, except for the 
coupling dependence entering the probability distribution fixing the initial 
conditions. Accordingly, the coupling does not enter the classical dynamic 
equations for correlation functions. All the A-dependence enters the initial 
conditions which are required to solve the dynamic equations. 

In contrast to the classical case, this reparametrization property is absent for 
the quantum theory: After the rescaling (14.154)1 one is left with S[(p', ip'] whose 
coupling dependence is given by the interaction part 




6N 


<Pa(x)<Pa(x)<p' b (x)(p'b(x)- 


A 2 


2dN 


<P'a(x) < P'a( X ) < P' b (xW b (x), 

(4.156) 
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according to (14.1181) . Comparing to (j4. 1 55j) one observes that the quantum vertex, 
which is absent in the classical-statistical theory, encodes all the A-dependence of 
the dynamics. 

The comparison of quantum versus classical dynamics becomes particularly 
transparent using the above rescaling. The rescaled macroscopic held and 
statistical correlation function are given by 

= V\M X ) , F ' ab (x, y) = ^F ab (x, y ), ( 4 . 157 ) 


while the spectral function p ab (x,y) remains unchanged according to (14.153[) . 
Similarly, we define for the statistical self-energy E^' (x, y) = AS ^ b (x,y). For 
instance, for the two-loop self-energies for vanishing macroscopic held given in 
(14.1581) and (14.1591) we obtain 


£' f (x,2/) 

E p (x,y) 


N + 2 
18N 2 
N + 2 
61V 2 


F\x,y ) 

p{x,y) 


F' 2 (x,y)- 7 


F“(x,y)- 


12 


A V(x,y) 
A V(x,y) 


(4.158) 

(4.159) 


The corresponding classical-statistical self-energies E f p(x,y) and E f(x,y) are 
given by the same expressions by dropping the A-dependent terms in (14.158f) 
and (I4.159p . which are proportional to p 2 . More precisely, one observes that the 
quantum evolution equations would be accurately described by the classical ones 
if the classicality condition 


F 2 (x,y ) » p 2 (x,y) 


(4.160) 


in terms of the non-rescaled correlation functions holds. A similar analysis of 
the loop-corrections in the presence of a nonzero macroscopic held given in 
section 13.1.21 or for the NLO 1/N expansion of section 13.2.31 yields the same 
condition. 

However, the requirement to fulfill (14.160jl for all space-time arguments is too 
restrictive and it can be typically only achieved for a limited range of time and 
relevant momenta. One expects that the classical description becomes a reliable 
approximation for the quantum theory if the number of held quanta in each mode 
is sufficiently high. The classicality condition (14.160jl entails the justification of 
this expectation. In order to illustrate the condition in terms of a more intuitive 
picture of occupation numbers, we employ the free-held theory type form of the 
spectral function and statistical propagator with mode frequency oj p as given 
in (13.331) . From this one obtains the following estimates for the time-averaged 
correlators: 


F 2 (M'; P) = ^ / df'F 2 (M'; P) 
J t—2ir/ujp 




(/ p (t) + 1/2) 2 


2uNt) 


p 2 (M';p) 


-+ 
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Inserting these estimates in (14.160j> for equal momenta yields 


fp{t) + \ > 1 or 


m »\ ■ 


(4.162) 


If the dominant momentum modes have occupancies much larger than the 
quantum-half, then a classical-statistical description can be an accurate 
approximation of the quantum dynamics. It is important to note that this is in 
general not the case in thermal equilibrium, where at temperature T the typical 
momenta p ~ T have an occupancy of order one. Consequently, the late-time 
approach to thermal equilibrium in quantum theories is beyond the range of 
applicability of classical-statistical approximations. Important examples, where 
a classical-statistical description is accurate, include continuous thermal phase 
transitions, since the relevant momenta for scaling behavior have p T with 
occupancy T/p> 1. Other examples are nonequilibrium instabilities, which yield 
high occupation numbers of characteristic modes, such as described in section O 
or wave turbulence and nonthermal scaling phenomena such as described in 
section 1 

4.4 Precision tests of quantum versus classical-statistical 
dynamics 

When the nonequilibrium quantum dynamics of a highly occupied system can 
be accurately mapped onto a classical-statistical held theory evolution, one can 
solve it without further approximations using lattice simulation techniques. This 
mapping is valid as long as the classicality condition (14.160(1 is fulfilled for typical 
momenta. The lattice held theory is then dehned on a spatial grid with spacing 
a s and side length N s a s in a box with periodic boundary conditions, as explained 
around equation (13.32jl at the end of section lB.1.31 Classical-statistical simulations 
consist of numerically solving the classical held equations of motion and Monte 
Carlo sampling of initial conditions according to (I4.142p . Thus, observables are 
obtained by averaging over the different classical trajectories that arise from the 
different initial held configurations. 

For the relativistic second-order differential equation of the scalar held (j4.135j) . 
the numerical integrations on a d -dimensional grid can be efficiently done using 
a standard leapfrog algorithm. For the nonrelativistic hrst-order equation (11.301) 
a conventional split-step method may be applied. These classical equations are 
then supplemented by suitable quantum initial conditions. For instance, typical 
Gaussian initial conditions represented by spatially homogeneous ensembles for 
a relativistic held equation as explained in section 12.41 can be specihed by a 
macroscopic held 0 (£q) = 0 o, its derivative 0(fo) = 0o and the distribution 
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function / p (i 0 ) as 


<p(to, x) = </>0 + 


d<i p / /p(^o) + l/2 ipx 

(27r) d y o/ p (* 0 ) P 


(4.163) 


with the initial frequency u; p (io) = \/p 2 + M 2 . For Gaussian initial conditions, 
the coefficients c p have to satisfy the relations 

(cpCp,) d = (27r) d <5(p - p') , (cpCp') c i = (CpCp,) d = 0 . (4.164) 

They can be realized by taking c p as Gaussian random numbers multiplied by 
complex random phase factors. This can be formulated as c p = 7L(p) e l27Ta (p) with 
a Gaussian distributed amplitude H(p) and uniformly distributed phase a(p) 
between 0 and 1. The random numbers c p have additionally to satisfy ck p = c p 
to ensure that ip(to,x) is real-valued. 

The conjugate momentum field tt (i 0 ,x) = <£>(t 0 ,x) is initialized in a very 
similar way: 

7 r(t o ,x)= 0 o + J AP_^f p ( to ) + l/ 2 )u p (t 0 )c p e^\ (4.165) 

with complex Gaussian random numbers c p satisfying the same relations as c p 
such as (14. 1641) . Since c p and c p are independent random numbers, one also has 
(c p c p /)ci = 0. Therefore, (cp(t 0 , x)7r(t 0 , y) + 7r(t 0 , x)ip(to, y)) c i vanishes automati¬ 
cally at initial time t = t 0 for the initial conditions discussed here. For instance, 
the macroscopic field or the statistical correlation function at times larger than t 0 
are then given by the averages (I4.136P or (j4.138j) . and equivalently for higher sta¬ 
tistical n-point correlation functions. Since p cl (t, t', x — x') is given by the Poisson 
bracket (14.1401) . spectral properties cannot be obtained from simple products of 
classical fields. However, in contrast to solving the equations (12.1111) for corre¬ 
lation functions, knowledge of the spectral function is not a prerequisite for the 
computation of correlation functions from classical-statistical simulations. 


In the following, we apply these simulation techniques to the scalar N- 
component field theory. Since the nonequilibrium evolution of classical-statistical 
correlation functions can be obtained numerically up to controlled statistical 
errors, the results include all orders in 1/N. Consequently, they can be used 
for a precision test of approximation schemes such as the 2PI 1/N expansion 
implemented in classical-statistical field theory. We emphasize that this compares 
two very different calculational procedures: the results from the simulation 
involve thousands of individual runs from which the correlators are constructed, 
while the corresponding results employing the 2PI 1/N expansion involve only 
a single run solving directly the evolution equation for the correlators. The 
accuracy of the simulations manifests itself also in the fact that the time-reversal 
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Figure 13: Left: Unequal-time two-point function F(t,0;p = 0) at zero 
momentum for N = 2,10, 20. The solid lines show results from the NLO classical 
evolution and the dashed lines from the full classical-statistical simulation. One 
observes a convergence of classical NLO and full results already for moderate 
values of N. Right: Damping rate extracted from F(t,0;p = 0) as a function 
of 1/N. Open symbols represent NLO and full classical evolution. The quantum 
NLO results are shown with black symbols for comparison. The initial conditions 
are characterized by low occupation numbers so that quantum effects become 
sizeable. One observes that in the quantum theory the damping rate is reduced 
compared to the classical theory. (All in units of rriR.) 


invariant dynamics can be explicitly reversed in practice for not too late times. 
We then compare to the corresponding NLO results of the quantum theory and 
check the classicality condition of section 14.31 

Here, we consider a system in d — 1 that is invariant under space translations 
and work in momentum space. Similar calculations have been performed also 
for d = 3 and for theories including fermions, which we refer to in the literature 
section 14.51 At to = 0, we choose the same initial conditions for the classical 
and the corresponding quantum theory: A Gaussian initial state with zero 
macroscopic field and a statistical propagator F(0,0;p) = [/ p (0) + l/2]/cap(0), 
where the initial f p ( 0) represents a peaked distribution around the momentum 
p = p ts as in section !?. 1.31 The initial mode energy is given by ca p (0) = \J p 2 + M 2 , 
where M is the one-loop renormalized mass in the presence of the nonequilibrium 
medium, determined from the corresponding one-loop gap equation as in (J3j33]). 
As a renormalization condition we choose the one-loop renormalized mass in 
vacuum rriR = M|y p ( 0 ) =0 as our scale. The results shown below are obtained 
using a fixed coupling constant A /m 2 R = 30. 

On the left of Fig. !T3l the classical-statistical propagator F(t,0]p = 0) is 
presented for three values of N. All other parameters are kept constant. The 
figure compares the time evolution using the 2PI 1/N expansion to NLO and 
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the classical-statistical simulation] 15 ! One observes that the approximate time 
evolution of the correlation function shows a rather good agreement with the full 
result even for small values of N. For N = 20 the exact and NLO evolution 
can hardly be distinguished. A very sensitive quantity for comparisons is the 
damping rate 7 , which is obtained from an exponential fit to the envelope of 
F(t,0\p = 0). The systematic convergence of the NLO and the Monte Carlo 
result as a function of 1/N can be observed from the right graph of Fig. [13] The 
accuracy of the description of far-from-equilibrium processes within the classical- 
statistical NLO approximation of the 2PI effective action is remarkable. 

The right graph of Fig. [13] also shows the damping rate from the quantum 
evolution, using the same initial conditions and parameters. One observes that 
the damping in the quantum theory differs and, in particular, is reduced compared 
to the classical result. The effective loss of details about the initial conditions 
takes more time for the quantum system than for the corresponding classical one. 
In the limit N —y 00 damping of the unequal-time correlation function goes to 
zero since the nonlocal part of the self-energies vanishes identically at LO large -N 
and scattering is absent. In this limit there is no difference between evolution in 
a quantum and classical-statistical held theory for same initial conditions. 

For finite N scattering is present and quantum and classical evolution differ 
in general. However, as discussed in section 14.31 the classical held approximation 
may be expected to become a reliable description for the quantum theory if 
the number of held quanta in each held mode is sufficiently high. We observe 
that increasing the initial particle number density leads to a convergence of 
quantum and classical time evolution at not too late times. In Fig. [14] (left) 
the time evolution of the equal-time correlation function F(t, t:p ) is shown for 
several momenta p and N = 10. Here the initial integrated particle density 
f dp/(2 tt) f p (0)/M = 1.2 is six times as high as in Fig. [13] At p — 2 p ts one 
hnds / 2 p ta ( 0 ) — 0.35 and a slightly larger value at this momentum of about 
~ 0.5 at later times shown. For these evolutions the classicality condition (14.1621) 
is therefore approximately fulfilled up to momenta p ~ 2 p ts , and one indeed 
observes from the left of Fig.[H]a rather good agreement of quantum and classical 
evolution in this range. For an estimate of the NLO truncation error we also give 
the full Monte Carlo result for IV = 10 showing a quantitative agreement with 
the classical NLO evolution during these times. 

From the left of Fig. [TT] one observes that the initially highly occupied modes 
“decay” as time proceeds and the low momentum modes become more and 
more populated. At late times the classical theory and the quantum theory 
approach their respective equilibrium distributions. Since classical and quantum 
thermal equilibrium are distinct, the classical and quantum time evolutions have 

15 The results presented below have been obtained from sampling 50000-80000 independent 
initial conditions to approximate the exact evolution of correlation functions. In general, far 
fewer configurations are required for d > 1 because of self-averaging. 
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Figure 14: Left: Nonequilibrium evolution of the equal-time two-point function 
F(t,t]p ) for N = 10 for various momenta p. One observes a good agreement 
between the full simulation (dashed) and the NLO classical result (full). The 
quantum evolution is shown with dotted lines. The integrated initial particle 
density is six times as high as in Fig. [T2J Right: A very sensitive quantity to 
study deviations is the time dependent inverse slope T(t,p ) defined in the text. 
When a Bose-Einstein distributed occupation number is approached, all modes 
get equal T(t,p) = T eq , as can be observed to high accuracy for the quantum 
evolution. For classical thermal equilibrium the defined inverse slope remains 
momentum dependent. 


to deviate at sufficiently late times. Figure fl4l shows the time dependent inverse 
slope parameter 

T(t,p) = ~fp(t)[f p (t) + 1] (Jj^j , (4.166) 

which has been introduced in section 13.2.51 to study the approach to a Bose- 
Einstein distribution]^] It employs the effective particle number f p {t) defined 
in (13.64p and mode energy e p (t) given by (I3.66p . Initially one observes a very 
different behavior of T(t,p) for the low and high momentum modes, indicating 
that the system is far from equilibrium. The quantum evolution approaches 
a Bose-Einstein distributed occupation number with a momentum-independent 
inverse slope T cq = 4.7 tur to very good accuracy. In contrast, in the classical 
theory the slope parameter remains momentum dependent since the classical 
dynamics does, of course, not reach a Bose-Einstein distribution. 

To see this in more detail we note that for a Bose-Einstein distribution, 
fp{e p ) — l/[exp(e p /T eq ) — 1], the inverse slope (I4.166P is independent of the mode 
energies and equal to the temperature T eq . During the nonequilibrium evolution 
effective thermalization can therefore be observed if T(t,p) becomes time and 

16 Note that dLog(/“ 1 (t) + 1 )/de p (t) = T -1 (t,p). 
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momentum independent, T(t,p) —> T eq . This is indeed seen on the left of Fig. [14] 
for the quantum system. If the system is approaching classical equilibrium at 
some temperature T c \ and is not too strongly coupled, the following behavior is 
expected. From the definition (13.64ft of f p (t) in terms of two-point functions, we 
expect to find approximately 


T(t,p)^T d (l~4/T“) , (4,167) 

i.e. a remaining momentum dependence with T(t,p ) < T(t,p') if e v > e'. Indeed, 
this is what one observes for the classical field theory result in Fig. [14] 

For a classical theory a very simple test for effective equilibration is available. 
An exact criterion can be obtained from the classical counterpart of the “KMS” 
condition for thermal equilibrium discussed in section 12.71 In coordinate space 
the classical equilibrium “KMS” condition reads 

(x-y) = -p { c i q) (x-y), (4.168) 

and in momentum space 

ddm = -i/ci^VS’A), Uk°) = §?■ (4.169) 

Differentiating equation (j4. 168j) with respect to y° at x° = y° = t gives 


T d dy° dx° cl 


(x-y) 


X 0—y0 = t 


9 (eq) 

W° Pd 


(x-y) 



(4.170) 


Combining this KMS relation with the equal-time condition (14. 141)1 for the 
spectral function leads to 

dtdt'F^ (t, t'\ x - y)\t=t' = F d 5(x - y). (4.171) 

In terms of the classical conjugate momentum fields 7T a (x) = d x oip a (x) this rep¬ 
resents the well-known equilibrium relation {n a (t, x.)^Tb(t, y))^ = T c i<5(x — y)J a6 . 
Out of equilibrium one can define an effective classical mode temperature 


T cl (t,p) = d t d t 'F cl (t,t';p)\ t=t '. (4.172) 


Effective classical equilibration is observed if T c \(t,p) becomes time and 
momentum independent, T c \(t,p) —> T c \. This is indeed the case for the results 
presented at sufficiently late times for given lattice regularization. 

Apart from the 2PI 1/N expansion, the late-time behavior can also be studied 
from the 2PI loop expansion. For the quantum theory this has been demonstrated 
in section 13.1.31 using two-loop self-energy corrections. We show that this can 
also be employed for the classical-statistical theory and use the corresponding 
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Figure 15: Nonequilibrium time evolution in 1 + 1 dimensions from the three- 
loop approximation of the classical-statistical 2PI effective action for one scalar 
held, N — 1. Shown is the effective mode temperature T c \(t,p). One observes 
the approach to classical equilibrium, T c i(t,p ) —> T c i, at sufficiently late times for 
given lattice regularization. 


loop approximation in the following to demonstrate the above statements about 
classical equilibration. 

In Fig. [15] the nonequilibrium evolution of the classical mode temperature 
T c i(t,p) is shown for various momentum modes in the (l + l)-dimensional classical 
scalar held theory for N — 1. The equations are solved by a lattice discretization 
with spatial lattice spacing m R a s = 0.4, time step a t /a s = 0.2, and N s = 24 
sites for A /m 2 R = 1. For the initial ensemble we take F c i(0, 0;p) = T 0 /(p 2 + m 2 R ) 
and d t dt'F c i(t,t'-,p)\t = t'=o = T 0 with T 0 /m R = 5. We have observed that at 
sufficiently late times the contributions from early times to the dynamics are 
effectively suppressed. This fact has been employed in Fig. [TS] to reach the 
very late times. One sees that at sufficiently late times the system relaxes 
towards classical equilibrium with a hnal temperature T c \/m R ~ 5.5. Though 
the typical classical equilibration times observed are substantially larger than 
the times required to approach thermal equilibrium in the respective quantum 
theory, this is not a regularization independent statement for the classical theory. 
In contrast to the quantum theory, statements about equilibration times are 
sensitive to the employed lattice regularization for the classical theory because of 
the Rayleigh-Jeans divergence. 
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5 Nonequilibrium instabilities 

5.1 Parametric resonance 

In classical mechanics parametric resonance is the phenomenon of resonant 
amplification of the amplitude of an oscillator having a time-dependent periodic 
frequency. In the context of quantum held theory a similar phenomenon describes 
the amplification of quantum fluctuations, which can be interpreted as particle 
production. It provides an important building block for our understanding of the 
(pre)heating of the early universe at the end of an inhationary period, and may 
also be operative at certain stages in relativistic heavy-ion collision experiments. 
The example of parametric resonance can lead to nonperturbative phenomena 
such as strong turbulence and Bose condensation far from equilibrium, even in 
the presence of arbitrarily small couplings, which will be discussed in section [6j 
Here we will consider the phenomenon as a “paradigm” for far-from-equilibrium 
dynamics following nonequilibrium instabilities. Much of the nonlinear physics 
turns out to be universal and thus independent of the details of the underlying 
mechanism that triggers the instability. 

We recall hrst the classical mechanics example of resonant amplitude growth 
for an oscillator with time-dependent periodic frequency. A physical realization 
of this situation is a pendulum with a periodically changing length as displayed: 

periodic 



In the linear regime, the amplitude y{t) is described by the second-order 
differential equation 

y(t)+u 2 (t)y(t) = 0 (5.173) 

with periodic uj(t + AT) = (j(t ) of period AT. Since the equation is invariant 
under t —> t + AT one expects periodic solutions 

y(t + AT) = Ay{t ). (5.174) 

Writing the time-independent amplitude as 

A = e aAT (5.175) 

in terms of the so-called Floquet index cc, these solutions can be expressed as 

y(t) = e at n (t) (5.176) 
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(5.177) 


with periodic II (t + AT) = II (i). This can be directly verified since 
y(t + AT) = e “ (t+AT) U(t + AT) = A e Qt U(t) . 

y{t) 

With y(t) also y(—t) is a solution and one concludes from (15. 176jl that for real 
Floquet index a ^ 0 there is an instability characterized by an exponential 
growth. 

In contrast to this mechanics example, in closed systems described by 
quantum field theory there will be no external periodic source. Below we will 
see that a large coherent field amplitude coupled to its own quantum fluctuations 
can trigger the phenomenon of parametric resonance. Mathematically, however, 
important aspects are very similar to the above classical example for sufficiently 
early times. There is even a precise mapping: The mechanical oscillator 
amplitude y plays the role of the statistical two-point function F in quantum 
field theory, and the periodic c o 2 {t) plays the role of an effective mass term 
M 2 (0(f)) whose time dependence is induced by an oscillating macroscopic field 
0(t). Simple linear approximations to the problem turn out to be mathematically 
equivalent to the above mechanics example. Accordingly, well-known Lame-type 
solutions of the mechanics problem will also play a role in the quantum field 
theory study. Substantial deviations do, however, quickly set in with important 
non-linear effects. 

5.1.1 Linearized classical analysis 

We consider the 0(N ) symmetric scalar field theory with classical action (I2.61|h 
Since the phenomenon of parametric resonance is essentially classical, for the 
purpose of comparison with the quantum treatment we start with a linearized 
classical analysis. The iV-component field is written as 

4>a{x) = 4>Sal + 8<t>a(x) , (5.178) 

where we will consider a spatially homogeneous but time-dependent “background 
field”, i.e. 0 = 0(x°). The classical equations of motion are obtained from the 
stationarity of the action, and in the following we will consider them in linear 
approximation in the “fluctuations” 50 a (x). These equation will also be seen to 
correspond to the evolution equations in the corresponding quantum field theory 
in the limit where all nonlinear or loop corrections are neglected. 

We denote the fluctuations in the “longitudinal” field direction as 50y(x) = 
(x) and in the “transverse” direction as 50 _l(x) = 50 a> i(x). The square of 
the fields appearing in the action (12.611) then reads 

0 a (x)0 a (x) = 0 2 + 2050||(x) + 50||(x) 2 + (N — l)50j_(x) 2 , (5.179) 
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whereas the quartic interaction term involves 

[4>a(x)4>a{x)} 2 = 0 4 + 40 3 8(j)\\ {x) + 20 2 [3 J0|| {x) 2 + (IV - 1)<5</>_l(x) 2 ] + (!) ((50 3 ) . 

(5.180) 

Here we expanded the quartic term neglecting cubic and quartic powers of the 
fluctuations, since this is sufficient to obtain the evolution equations to linear 
order. 

The classical held equation for the longitudinal component is given by the 
stationarity condition 5S[(j) + 5(j)]/65(j)\\(x) = 0. Using (15.1791) and (j5. 180j) one 
observes that to lowest order in an expansion in powers of 8^>\\(x) it leads to the 
background held equation 

/ d 2 \ \ 

+ m2+ 6N ^ Xo ) = 0 • ( 5 - 181 ) 

To next-to-leading order in powers of 8<f>\\(x) the same stationarity condition gives 

(□x + m 2 + ^ (fi 2 (xo)^J 8</)\\(x) = 0 . (5.182) 

For the transverse components one hnds from 8S[(f> + 8<fi]/8 8(f)±(x) = 0 the 
equation 

(□.X + m2 + 0 2 ( x o)^ 5<j)±(x) = 0 . (5.183) 

In the following we will choose m — 0 as motivated in section 11.31 and denote 
t = x o- It is convenient to introduce the rescaled background held 

°{t) = (5.184) 

We will consider the relevant case of a parametrically large initial held amplitude, 
4>(t = 0) ~ yj6N/X with A -C 1, such that the initial cr 0 = a(t = 0) is of order 
one. In terms of the rescaled held the evolution equation (15. 181ft reads 

cr(t) + & 3 (t) = 0. (5.185) 

The solution of this equation can be given in terms of Jacobi elliptic functions. 
For an initial amplitude Uo and derivative a(t — 0) = 0 one has 

a(t) = cr 0 cn((j 0 t). (5.186) 

Here the Jacobi cosine cn(z) = cn (z;n = 1/2) is a doubly periodic function 
in z with periods 4 K n and 4 iKi_ n , where K n = dO/\/l — n sin 2 0 is the 
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complete elliptic integral of the first kind. We denote the characteristic frequency 
of the oscillations of a(t) as 


= 77 ^ - 0.847 cro, (5.187) 

2 A1/2 

where Ki/ 2 — 1.854. Since the function will enter quadratically the equations for 
fluctuations, it will also be important that it is a quasi-periodic function with 
period cn(z + 2K 1 / 2 ) = —cn(z). Averaged over one period, the square of the field 
amplitude gives a 2 = Jq Ki/2 dt a 2 (t)/(2K 1 / 2 ) — 0.457 ctq. 

In order to study the fluctuations, we will consider the products 

F\\{x,y) = 6<j)\\(x)5(j)\\(y) , F ± (x, y) = 5<j>±(y ), (5.188) 

which will be the relevant quantities to compare to in the quantum 
treatment below. For spatially homogeneous systems they only depend on 
the relative spatial coordinates, F^±(x,y) = F|| ) j_(x°, y°; x — y). According 
to (15.182)1 and (15.183)1 their spatial Fourier transforms, F^±(t,t'',p) = 

f d d sexp(—ips)Fjj t ±(t,t';s), obey in this linearized classical approach 

[d 2 + P 2 + 3Up cn 2 (<j 0 t)] F\\(t, t'\ p) = 0 , (5.189) 

[d 2 + p 2 + ag cn 2 (cr 0 t)] F±(t, t'] p) = 0 . (5.190) 

We will consider initial conditions characterizing a pure-state initial density 

matrix, where we take d t F^±(t, 0; p)|t= 0 = 0 and d t d t >F^±(t, F; p)| t=t / =0 = 

F|P_{^(0,0; p)/4 = o;|| i x(p)/2 in accordance with the results of section (12.4p . 

The frequency of the longitudinal modes is o;||(p) = \/p 2 + 3crg, and cax(p) = 
\/p 2 + (Tq for the transverse modes. 

Of course, in these linear equations the two-point functions can be factorized 
as products of momentum-dependent single-time functions /||,j_(t; p) with 

t p) = ^ [/||,±(t; p)/||,±(i / ; p) + p)f\\,±(t; p)] , (5.191) 

where ,fu ± {t] p) denotes the complex conjugate of /||,x(t;p). In terms of these 
so-called mode functions the equations of motion read, e.g., for the transverse 
modes 


[d 2 + p 2 + o-q cn 2 (<7 0 f)] f ± (t] p) = 0 . (5.192) 

Up to an overall arbitrary phase, the above initial conditions for the two-point 
functions translate into 

f±( 0,p) = 1 /a/ 2 ^±(p) , dtf±(t,p)\ t=0 = -i >J U± !^ (5.193) 
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and equivalently for f\\. For given momentum p, equation (15. 192)1 is 
mathematically equivalent to the classical mechanics oscillator described above. 
Before we summarize the analytical solution of the Lame-type equation (15. 19211 . 
we discuss some general properties. In order to be able to observe a significant 
resonance, the frequency ujq of the field a(t) and the characteristic frequency of 
the fluctuations /n (f; p) should be similar. Therefore, for transverse modes the 
dominant resonances should fulfill the approximate condition 

w±(t, p) = Vp 2 + — \Jp 2 + ^o/ 2 ~ w o, (5.194) 

where we replaced the held square by its approximate time average. Since 
cj 0 ~ cro, this condition could be met for not too high momenta p 2 < crg/2. 
In contrast, for longitudinal modes the same analysis leads to the approximate 
condition ^/p 2 + 3 <Tq/2 ~ ujq which even for vanishing momenta is not fulfilled 
such that only subdominant resonances may be expected. The results of this 
qualitative analysis are indeed validated by the analytic solutions of the Lame- 
type equation which we turn to next. 

In the following, we concentrate on the dominant transverse modes to analyze 
the instability dynamics at early times. We choose U p (t) and U p (—t ) as an 
independent set of solutions of (15.1921) such that p) is a linear combination of 
those. Here we will be interested in the unstable modes which show exponential 
amplification and state the corresponding growth rates. Similar to the above 
classical oscillator example, we consider a Floquet analysis of the periodic solution 

U p (t + AT) = e apAT U p (t) , (5.195) 

where the time-independent index a p is a function of spatial momentum. 
Correspondingly, we can write 


U p (t) = e apt Tl p (t) (5.196) 

with the periodic function n p (t + AT) = n p (i) of period AT = 2K 1 / 2 a 0 . Again, 
this can be directly verified with U p [t + AT) = exp{a p (t + AT)}n(i + AT) = 
exp{a p AT}U p (t). The solutions of Lame equations can be found in textbooks 
on differential equations and we only give here the relevant properties of the 
Floquet index a p for (15.1921) . In accordance with our qualitative discussion above, 
resonant amplification of modes can be found for not too large momenta. More 
precisely, for p 2 > <Tq/ 2 purely oscillating solutions are obtained. In contrast, for 
0 < p 2 < ctq/ 2 the Floquet index 


a p = iujo — 'fp (5.197) 

has a real part and an imaginary part cuo, given by the characteristic frequency 
(15.187[) . The real part implies exponential solutions with a growth or decay rate 
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for either U v {t) or U p (—t ) given by 


7 P = cr 0 Z [ cn 1 


'2p 2 

Cn 


(5.198) 


Here, the Jacobi zeta function Z(x) has the series expansion 


Z(x) = 


2vr 
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Y — 

!/ 2 n=l 1 e 


2vr 
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e _7r sin 


1/2 


2ri7r 


7T;C 


sin 


T171 X 
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1/2 


AT 


1/2 


(5.199) 


where the latter approximation gives an excellent description for our purposes. 
Taking into account that cn” 1 (0) = AT L / 2 and cn _1 (l) = 0, we note that the rate 
7 p vanishes at the boundaries p 2 = 0 and P 2 = Cq/ 2 of the instability band. 

In order to get a further analytic understanding we may use the somewhat 
crude approximation 

9 K, 

arccos(x), (5.200) 


-i/ \ 2A 1/2 

cn (x ~ — 


7T 


which in our case will be good at the few percent level, and the identity 


sin (2 arccos [y/x )) = 2a/2x(1 — x). 
Then the rate (j5. 198j) can be expressed as 


(5.201) 


4vt(Jo 

7 P — —— e 
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1/2 


2p 2 

Cn 


1 - 


2P 2 

Cn 


(5.202) 


and the maximally amplified growth rate 70 = 7(p 2 = Pq) occurs for Pq = Cq/4 
with 

70 ~ ^ a ° e -7r ~ 0.146 cr 0 . (5.203) 

-ft 1/2 

Therefore, there is a separation of scales between the characteristic frequency of 
oscillations, (u 0 , and the characteristic growth rate, 7 0 : 


70 < w 0 • (5.204) 

When analyzing the growth in the nonlinear regime below, we will often exploit 
this separation of scales by considering suitable time averages to smooth out the 
rapid oscillations. 
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5.2 Nonlinear regime: secondary instabilities 

5.2.1 Dynamical power counting 

Above we obtained the set of evolution equations (15.1851) and (15.1901) for the field 
4>(t) = a (t) y/6N/\ and the fluctuations F^± from linearized classical dynamics. 
According to section 12.61 the very same set of equations can also be obtained 
starting from the quantum evolution of the rescaled field expectation value and 
the statistical two-point functions by neglecting all loop corrections: Setting E|j_, 

and to zero in the field equation (12.1171) and the equations for the two- 
point functions (12.1111) leads precisely to (15.185ft and (15.190ft . We note that the 
evolution equation for the spectral function ^ decouples from the rest in this 
linear approximation. It is important to point out that this linearized classical 
approximation is only valid at early times for A <C 1. For strong coupling loop 
corrections cannot be neglected in this case. 

In turn, we may use the quantum evolution equations (12.117ft and (12.111ft to 
determine the range of validity of the above classical approach. Parametrically 
the initial statistical propagators are of order one, i.e. 

iV(0,°;P) ~ 0(N°X 0 ), (5.205) 

and the initial field is 

0(0) ~ O (N 1/2 X~ 1/2 ) . (5.206) 

As a consequence, all loop corrections are initially suppressed by powers of the 
coupling constant A C 1 and the classical approximation represents an accurate 
description of the early quantum dynamics. However, parametric resonance leads 
to an exponential amplification of modes, where the dominant growth behavior 
is parametrically given by 

F ± (t,t'\p) ~ e lp{t+t,) . (5.207) 

Accordingly, the maximally amplified mode F±(t, t,\ p 0 ) grows with rate 2yo. As a 
consequence of the exponential amplification of the statistical propagator, there 
is a characteristic time when F±(t, t: p 0 ) turns out to be no longer parametrically 
of order one. This is the time where the linearized classical approximation 
breaks down. To make further analytical progress, the nonlinear regime is most 
efficiently described using the evolution equations of section 12.61 This allows us 
to study what happens when loop corrections start to become relevant. Indeed, 
we will find that there is a well separated hierarchy of times, where different loop 
integrals start to become of order one at corresponding separate times. At some 
point an infinite number of diagrams with an arbitrary number of loops becomes 
of order one and no power counting based on a small coupling can be performed. 
We will address this later stage in a subsequent section on nonthermal fixed points 
using nonperturbative large-At techniques. 
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A general loop contribution to the evolution equations contains powers of 
A, the held 0, the propagators Fy j_ and the spectral functions p\\ ,j_. Here it is 
important to note that the “weight” of the spectral functions in loop integrals 
remains parametrically of order one at all times as encoded in the equal-time 
commutation relations (I2.98jh It is also important to take into account the 
fact that transverse fluctuations (F±) exhibit the dominant growth in the linear 
regime. Consequently, contributions containing more transverse propagators 
(Fj_) can become important earlier than those diagrams containing longitudinal 
propagators (Fj|) instead. For instance, an expression containing powers \ n F™(j) 21 
with integers n, m and l may be expected to give sizable corrections to the 
linearized evolution equations once F± ~ 1/A ( n ~ l )/ m for typical momenta. Here 
n yields the suppression factor from the coupling constant, whereas m introduces 
the enhancement due to large fluctuations for typical momenta and / due to 
a large macroscopic held. The power counting can become more involved as 
time proceeds, and it is remarkable that one can indeed identify a sequence of 
characteristic time scales with corresponding growth rates. 

Before we consider the calculations in more detail, the situation is 
schematically summarized in Fig. [161 The time when Fj_(t,t;po) ~ O(N 0 X -1 / 2 ) 
is denoted by t — t\. At this time the one-loop diagram with two held insertions 
indicated by crosses as depicted in Fig. [[6] will give a contribution of order one 
to the evolution equation for F||(£, £;p). For instance, the two powers of the 
coupling coming from the vertices of that diagram are canceled by the held 
amplitudes and by propagator lines associated to the amplified F±(t, t'\ p 0 ). 
Similarly, at the time t = t 2 the maximally amplified transverse propagator mode 
has grown to F±(t, p 0 ) rs_/ £>(AT/ 3 A- 2 /3). As a consequence, the “setting sun” 

diagram in Fig. [12] becomes of order one and is therefore of the same order as 
the classical contributions. Though the loop corrections become of order one 
later than the initial time, they induce amplification rates that are multiples of 
the rate 2y 0 which lead to a very rapid growth of modes in a wide momentum 
range. Finally, when the fluctuations have grown nonperturbatively large with 
F±(t,t'] p 0 ) O(N 0 A x ) any loop correction will no longer be suppressed by 
powers of the small coupling A. In this case the nonperturbative 1/N expansion 
of the 2PI effective action, which will be discussed later, becomes of crucial 
importance for a quantitative description of the dynamics. 


5.2.2 Nonlinear amplification 

The O(A 0 ) approximation for longitudinal modes breaks down at the time 

t ~ t' = h : F±(t, t'; p 0 ) ~ e>(N 0 A" 1/2 ) . (5.208) 


This can be derived from the 0( A) evolution equations to which one-loop self¬ 
energies contribute, which diagrammatically are given by O ,~Q~. The 
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Figure 16: Schematic overview of the characteristic time scales and the respective 
relevant diagrammatic contributions. 


approximate evolution equation reads 


where 


(dt + p 2 + M 2 (t) + 3a 2 (t))F||(t,t'; p) ~ 

2M '3n 1> g(t) { j‘ (*■ t "’ p)^ii (*".p) 

-\f d<V(t")n£(«,t"i p)pii(«", i';p)| = ms. 


M 2 {t) 

T\\At) 


m 2 + 


A 


61V L 
A d 3 p 


(2tt) 


3T,|(t) + (IV — l)Tj_(t) 
*ii,±(M;p) • 


(5.209) 

(5.210) 

(5.211) 
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The “tadpole” contributions Tj| and T± from the longitudinal and transverse 
propagator components, respectively, are regularized by some A>j ) 0 whose 
precise value is irrelevant as long as the integral is dominated by momenta 
much smaller than the cutoff. Here we have abbreviated n^(t, t";p) = 
I (27t) 3 F±(t, t"; p — q) X±(t, t"; q) with X = {F, p}, and we used that F]_ 
p^O One observes that indeed for F± ~ O(N 0 X~ 1 ^ 2 ) the RHS of (I5.209P becomes 
~ 0(1) and cannot be neglected. 

In order to make analytical progress, one has to evaluate the “memory 
integrals” in the above equation. This is dramatically simplified by the 
fact that the integral is at this stage approximately local in time, since the 
exponential growth lets the latest-time contributions dominate the integral. For 
the approximate evaluation of the memory integrals we consider 



(5.212) 


with c ~ 0(1). As long as we are only interested in exponential growth rates 
and characteristic time scales, the value for the constant c will turn out to be 
irrelevant to “leading-log” accuracy. We then perform a Taylor expansion around 
the latest time t (£'), 


P||,±(M";P) ^ d t "Pi±(t,t"-,P)\t=t"(t"-t) = (t-t"), 

Fi±(t,t";p) ~ F l± (t,t;p), (5.213) 


where we have used the equal-time commutation relations (12.981) . With these 
approximations the RHS of (15.209(1 can be evaluated as: 


RHS ~ A a\t) 


2^(iV-l) 


+ Xa(t)a(t) -3 


2 T x (t) F\\(t,t'-, p) (mass term) (5.214) 

l c2 


Un 


6N 


-ni(t,t';p) (source term) (5.215) 


The first term is a contribution to the effective mass, whereas the second term 
represents a source. Note that both the “tadpole” mass term and the above 
correction to this mass are of the same order in A, however, with opposite 
sign. To evaluate the momentum integrals, we use a saddle point approximation 
around the dominant p ~ p 0 , valid for t,tf y^ 1 , with F±(t,t',p) ~ 
F±(t, t', p 0 ) exp[— | 7 o|(t + t')(p — Po) 2 /2]@ From this one obtains for the above 
mass term: 


A T ± (t) 


, Pi F ± (t,t] p 0 ) 
2(tt 3 | T "| ty /2 ■ 


(5.216) 


17 This corresponds to the classicality condition of section 1Q1 
18 Here y(p) ~ 70 + 7o (? - Po) 2 /2 with 7 ^ < 0. 
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The result can be used to obtain an estimate at what time t this loop correction 
becomes an important contribution to the evolution equation. Note that to this 
order in A it is correct to use F±(t,t'] p 0 ) ~ e 7°(<+<') on the RHS of (15.2161) . The 
condition XT±(t = i 3 ) ~ 0(1) can then be written for A 1 as: 


t 3 ~- In A 1 

2 7o 


(5.217) 


The same saddle point approximation can be performed to evaluate the source 
term (15. 215ft : 


An£(t,t';0) 


N Po F±(t,tf ]p 0 ) 

~ 4(7T 3 |7"|(f+ f')) 1/2 ' 


(5.218) 


Here we only wrote the source term for p = 0 where it has its maximum, although 
it affects all modes with p < 2 p 0 . Again this can be used to estimate the time 


t — t\ at which An^ ~ 0( 1): 

t^\t 3 


(5.219) 


One arrives at the important conclusion that the source term (j5.215j) becomes 
earlier of order one than the mass term (15.214)1 : For H < f < f 3 the source term 
dominates the dynamics! Using these estimates in (j5.209jl one finds that the 
longitudinal modes with 0 < p < 2po get amplified with twice the rate 2yo: 


F||(f,f;p) ~ A F]_(t, t; p 0 ) ~ A e 470< . (5.220) 

Though the non-linear contributions start later, they grow twice as fast. The 
analytical estimates for t\ and rates agree well with the numerical solution^ of 
the evolution equations without memory expansion as shown in Fig. [17j where 
we plot the effective particle number distribution 

P) = \jFi ± (t,t-,p)d t d t/ F l± (t,t';p)\ t=t , - i (5.221) 


and Mq = M 2 (t = 0). 

A similar analysis can be made for the transverse modes. Beyond the Lame- 
type C?(A°) description, the evolution equation for F± receives contributions from 
the feed-back of the longitudinal modes at 0( A) as well as from the amplified 
transverse modes at 0( A 2 ). They represent source terms in the evolution equation 

19 The plots are obtained for the IV-component field theory at NLO in the 1/N expansion 
of the 2PI effective action with nonzero mass parameter m. The latter is not relevant for our 
purposes. 
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Figure 17: Effective particle number for the longitudinal inodes as a function 
of time for various momenta. For t > t\ the nonlinear corrections trigger an 
exponential growth with rate 470 for p < 2po- The thick line corresponds to a 
mode in the parametric resonance band, and the long-dashed line for a similar 
one outside the band. The resonant amplification is quickly dominated by source- 
induced amplification. 


for F±(t, t'\ p) which are both parametrically of the form ~ \ 2 F^_/N as is depicted 
below: 



Following along the lines of the above paragraph and using (15. 217ft this leads to 
the characteristic time t — t 2 at which these source terms become of order one: 


2 

In N 

t2 ^ 2 h + 

670 


For t ~ t' ~ t 2 the dominant transverse mode has grown to 


Fj.(*,*';P o) ~ 0(N 143 A- 2 / 3 ). 


(5.222) 


(5.223) 


Correspondingly, for t 2 < t < £3 one finds a large growth rate ~ 670 in 
a momentum range 0 < p < 3po, in agreement with the numerical results 
shown in Fig. [18l In this time range the longitudinal modes exhibit an 
enhanced amplification as well (cf. Fig. [T71) . It is important to realize that the 
phenomenon of source-induced amplification repeats itself: the newly amplified 
modes, together with the primarily amplified ones, act as a source for other 
modes, and so on. In this way, even higher growth rates of multiples of 70 can 
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Figure 18: Effective particle number for the transverse modes as a function of 
time for various momenta p < 5po- At early times, modes with p ~ p 0 are 
exponentially amplified with a rate 2yo. Due to nonlinearities, one observes 
subsequently an enhanced growth with rate 670 for a broad momentum range. 


be observed and the amplification rapidly propagates towards higher momentum 
modes. We note that this regime is present if t 2 < t 3 , i.e. as long as N < A -1 . 

Around t < t 3 is the earliest time when sizeable corrections to the maximally 
amplified mode F±(t,t; po) and to the field appear. Around t 3 there are 
corrections of order one which come from diagrams with an arbitrary number 
of loops. As a consequence, the dynamics is no longer characterized in terms of 
the small expansion parameter A. This is the nonperturbative regime, where the 
1/N expansion may be used to describe also the subsequent evolution governed 
by a nonthermal fixed point as is explained in section |6] 

The accuracy of the above description based on the 2PI effective action can 
be tested using numerical simulations on a space-time lattice. This exploits 
the fact that classical-statistical field theory descriptions of the dynamics have 
an overlapping range of validity with the underlying quantum field theory, as 
described in section [2 As an example, Fig. [T9] shows F±(t,t, p) for times t = 10, 
t = 40 in the nonlinear regime, and t — 90 in the nonperturbative regime. The 
solid lines give the results for the quantum evolution for A = 0.01 and IV = 4 
by solving the NLO equations from the 2PI effective action of section 13.2.31 in 
d = 3 numerically. For comparison, the dashed lines give the same quantity as 
obtained from simulations of the corresponding classical-statistical field theory on 
a lattice with same initial conditions. The level of agreement between the different 
results is remarkable. Quantum fluctuations are expected to be suppressed if the 
classicality condition F 2 p 2 of section l4~3l is fulhllcd, which is very well realized 
for typical F ~ 1/A with p being order unity. In turn, contributions beyond NLO 
in the 1/N expansion seem to play an inferior role even for the nonperturbative 
regime in this case. 
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Figure 19: Two-point function F±(t,t',p) as a function of momentum |p| for 
three different times. The quantum evolution (solid) from the 2PI 1/N expansion 
to NLO and the full classical-statistical simulation using Monte Carlo sampling 
(dashed) agree remarkably well. 
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6 Nonthermal fixed points and turbulence 

From the example of parametric resonance in section [5] we have seen that once the 
exponential growth of fluctuations stops, the systems slows down considerably. 
At this stage there is an infinite number of loop corrections of order unity. The 
coupling “drops out” of the problem and no power counting in terms of a small 
coupling parameter can be given. In this section we want to study the physics 
of this slow evolution for the examples of relativistic TV-component scalar field 
theory as well as nonrelativistic (Gross-Pitaevskii) held theory. They serve here 
as a paradigm to investigate universal behavior of isolated many-body systems 
far from equilibrium, which is relevant for a wide range of applications from high- 
energy particle physics to ultracold quantum gases as pointed out in section |T| 
The universality is based on the existence of nonthermal fixed, points, which 
represent nonequilibrium attractor solutions with self-similar scaling behavior. 
The corresponding dynamic universality classes turn out to be remarkably large, 
encompassing both relativistic as well as nonrelativistic quantum and classical 
systems. 

We start with a standard introduction to the phenomenon of wave turbulence 
in section 16.11 where perturbative kinetic theory can describe the stationary 
transport of energy from a cascade towards higher momenta. Perturbative kinetic 
theory has also been employed in the literature to describe infrared phenomena 
such as Bose condensation, however, such an approach neglects important vertex 
corrections since the large occupancies at low momenta lead to strongly nonlinear 
dynamics. Here we apply the vertex-resummed kinetic theory based on the 2PI 
1/N expansion to NLO, which is described in section 13.3.21 This allows us to gain 
analytic understanding of the formation of a dual cascade and the phenomenon 
of far-from-equilibrium Bose condensation. 

However, isolated systems out of equilibrium have no external driving forces 
that could realize stationary transport solutions. Instead, the transport in 
isolated systems is described in terms of the more general notion of a self-similar 
evolution, where the physics is described in terms of universal scaling exponents 
and scaling functions as outlined in section [H We can use the vertex-resummed 
kinetic theory for an analytic description of the phenomena, which we present 
in section 16.21 Because of the typical high occupancies, this can be compared to 
results from classical-statistical lattice simulations using the methods of sectional 

6.1 Stationary transport of conserved charges 

While many aspects of stationary transport of conserved charges associated to 
turbulence have long reached textbook level, there is still rather little known 
about turbulent behavior in nonperturbative regimes of quantum field theories. 
Here the strong interest is also driven by related questions concerning the 
dynamics of relativistic heavy-ion collisions. In this section we consider stationary 
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solutions of transport equations with a net flux of energy and particles across 
momentum scales, thus violating detailed balance. We start with standard 
perturbative results relevant for high momenta in section 16.1.11 and consider the 
nonperturbative infrared regime in section 16.1.21 


6.1.1 Perturbative regime: Weak wave turbulence 

In perturbative kinetic theory, when two particles scatter into two particles, the 
time evolution of the distribution function f p (t ) for a spatially homogeneous 
system is given by 

= C 2 " 2 [/] (t, p). (6.224) 

The collision integral is of the form 


C^ 2 [/](p) = J dfl 2 ^ 2 (p, 1, q, r) [(/ p + 1)(/! + l)/q/ r - /p/l(/q + l)(/r + 1)] , 

(6.225) 

where we suppress the global time dependence to ease the notation. The details 
of the model enter J dfi 2 ^ 2 (p, 1, q, r), which for the example of the relativistic 
IV-component scalar field theory with quartic A/(4!lV)-interaction reads: 


dft ** (p, 1, q, r) = A' 


2 N + 2 
61V 2 


(2n) d+1 <5(p + 1 - q - r) 


' lqr 


X ^(tJp + UJ\ — CJq — U r ) 


2u p 2ui2uq2u r 


(6.226) 


with ojp = y/p 2 + m 2 and f = f d d p/(27r) d . This kinetic equation can be 
obtained from quantum field theory using the two-loop self-energies (j3.24j) and 
(13.25(1 for the lowest-order gradient expansion along the lines of section l3.3.21 and 
taking the on-shell spectral function (I3.36P of the free theory. The expression 
(I6.225P with (16. 226 p for the collision integral represents a standard Boltzmann 
equation for a gas of relativistic particles. 

Clearly, this approximation cannot be used if the occupation numbers per 
mode become too large such that higher loop-corrections become sizeable. 
Parametrically, for a weak coupling A a necessary condition for its validity is 
f p <C 1/A since otherwise any loop-order contributes significantly as explained in 
detail in section [5j On the other hand, the phenomenon of weak wave turbulence 
is expected for not too small occupation numbers per mode as is explained in the 
following. For the corresponding regime 1 <C f p 1/A one may use the above 
Boltzmann equation, which approximately becomes 

-vjp- - J d0 2 ^ 2 (p, 1 , q, r) [(fp + /i)/ q / r - /p/i(/ q + / r )] • (6.227) 
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For number conserving 2 -H- 2 scatterings, apart from the energy density e also 
the total particle number density n is conserved, which are given by 


e = [ w p/p > n — j fp- (6-228) 

J p J p 

The fact that they are conserved may be described by a continuity equation in 
momentum space, such as 

^K/ P ) + V p .j p = 0 (6.229) 

for energy conservation. Similarly, particle number conservation is described by 
formally replacing co p —> 1 in the above equation and a corresponding substitution 
of the flux density. For the isotropic situation we can consider the energy flux 
A{k) through a momentum sphere of radius k. Then only the radial component 
of the flux density j p is nonvanishing and 



' Ok 


jp ■ dAp 


(2v r) d A{k). 


(6.230) 


Since in this approximation u p is constant in time, we can write with the help of 
(16.2291) and the kinetic equation (16.2241) : 


m 


i 

2 d 7r d / 2 T(d/2 + 1) 


d p |p 


d -1 


ujpC 


2-H-2 


(P) 


(6.231) 


Stationary wave turbulence is characterized by a scale-independent flux A(k ), for 
which the respective integral does not depend on the integration limit k. To this 
end, we consider scaling solutions 


fp 


7 


sp 


CJp 


S 1 CO 


S p > 


(6.232) 


with occupation number exponent k and assuming a linear dispersion relation 
relevant for momenta |p| m. Since the physics is scale invariant, we can 
choose s = 1/|p| such that f p = |p| -K /i and uj p = |p|cui. 

Using these scaling properties, one obtains for the collision integral (16.226(1 
and (16.227(1 of the theory with quartic self-interaction: 

U 2 ^ 2 (p) = s- /i4 U 2o2 (sp), (6.233) 


where the scaling exponent is given by 


p 4 = (3d - 4) - (d + 1) - 3k = 2d - 5 - 3k . (6.234) 


The first term in brackets comes from the scaling of the measure, the second from 
energy-momentum conservation for two-to-two scattering and the third from the 
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three factors of the distribution function appearing in (16.2271) . Apart from the 
4-vertex interaction considered, it will be relevant to investigate also scattering 
in the presence of a non-vanishing field expectation value such that an effective 
3-vertex appears according to section 12.51 

To keep the discussion more general, we may write for the scaling behavior of 
a generic collision term 

C(p) = |p| w C(l) (6.235) 

in terms of the scaling exponent pi for /-vertex scattering processes. In the general 
case of an /-vertex one obtains along these lines 


Hi= (/ — 2 )d — (/ + 1) — (/ — 1)k . 

For the scaling properties of the energy flux we can then write 

1 


A(k) = 


dp|p| d+w uqC'(l). 


2 d n d / 2 T(d/2 + 1) 

If the exponent in the integrand is nonvanishing, the integral gives 

fii 


m 


■uiC{ 1 ). 


d + 1 + fii 

Thus, scale invariance may be obtained for 

d + 1 + fii = 0 . 

This gives the scaling exponent for the perturbative 


relativistic energy cascade: k = d — 


l-l 


(6.236) 


(6.237) 


(6.238) 

(6.239) 

(6.240) 


One observes that stationary turbulence requires in this case the existence of the 

limit 

C(l) 

lim —- = const ^ 0, (6.241) 

d+i+n-x) d + 1 + in 

such that the collision integral must have a corresponding zero of first degree. 
Similarly, starting from the continuity equation for particle number one can study 
stationary turbulence associated to particle number conservation. This leads to 
the perturbative 


relativistic particle cascade: k 


/ + 

l-l' 


(6.242) 


Accordingly, for quartic self-interactions we get k = d— 4/3 for the energy cascade 
and k = d — 5/3 for the particle cascade. In the presence of a coherent held, when 
a 3-vertex can become relevant, the associated scaling exponents are k = d — 3/2 
for the energy cascade and n = d — 2 for the particle cascade. 
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6.1.2 Nonperturbative regime: Strong turbulence 

The above perturbative description is expected to become invalid at sufficiently 
low momenta. In particular, the occupation numbers / p ~ |p|” K for ft > 0 
would grow nonperturbatively large in the infrared such that the approximation 
(16.2251) becomes questionable. This concerns, for instance, the relevant case 
of d = 3 for the description of the early universe dynamics outlined in 
section 11.31 To understand where the picture of weak wave turbulence breaks 
down and to compute the properties of the infrared regime, we have to consider 
nonperturbative approximations. For this purpose, we employ the vertex- 
resunnned kinetic theory of section 13.3.21 based on the expansion of the 2PI 
effective action in the number of field components to NLO. 

For a relativistic theory with dispersion (Up = y/ p 2 + m 2 , the scaling 
assumption (I6.232|) should be valid for sufficiently high momenta |p| m such 
that the dispersion is approximately linear with u p ~ |p|. However, this is more 
involved at low momenta if a mass gap exists. An effective mass gap is typically 
expected because of medium effects even if the mass parameter in the Lagrangian 
is set to zero. In that case, the infrared modes behave effectively nonrelativistic 
as explained in section 12.81 In the following, we will analyze the two cases of a 
relativistic theory without mass gap and a nonrelativistic theory separately for 
comparison. 

Following similar lines as in section 16.1.11 we first look for relativistic scaling 
solutions. To be able to cope with occupancies of order ~ 1/A, we replace the 
perturbative collision term (I6.227P by the vertex-resummed expression (I3.106p . 
which can be approximated for f p ^> 1 accordingly by 


<^ NLO [/](p) - / dfi NLU [/]0,/,g,r) [(/p + fi)f q f r - f p fi(f q + fr)} , (6.243) 


%NLO 


with 


' jr'iNLO/ i \ A 2 r°° dp°dl 0 dq°dr 0 f 

d!J "ml (27r) 4 -( d+1 ) l v Hp+l - g ~ r) 

x pppipqpr [v e s(p + l) + v eS (p - q) + v eS (p - r)] (6.244) 


We emphasize that this still involves integration over frequencies as well as 
spatial momenta since no free-held form for the spectral function p p = p(p°, p) 
is used so far, and f p = /(p°, p). A crucial difference to the perturbative kinetic 
equation is the appearance of the vertex function v e s(p°,p) given by (I3.100p . 
which encodes the emergence of a momentum-dependent effective coupling from 
the NLO corrections of the 1/N expansion. It should be emphasized that u e ff 
depends itself on the distribution function via the retarded self-energy (I3.10ip . 
By writing down (16. 243p . we neglected off-shell processes included in (I3.106p . 
However, their contributions are expected to be small for the scaling behavior 
considered and we discard them in the following. 
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In principle, nonperturbative scaling phenomena may involve an anomalous 
scaling exponent for p(p°, p). Using isotropy we write 

p(p°, p) = s 2 ~ v p{s z p°, sp), (6.245) 

with a nonequilibrium “anomalous dimension” 77 . The dynamical scaling 
exponent z appears since only spatial momenta are related by rotational 
symmetry and frequencies may scale differently because of the presence of medium 
effects. Accordingly, the scaling behavior of the statistical correlation function 

F(p°, P) — s 2+Ks F(s z p °, sp) (6.246) 

will be described using a scaling exponent k s . This translates with the definition 
(13.89P for / p >l into 

/(p°,p) = s- + "f( S y,sp). (6.247) 

Using these definitions, one can infer the scaling behavior of u e ff(p 0 ,p)? which 
is given in terms of the “one-loop” retarded self-energy I1 r(i» 0 ,p) according to 
flihTOOl) . From fliUOB follows 

n R (p°,p) = s A U R (s z p°, sp) (6.248) 

with 

A = 4 — d — z + k s — 77 . (6.249) 

If A > 0 one finds from (13.1001) the infrared scaling behavior 

v eS (p°,p) = s~ 2A v eS (s z p°,sp). (6.250) 

(For A < 0 the effective coupling would become trivial with w e ff — 1, on which 
we comment below.) Employing these scaling properties, (I6.244P gives 


dU NLO (p, l, q,r) = s 


„ —2k s — 2 —2r) 


dU 


nlo ^sV, s z l°, s z q°, s 2 r°; sp, si, sq, sr) 


(6.251) 


Following the procedure of section 16.11 for any conserved quantity we can 
compute the flux through a momentum sphere k. Stationary turbulence solutions 
then require that the respective integral does not depend on k. Energy 
conservation, expressed in terms of the effective particle number distribution 
(13.941) . corresponds to the fact that 

<6 - 252) 
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is a constant of motion in this description. Also the effective particle number 
density 

n = / %Swf 2v ° Mr (6 ' 253) 

is constant for the collision integral (16. 2431) . Similar to (16. 231 1) . the flux for this 
effective particle number reads 

W = - 2W (, /2 + 1) /*> IPl ^- 1 CNLO (P)' 

The momentum integral can be evaluated along similar lines as before using the 
above scaling properties such that 

Ud—K, s +z—ri 

■4w~ d . K , +z _A NLO(1) - (6 - 255) 

Therefore, scale invariance may be obtained for the 

particle cascade: k s = d + z — rj (6.256) 


in the nonperturbative low-momentum regime. Similarly, for the scaling solution 
associated to energy conservation one finds, taking into account the additional 
power of p° in the integrand of (16.2521) . the exponent for the 

energy cascade: n s = d + 2 z — p . (6.257) 


With these solutions, we can now reconsider the above assumption that A > Oby 
plugging (I6.256P or (I6.257P into (I6.249p . Indeed, it is fulfilled under the sufficient 
condition that r] < 2. Taking into account that the anomalous dimension for 
scalar field theory is expected to be small for not too low dimension, with p 
maybe of the order of a few percent for d = 3, and employing a relativistic 2 ~ 1, 
we find k s ~ d + 1 for the particle cascade and k s ~ d + 2 for the energy cascade 
solution. 

The scaling of the effective occupation number distribution /(p), which 
depends only on spatial momentum, can finally be obtained from the definition 
(13.94p . This we write as 


/(p) + \ 



dp 0 


2 pV(p°,p) 


/(p°.p) + 5 



dp 0 

27 


2 p°F(p° : p), 


(6.258) 

using that ff° dp 0 / (2tt)p° p(p°, p) = 1/2 from the commutation relation (I2.98P 
in Fourier space. Of course, the scaling behavior (I6.246P for F(p°, p) may only 
be observed in a momentum regime with sufficiently high occupancies /(p) as 
implied in the above derivation. Then (I6.258P yields: 


/(P) = 


0 /c s +2— 2z 


f(s P> 


p| ^ +2 2 ^ relativistic particle cascade 
p|-(<2+2-?7) relativistic energy cascade 


(6.259) 
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These estimates show that vertex corrections can lead to a strongly modified 
infrared scaling behavior as compared to the perturbative treatment of 
section 16.1.11 

Anticipating that there is a mass gap in the relativistic theory, and because 
of the very interesting applications such as to the physics of ultracold atoms, we 
would like to compare this to the nonrelativistic case. The nonrelativistic limit 
is outlined in section (12. 8 p for the complex Gross-Pitaevskii field theory and we 
present here the relevant changes as compared to the relativistic case. To this 
end, we consider a nonrelativistic IV-component complex scalar field theory and 
perform again the 1/N expansion to NLO. Since already the relativistic scaling 
exponents (16.259j) indicate no explicit dependence on N at NLO - it can only 
enter indirectly via 77 and z - this seems to be a very good starting point to 
understand also the single complex field case of Gross-Pitaevskii. 

To proceed with the analytic estimate, we first note that for the nonrelativistic 
theory (I2.17ip with quartic (g/2)(x*x) 2 interaction the RHS of the evolution 
equation (12.176P for F (nr ^ has the same functional dependence on p (nr> and F 
as its relativistic counterpart. Therefore, proceeding in the same way as for the 
relativistic theory, with the corresponding scaling ansatz (I6.245P for p ( n d and 
(16.246ft for leads to the very same solutions (16.256ft and (16.2571) . A crucial 

difference arises when the occupation number distribution / nr (p) is determined. 
Using the nonrelativistic definition (II ,33ft for the distribution function, we have 
with the notation (12.174p in the absence of a condensate: 

1 f°° Hu 0 , 1 r 

/nr(p)+- = J' —F^ r \p°, p) = -J d 3 xe~ ipx ( X (t, x)x*(t, 0)+x(t, x)). 

(6.260) 

Comparison with the relativistic case shows that a difference in the scaling 
behavior is caused by the additional factor of ~ p° in the integrand of (I6.258p . 
Therefore, we find that 


/nr(p) = 


^k s +2—z 


/nr(sp)' 


p| ^ +2 ^ nonrel. particle cascade 
p|-(d+ 2 +z-r/) nonre i energy cascade 


(6.261) 


scales with one “z” difference than the relativistic solution (16.259ft . This can 
have important consequences, such as the fact that the scaling exponent for the 
nonrelativistic particle cascade is now independent of the dynamic exponent z 
describing the dispersion l>j p ~ |p| z . As a consequence, the same particle cascade 
scaling solution is found also in the presence of a condensate, where for the Gross- 
Pitaevskii theory the approximate (Bogoliubov) dispersion is given by 


^=V^(^ +29W2 )' (6 - 262) 

At large momenta, or in the absence of a condensate, one recovers u p = p 2 /(2m), 
while for low momenta one has u p ~ |p|. The particular importance of the 
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particle cascade and the associated phenomenon of Bose condensation far-from- 
equilibrium, as outlined in the introductory section [T] will be addressed in more 
detail next. 

6.2 Nonthermal fixed points 

6.2.1 Self-similarity 

In contrast to the stationary turbulence described above, we now turn to 
time evolution. In general, isolated systems out of equilibrium cannot realize 
stationary transport solutions. Instead, we have to consider the more general 
notion of a self-similar evolution, where the dynamics is described in terms of 
time-independent scaling functions and scaling exponents. 

A self-similar evolution of the distribution function f(t, p) for a spatially 
homogeneous and isotropic system is characterized as 

f(t,p) = s a/ Pf( S - 1/ Pt, S p) (6.263) 

with the real scaling exponents a and f3. Again, all quantities are considered to 
be dimensionless by use of some suitable momentum scale. Choosing s -1 ^ = 1 
we recover (11.251) . i.e. 

/(CP) =t a fs(t 0 p), (6.264) 

where the time-independent scaling function /s(t /3 p) = f(l,t^p) denotes the 
fixed point distribution. This scaling form represents an enormous reduction of 
the possible dependence of the dynamics on variations in time and momenta, since 
t~ a f(t, p) depends on the product £^|p| instead of separately depending on time 
and momenta. Therefore, an essential part of the time evolution is encoded in 
the momentum dependence of the fixed point distribution /s(p). Moreover, the 
values for a and f3 determine the rate and direction of transport processes, since 
a given characteristic momentum scale K(tf) = Ad evolves as K(t) = Ki(t/t\)~^ 
with amplitude f(t,K(t )) ~ t a . This aspect is also further discussed in the 
introductory section [I] 

For the self-similar distribution (I6.264p . the scaling behavior of a generic 
collision integral is then given by 

C[f](t, p) = s~^C[f](s~ 1/fi t,sp) = C[f s \{l,t p p) , (6.265) 

where fi is a function of scaling exponents similar to the discussion above. 
Substituting this scaling form into the kinetic equation leads to the time- 
independent fixed point equation for f s ( p), 

[a + /3p • V p ] f s ( p) = C[f s ]( 1, p) , (6.266) 

and the scaling relation 

a-l = -(3fi. (6.267) 
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This follows from comparing the LHS of the kinetic equation, 

§j:\ ta fs(t P P)] =t a ~ 1 [a + /3q- V q ]/ 5 (q)| q=t/3p , (6.268) 

to its RHS given by (16. 2651) . 

Further relations can be obtained by either imposing energy conservation or 
particle number conservation if applicable. For constant 

” = / f(t ' p> = t “~ w / (1^ /s(q) (e 269) 

one obtains the relation for 

particle conservation: a = /3d. (6.270) 

Similarly, one obtains from 

energy conservation: a = f3(d + z) . (6.271) 

One observes that there is no single scaling solution conserving both energy and 
particle number. As outlined already in the introductory section CD in this case a 
dual cascade is expected to emerge such that in a given inertial range of momenta 
only one conservation law governs the scaling behavior. 

Perturbative scaling behavior: We start again with a perturbative analysis 
of the relativistic theory, which is relevant at momenta above a possible mass gap 
such that cn p = s _1 o; S p, and afterwards we consider the highly nonlinear infrared 
regime. For the self-similar distribution (I6.264p . the scaling behavior of a generic 
collision integral for /-vertex scattering is then given by 

C(t, p) = s (i-iW0+*+i-(*-2)d C(s~ 1/f) t, .sp). (6.272) 

The reasoning is the same as for (16.2361) . with the only difference that now the 
distribution function scales according to (16. 2631) instead of (16.2321) . Consequently, 
the factor (/ — l)a/(3 appears, replacing the term (/ — 1 )k of (I6.236|) . Using 
the scaling relation (I6.267j) and (16. 270 j) from particle conservation gives the 
perturbative solution for 

relativistic particle transport: a = — 

Similarly, one finds the perturbative solution for 

, d +1 

relativistic energy transport: a = —— -- 


/3= -_ (6 . 27 3) 

< 6 ’ 274) 
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For instance, for quartic self-interactions the perturbative energy transport is 
characterized by a = —(d+ l)/7, and (5 = —1/7, where the latter is independent 
of the dimensionality of space d. Likewise, for the 3-vertex in the presence of a 
coherent field one has for the energy transport a = —(d + l)/5 and /3 = —1/5. 
The latter is used for d = 3 to describe the direct energy cascade of the inflaton 
dynamics in section [L3l 


Infrared regime: In order to go beyond the perturbative estimates, which are 
not applicable in the infrared, we again consider the 1/N expansion to NLO. 
At the lowest-order gradient expansion employed, the spectral function is time- 
independent for a homogeneous system according to (13.871) . Consequently, its 
scaling properties are still described by (16.245j) also for the self-similar evolution. 
However, for /(t,p°, p) we have to take into account the scaling with time: 

fit, P°, P) = s Qs//3 f(s- 1 /p t, s z p°, sp), (6.275) 

which increases the number of scaling exponents as compared to (I6.247P for the 
more restrictive case of stationary turbulence. With these definitions we can 
analyze the scaling of the collision integral (16. 243 ft and (16. 244 j) . which gives 

C NLO (t, p ) = S as/0 ~ z C NLO {s~ 1/f} t, sp) = t as ~P z C NLO (l, t p p) (6.276) 


following the corresponding steps done in section 16.1.21 Using the scaling relation 
(I6.267P we infer 


a = a s — (3z + 1 . 


(6.277) 


In particular, with the same reasoning as before we note that (I6.276P applies as 
well to the corresponding nonrelativistic collision integral. 

Again, the crucial difference emerges when the occupation number 
distributions f{t, p) for the relativistic and f nT (t, p) for the nonrelativistic theories 
are determined. For the former, we have according to (16.2581) 


f(t, P ) = s“ s//3+ 2-2z-r, 


2 (sV )p(s z p° , sp) /(s- 1 ^, s z p° , sp) 


= s“ s/ ^ +2 - 2z ”V(s _1//3 t,sp) = t a °+P( 2 ~ 2z -^ fs(tP p) 


(6.278) 


such that relativistically a = a s +/3{2—2z—ri) by comparison to (16.2641) . Together 
with (I6.277p one hnds in the absence of a mass gap for 

I f particles: a = d/{2 — z — p) 

relativistic transport: p =- of < 

2 — z — p [ energy: a = (d + z)/(2 — z — g) 

(6.279) 

Remarkably, the scaling exponent fd is obtained without using in addition energy 
or particle conservation, whereas the different solutions for a arise from imposing 
(I6.270P or (I6.27ip . respectively. 
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Similarly, the nonrelativistic case yields according to the definition 116.2601) : 


/ nr (t,p)=t“^ 2 -^)/ s , nr (^p). 


(6.280) 


Therefore, a = a s + /3(2 — z — rj) by comparing to (16. 264ft . which scales again with 
one “z” difference than in the relativistic case. Correspondingly, one finds for 


nonrel. transport: /3 = - of 

2 — rj 


particles: a = d/(2 — rf) (6 281) 

energy: a = (d + z)/{2 — rj) 


The nonrelativistic particle transport solution with taking 77 —> 0 is referred to in 
(I1.26P for the description of the inverse cascade in the nonequilibrium evolution 
of an ultracold Bose gas in section 11.41 


6.2.2 Inverse particle cascade and condensate formation 

In order to check the above analytic estimates, one may use the fact that 
the nonequilibrium quantum dynamics of the highly occupied system can be 
accurately mapped onto a classical-statistical field theory evolution as described 
in section [4] This mapping is valid for the system considered, as long as 
/ nr (t, p) S> 1 for typical momenta p according to (14.1621) . We employ initial 
conditions with high occupation numbers as motivated in section [T] where it is 
also described that the nonrelativistic and the relativistic theories show the same 
universal behavior for low momenta. 

In the following, we consider the nonrelativistic Bose gas described by (12.171ft 
in the dilute regime (( C 1). The initial distribution function at t — 0 is taken 

as 

/ nr ( 0 , p) ~ ^ @(Q - |p|), (6.282) 

which describes overoccupation up to the characteristic momentum Q. The initial 
condensate fraction is taken to be zero, i.e. |" 0 o 1 2 (^ = 0 ) = 0 in the defining 
equation ()1,33j) . with an initial / nr (0, p) = 50/(2 mgQ) 0(Q— |p|). We always plot 
dimensionless quantities obtained by the rescalings f nr (t, p) —> / nr (i, p) 2mgQ } 
t —$■ tQ 2 /(2m) and p —> p/Q- This reflects the classical-statistical nature of the 
dynamics in the highly occupied regime, which has the important consequence 
that if we measure time in units of 2 m/Q 2 and momentum in units of Q then the 
combination f m (t. p) 2mgQ does not depend on the values of m, g and Q. 

The initial mode occupancies (I6.282P are expected to get quickly redistributed 
at the beginning of the nonequilibrium evolution and then a slower behavior sets 
in. The latter reflects the dynamics near the nonthermal fixed point, where 
universality can be observed. We concentrate on the low-momentum part of the 
distribution and analyze its infrared scaling properties. 

Fig. fZT)l shows the results for the rescaled distribution (t/t ie {)~ a f nr (t, p) as a 
function of (f/f re f)^|p|, where the reference time t re {Q 2 /(2m ) = 300 after which 
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Rescaled momentum: F p 


Figure 20: Rescaled distribution function of the nonrelativistic theory as a 
function of the rescaled momentum for different times. The inset shows the 
original distribution without rescaling. 


self-similarity is well developed is the earliest time shown. In contrast, the inset 
gives the curves at different times together with the initial distribution without 
rescaling for comparison. With the appropriate choice of the infrared scaling 
exponents a and /3, all the curves at different times lie remarkably well on top 
of each other after rescaling. This is a striking manifestation of the self-similar 
dynamics (11.251) near the nonthermal fixed point. The numerical estimates for 
the scaling exponents obtained are 

a =1.66 ±0.12, 0 = 0.55 ± 0.03. (6.283) 

Comparing these values to the analytic estimates for the particle cascade in (11.261) 
or (16.28 1 [) . one observes that the numerical results (I6.283P agree rather well with 
the NLO approximation for a vanishing anomalous dimension 77 . Furthermore, 
the simulation results confirm that a = 3/5 to very good accuracy as expected for 
d = 3 from number conservation (16.2701) in the infrared scaling regime. 

The positive values for the exponents imply that particles are transported 
towards low momenta, which has important consequences. For the initial 
conditions (I6.282p . there is no condensate present at t — 0. However, the inverse 
particle cascade continuously populates the zero-mode, which will lead to the 
formation of a condensate far from equilibrium. 

Of course, the time needed to fill the entire volume with a single condensate 
increases with the volume V in which the simulations are done. Using that the 
parametrically slow power-law dynamics dominates the time until condensation 
is completed, we can use the scaling exponent a describing the growth in the 
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Figure 21: Evolution of the condensate fraction for the nonrelativistic Bose gas 
for different volumes V. The different curves become approximately volume 
independent after rescaling of time by V~ l t a , in agreement with (16.284)1 . 


infrared to estimate this condensation time. Taking the value of the zero- 
momentum correlator V -1 F( nr \t : t , p = 0) at the initial time t 0 of the self-similar 
regime as V^ 1 f nr (t o ,0) and its final value at the time tf as |'0 o| 2 (^/)j we can 
estimate from V~ l F <ynT \t 1 t, P — 0) t a the condensation time as 


t f ~ t 0 


( iVbl 2 (t f ) \ l/a 1/( 

\fnv(t O ,0)J 


(6.284) 


The time-dependent condensate fraction at zero momentum is given by 


Np(t) _ |^o| 2 {t) _ 

^totai /d 3 p/(2vr) 3 / nr (t,p) + |^o| 2 W ’ 


Accordingly, we may define the condensate fraction for the case of finite volumes 
as Nq/N total -t E -1 F (nr) (t,f,p = 0)/F (nr) (f, t, x = 0), using that N total = 
F^ nr '(t, t, x = 0). In Fig. ED we show the evolution of the condensate fraction 
for different volumes and rescale the time axis by Indeed, as predicted by 

(16.284)1 . the different curves are approximately volume independent. One finds 
that the condensate fraction saturates at A^/lVtotai — 0.8. 
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